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Scope of JET PROPULSION 


This Journal is a publication of the Ameri- 
can Rocket Society devoted to the advance. 
ment of the field of jet propulsion through 
the dissemination of original papers disclosing 
new knowledge or new developments. Ag 
used herein, the term “jet propulsion” 
embraces all engines that develop thrust by 
rearward discharge of a jet through a nozzle 
or duct; and thus it includes air-consuming 
—— and underwater systems as well as 
rockets. JET PROPULSION is open to con- 
tributions dealing not only with propulsion 
but with other aspects of jet-propelled flight, 
such as flight mechanics, guidance, telemeter- 
ing, and research instrumentation. In- 
creasing emphasis will be given to the scie tific 
problems of extraterrestrial flight. 
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Manuscripts must be as brief as the proper 
presentation of the ideas will allow. Ex- 
clusion of dispensable material and co):cise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract, and 12 illustrations, 
Fewer illustrations permit more text, anc vice 
versa. Greater length will be acceptable 
— in exceptional cases. 

hort manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes or Technical Comments. 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are usually published within two 
months of the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication according to the above-stated 
requirements as to subject, scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as & 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 

word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows: for journal articles: authors first, 
then title, journal, volume, year, page numbers; 
for books: authors first, then title, publisher, 
city, edition, and page or chapter numbers. 
Line drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under governt- 
ment sponsorship. Full responsibility rests 
with the author. 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
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are automatically considered for publication. 
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Dr. Hall, a research aerodynamicist at the Cornell Aero-_ 


nautical Laboratory, has had many years of experience in 
He was educated in Canada 
(B.A., 1950, University of British Columbia; M.S., 1951 and 
Ph.D., 1954, University of Toronto). After receiving his 
doctorate he was a research associate and assistant pro- 
fessor at the Institute of Aerophysics, University of To- 


shock tube investigations. 


ronto. 
this vear. 


He moved to the Cornell Aeronautical Laboratory 


Mr. Hertzberg, presently 
of Aerodynamic Research 


Assistant Department Head 
at the 
Laboratory, is well known for his work and publications on 


Cornell Aeronautical 
the development of the shock tube as a tool for various 
hasic Trained in aeronautical engineering at 
Virginia Polytechnic Institute (B.Ae.E., 1943) and Cornell 
(M.Ae.E., 1949), he has been associated with ( ornell \ero- 
nautical Laboratory most of his career. 


studies. 


Nomenclature 


( = specific heat capacity 

e = film voltage 

{ = film current 

k = heat conductivity 

K = thermal diffusivity, k/pc 
= film thickness 
= rate of heat transfer per unit area — 
= time 

T = film temperature 

T = average temperature through film 

R = film resistance —_ 

R; = film resistance at temperature 7; 

r = distance normal to film measured from exposed surface 


» = temperature coefficient of resistance 

6 = 

p = mass density 

i = VKt = characteristic thermal diffusion depth 
+ = variable of integration; time 

Subscripts 


F refers to film 


Brefers to backing material 


ETHODS of determining transient surface temperatures 

with a time resolution of the order of 1 microsee or 
better are of interest today in various fields of scientific and 
tigineering research. Until about 1953, the principal ad- 
vances toward reliable techniques for such measurement were 
made in connection with two problems: Determination of 
‘ylinder wall temperature variations in reciprocating en- 
sines, and determination of the surface temperature history 
of gun barrels subjected to continuous firing. In 1939, 
Meier in Germany (1),! following an earlier suggestion by 
Piriem (2), successfully used a 2-u thick gold film as a 


Re abel Sept. 19, 1958. 
' Numbers in parentheses indicate References at end of paper. 
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Recent Advances in Transient Surface Temperature 
| Thermometry 


J. GORDON HALL and A. HERTZBERG 


Cornell Aeronautical Laboratory, Inc., Buffalo, N. Y. 


shock tunnel research, 


surface resistance thermometer to record quasi-steady cylin- 


der wall temperatures. The gold film was evaporated in 
spiral form onto a thin insulating layer of aluminum oxide 
electrolytically deposited on an aluminum base plug. Around 
1941, Hackemann in Germany developed a rugged thermo- 
couple for gun barrel application (3). This consisted of a 
2-u thick nickel film, plated and ground on the end 
an encased steel wire. Imbedded in the gun barrel with the 
nickel film flush with the bore, the thermocouple had very 
rapid response and sensed a temperature close to the bore sur- 
face temperature. In subsequent years, American workers 
developed improved versions of this instrument (4 to 8). In 
a design described by Bendersky (6), the film thickness was 
reduced to ly thereby lowering the response time to 
+ Microsec. 

Around 1953, a new demand for a microsecond response 
instrument arose in the field of shock tube and hypersonic 
where it was necessary to determine 


quasi-steady gas flows of durations ranging from 1 millisec 


to 20 microsec. One approach taken here was to measure 


model surface temperature change with time and from this 
calculate the corresponding heat transfer. 
Bendersky-type thermocouple had sufficiently fast response 


Although the 


for the measuring times involved, it was considered too in- 
sensitive for general application. This difficulty was over- 
come by depositing thin films directly on dielectric models or 
inserts which experienced measurable surface temperature 
rise in the brief flow times available. Films of sub-micron 
thickness were tried in this way as both thermocouples and 
resistance thermometers. The latter method proved gener- 
ally superior and has been subsequently developed with a high 
degree of success. 

Another successful approach to shock tube heat transfer 
determination involved the development of relatively thick 
film surface calorimeters. These retain essentially all the 
heat absorbed during the short flow time and provide a direct 
indication of surface heat transfer rate. 

The purpose of the present article is to review the foregoing 
advances in transient surface temperature thermometry and 
calorimetry which have been made over a period of about 
the past five years. As the developments during this period 
have been made largely in the field of shock tube research, 
the discussion of the thin and thick film techniques whic ‘h 
follows makes frequent reference to shock tube heat transfer 
application. 


Thin Film Surface Thermometer Techniques 


The possibility of using very thin metallic films on dielec- 
tric bodies for transient surface temperature measurement in 
shock tube flows was apparently considered in several labo- 
ratories at about the same time. Early development of the 
techniques involved (1953-1954) was done at Princeton and 
Lehigh Universities (9,10) and at the Cornell Aeronautical 
Laboratory, Inc. (11). The early studies (e.g., (11)) achieved 
greater success using the surface film as a resistance ther- 
mometer rather than as an evaporated thermocouple (requir- 
ing two overlapping films). Erosion difficulties were experi- 


719 


‘S102 
E 
t 
of 
he 
one. 
ab 
one 
— 
nu- 
twi 
us 
a 
ted 
read 
e 
f 
or 
tof 
4 
¥ 
rst 
er 
et 
vi 
whit 
be 
tion. 
i 
1 
by 
ce i 
sified 
rests 
gi 
| 
i 
.00 
of 
the 
are 
ct the 
fon 
a 


enced with the latter, and, in addition, its sensitivity was 
~ relatively low. Accordingly, there was an initial concentra- 
tion of effort on the resistance thermometer, and its develop- 
- ment proceeded rapidly. While some success has been ob- 
tained with the evaporated thermocouple (12f), it appears 
to have been largely abandoned for shock tube use. In what 
follows, primary consideration is given to use of the thin film 
as a resistance thermometer, although much of what is said 

—_— with little or no change to the film thermocouple. 
The thermal response characteristics of the thin film can be 
determined from the classical theory of one-dimensional heat 
flow (13), and detailed analyses have been given by several 
authors (11,14,15). Film thicknesses used, of the order of 
0.1 yw, are ordinarily much smaller than the character- 
istic thermal diffusion depth 6 of the film metal for the shortest 
times of interest.2 Thus temperature gradients in the film 
may be neglected. The film senses the instantaneous sur- 
face temperature of the dielectric backing, but some response 
lag relative to zero film thickness occurs because of the small 
but finite film heat capacity. The magnitude of this response 
lag is illustrated by considering a suddenly applied heat trans- 
fer rate varying inversely as the square root of time. In the 
ideal case of zero film thickness / the surface temperature in- 
stantly jumps to a new constant value. For a finite film 
thickness, the backing surface temperature attains 94 per cent 

of the ideal value in a time t = 100 72, where 


2 (1,2 ]2 


is a characteristic time parameter of the film plus backing 
(15). With a 0.1 yw platinum film on quartz, for ex- 
ample, tr is about 3.5 X 1078 sec, and the response is ex- 
tremely fast. In the shock tube or shock tunnel, a suddenly 
applied constant rate of heat transfer is of interest. In this 
case, the surface temperature for zero film thickness increases 
as the square root of time. With a finite film thickness, 90 
per cent of the ideal temperature rise is attained in a time of 
100 vz. For testing times large compared to 100 72 the film 
heat capacity effects may be neglected. 

The technology of the preparation of thin films has been 
fairly well developed and reported in the literature. The 
methods of applying films to pyrex, glass or quartz surfaces 
include sputtering and evaporation under vacuum (16 to 
20), the painting technique using commercial metallic paints 
such as Hanovia preparations (11,21) and the use of metallic 
foils (10,21). The sputtering and evaporation methods 
provide good control of film thickness and uniformity on 
plane surfaces but are somewhat difficult to apply with 
sharply curved surfaces. After deposition of the metal on 
glass by sputtering or evaporation, it is common to bake the 
element at a temperature close or equal to the softening tem- 
perature of the backing (around 1100 to 1200 F for pyrex 
or quartz), after which it is allowed to cool slowly. The 
baking ensures some penetration of metal into the backing 
which provides a strong metal-to-glass bond and greatly in- 
creases the resistance of the film to abrasion. Electrical leads 
may be soft-soldered to the film with suitable flux. Plati- 
num, gold and rhodium have been the metals most used for 
sputtering or evaporation followed by baking. Reference 
(20) reports that evaporated platinum films on glass without 
baking are easily removed, whereas chromium, with its affinity 
for oxygen, which is believed important for bonding, forms a 
film very resistant to corrosion and abrasion. Similar results 
were observed in tests of (21) where, in addition, evaporated 
rhodium films (unbaked) were also found very durable. 

In the painting method the metallic solution is applied 
with a brush, pen or air brush. Then, the element is baked 
in a vented oven at a temperature near the softening point of 
the backing. During the baking process, the vehicle holding 


? The thermal diffusion depth 6 = ~/Krt is a measure of the 
depth of heat penetration in time ¢. 
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the metallic particles in suspension is driven off, and the par- 
ticles react with a reducing agent to produce a smooth film, 
As with sputtered or evaporated elements, the baking pro- 
duces a film which is very durable. An important advantage 
of the painting method is its simplicity, but film thickness 
and uniformity cannot be controlled to the degree possible 
with evaporation or sputtering. Provided the film thickness 
is small compared to the thermal diffusion depth 6, precise 
control or knowledge of thickness is unimportant. Uni- 
formity of thickness is desirable to provide uniform he: ting 
if the electrical pulse method of calibration is used. Refer- 
ence (21) reports that film resistance (thickness) is easil. re- 
peatable to +50 per cent, and temperature coefficients of film 
resistance are repeatable to +5 per cent using Hanovia } lati- 
num 05-X paint. 

The use of bonded metallic foils becomes attractive w here 
the film cannot be sputtered, evaporated or baked. Cold, 
copper and nickel foils are available commercially in ti:ick- 
nesses ranging upward from 0.1 micron. Gold foils ‘iave 
proved very satisfactory where erosion is not a problem 10), 
In shock tunnel tests with foils mounted at the stagn: tion 
point of a sphere, bonded nickel foil was found very durable 
and gold foil extremely fragile (21). 

Electrical stability of the film is necessary for its use as a 
resistance thermometer. The natural and artificial «ging 
of thin metal films has been extensively studied over the 
years (see for example (22)). With sputtered, evaporated 
and painted films baked at high temperatures as described, 
the films produced are usually quite stable. A slow increase 
of resistance with time is sometimes observed (e.g., (11)) 
which, while not a problem for short measuring times, is 
bothersome if it necessitates too frequent calibration. 

Because the physical properties of thin films are dependent 
on film thickness and can differ markedly from the bulk 
metal values (e.g., (17,20,38)), calibration of the thin film 
resistance thermometer is necessary for its quantitative use. 
Over the usual range of application, involving maximum 
temperature increases up to the order of a few hundred deg C, 
the film resistance-temperature relation can be assumed linear. 
The film resistance FR is then given by 


R=R,{1 + a(T — 7;)) 


where R; is the film resistance at initial temperature 7’;, and a 
is a constant coefficient which is determined by standard 
laboratory methods. Measured values of @ are always less 
than for the bulk metal, and typical values for sputtered or 
evaporated platinum and gold films might lie in the range 
0.0015 to 0.0025 per deg C. 

The principal aerodynamic application of the thin film as 
a resistance thermometer has been to heat transfer measure- 
ments in the shock tube and shock tunnel. Local surface 
heat transfer rates are calculated from the measured surface 
temperature histories by means of classical heat conduction 
theory. In the short flow times of the shock tube/tunnel, 
the depth of heat penetration into the film backing material 
is small, and the local heat flow may be assumed one-dimen- 
sional. Neglecting the response lag due to thermal capacity 
of the film, the instantaneous surface heat transfer rate ¢(t) is 
then given by (see for wr (18)) 


_ aT 
> q(t) = 
where B = V/ pakeCe is a constant which characterizes thé 


backing material. In terms of film voltage e Equation [I 
becomes 


. 


(t) B fi de 

q(t) = = | 

3In tests reported in (23) the linear resistance-temperatul 

relation was observed to hold to within +2 F from 70 F up ‘o the 

highest test temperature used, 850 F. 
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Thus, calculation of ¢g from the film voltage history requires 
the value of B/a. 

If there is some penetration of film particles into the back- 
ing material during the film baking process, values of 8 deter- 
mined from bulk properties of the backing may be unreli- 
able. Calibration is therefore desirable, and several methods 
for determining 6 under known heat flow conditions are 
used. A common technique (18,24) is to apply a sudden 
and known constant rate of heat input to the backing using 
a step function electric current through the film. The cur- 
rent is obtained with a long time constant, capacitor dis- 
charge. From the theoretical solution for this case (T= 
Vi) and an oscilloscope record of film voltage variation with 
time. the overall constant 8/a may be determined. For heat 
transfer measurement only, a separate determination of a@ 
is unecessary. The repeatability of this calibration tech- 
nique is reported to be +5 per cent (18,25). 

A second calibration procedure to determine 8/a@ involves 
flusl: mounting the element in a suitable aerodynamic model 
whic is then subjected to a known heat transfer rate by 
aero: lynamic heating in a shock tube flow (11). The known 
theoretical heat transfer rate and recorded film voltage again 
determine B/a. A repeatability of +5 per cent is reported 
also with this method (11). 

Experience has shown that with some glass backing mate- 
rials. such as fused quartz and pyrex, measured values of 8 
agree reasonably well with bulk values and do not vary 
greaily from one sample to the next. Consequently, for some 
applications assumption of an average value for 6 (about 
0.074 BTU/ft? F sec'/: for quartz and pyrex) may be. suffi- 
ciently accurate, and calibration can be avoided. In this 
case only the straightforward determination of @ is required 
toobtain heat transfer data. 

A preferable form of Equation [2] for numerical evaluation 
of heat transfer from oscillographs of film voltage variation 
with time is 


B et) ft elt) — er) 


which does not involve measurement of slopes. Reference 
(11) gives an expression for q which includes a first order cor- 
rection for film response lag where such is required. 

Shock tube and shock tunnel studies of hypersonic heat 
transfer distributions over models instrumented with thin 
films are currently being carried out by a number of research 
groups (19,25 to 33). Typical film geometry for these appli- 
tations is rectangular with dimensions of the order of 1 to 
2em by 1 to 2 mm. With a typical film resistance of 50 
ohms, constant current of 20 ma and resistant coefficient a of 
0,002 per deg C, a film output voltage of 2 mv per deg C o 
temperature rise is obtained. The heat transfer rates meas- 
wed in shock tunnels range roughly from one to several 
hundred BTU/ft?/seec. With flow durations of the order of 
| millisee, film temperature rise at the higher heating rates 
does not usually exceed 200 C. The overall accuracy of tem- 
perature rise measurement as estimated by various investi- 
gators at the present time varies from +5 to +10 per cent, 
ind that for heat transfer rate from +5 to +15 per cent. 
These limits indicate that heat transfer rates can now be 
teasured in the shock tube and shock tunnel with an accu- 
ey which possibly exceeds that attainable in conventional 
hypersonic tunnels. 

Another application of the thin film in shock tube research 
Sas a shock wave detector (9). The film is mounted flush 
it the shock tube wall and provides a very fast electrical 
wulse for trigger purposes on passage of the shock wave. 

One limit to quantitative accuracy of the film thermometer 
‘8 described above is imposed by the onset of thermal ioniz- 
ition of the gas at high temperatures (19,26,28,34). When 
oniztion occurs, it appears that the film resistance is par- 
ally short-circuited by the electrically conducting gas, and 
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the film voltage is no longer a known and simple function 
of the film temperature. Attempts have been made to 
eliminate the effects of ionization by evaporation of a thin 
insulating layer over the metallic film. In one instance, a 
1-4 silicon monoxide layer was evaporated and subse- 
quently baked to obtain a more durable and better insulating 
layer of silicon dioxide (19). Further experimentation with 
this technique is required before the usefulness of thin films 
in ionized flows is fully established. 

With its fast response, high sensitivity and ruggedness, it 
is to be expected that the thin film surface resistance ther- 
mometer will find useful application in fields other than shock 
tube and shock tunnel research. One example of this is its 
application as a total temperature transducer to study stag- 
nation temperature fluctuations in a jet engine inlet during 
flight (23). The transducer consists of a thin film (painted) 
on the spherical head of a pyrex probe. Stagnation point 
heat transfer rate is deduced from the measured surface tem- 
perature, and from the heat transfer rate the stagnation tem- 
perature can be calculated. The instrument is designed to 
sense total temperature fluctuations in the range —65 to 
1000 F lasting for 1 millisee. The probe operatcs continu- 
ously at 500 F and can withstand 1000 F for 4 min. 

Another example of thin film application is the hot film 
anemometer (35) which operates in a similar way to the well- 
known hot wire instrument. Reference (35) describes a film 
anemometer consisting of a pyrex probe with a 30 deg wedge 
tip on which a 0.2 X 1 mm platinum film is deposited. 
The film is heated by electric current to a temperature above 
the stream recovery temperature, and its resistance serves to 
indicate the rate of heat loss to the stream. Advantages 
claimed for the film anemometer over the conventional hot 
wire are increased signal/noise ratio, improved high-frequency 
response and excellent mechanical characteristics. 


Thick Film Surface Calorimeter Techniques 


Although the thin film resistance thermometer has a rather 
wide range of application in aerodynamic heat transfer meas- 
urement, its successful use becomes difficult under certain 
extreme conditions. These conditions include very high heat 
transfer rates and/or long flow times giving excessive surface 
temperature rise, high electrical interference such as may be 
present with are driven tunnels, flow ionization as mentioned 
above and excessive erosion by foreign particles in the gas 
stream. The shortcomings of the thin film in these respects 
may be reduced by increasing the film thickness. With the 
thin film, the thickness / is much smaller than the character- 
istic metal diffusion depth 6, and the film retains essentially 
no heat. If / is increased to a value equal to or greater than 
6, practically all the heat entering the exposed film surface 
is retained, and the film functions as a surface calorimeter 


(36). Practical upper and lower limits on / are set by the 
need for measurable film temperature rise and the need for 
minimum heat loss to the dielectric backing, respectively. 
6 is a practical compromise 


The 


Reference (36) shows that 1 = 
such that the heat loss may be neglected with small error. 
surface heat transfer rate may then be expressed as 


dT t) 


dt 


where 
1 l 
T(t) = t) dx 


is the average temperature through the film at time ¢. Hence, 
if the average film temperature 7 is measured as a function of 
time, the instantaneous heat transfer rate is given by the 
local slope of 7 vs. t. Although the mechanics of data 
reduction are thus simpler for the thick than for the thin 
film, the measurement of slopes is inherently inaccurate. It 
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might be noted that with / = 6, the thick film may still be 
thin by usual standards. With platinum, and a measuring 
time of 100 microsec for example, the value of 6 is about 50 yp. 

One method of measuring T is to use the film as a resistance 
thermometer. To a first order approximation, it may be 
shown (36) that the film resistance R is directly proportional 
to T, in which case q is given in terms of film voltage e by 


Second and higher order corrections to R due to the initially 
nonlinear temperature distribution through the film become 
negligible in a short time. A particular resistance calorim- 
eter for stagnation point heat transfer measurements described 
in (36) consists of a curved platinum strip (bulk metal) of 
dimensions 0.636 X 0.475 em soldered to a thin platinum film 
painted and baked on the nose of a pyrex model. With the 
strip thickness used, about 30 yw, the resistance is of 
the order of 4 X 107* ohms, and relatively large currents of 
the order of 5 to 10 amp are necessary for adequate voltage 
output. Bulk metal properties are used to determine the 
value of the constant e-Crl/aR; in Equation [3], and a di- 
rect calibration check is also made with electrical heating 
methods. The maximum calorimeter resistance change is 
limited to about 25 per cent to avoid possible significant 
error from variation of erC rl/aR; with temperature. 

This instrument has been used to measure stagnation point 
heat transfer rates in the shock tube ranging from 1000 to 
35,000 BTU/ft?/see with flow durations of the order of 25 
< to 50 microsec (33,36). The maximum surface temperature 

rise on pyrex with these conditions is around 2800 C, well 
above the melting point. The low resistance of the calorim- 
eter makes it less susceptible to ionization effects than the 
thin film, while its increased thickness gives it greater dur- 
ability. The latter is particularly important at model stag- 
nation points if foreign particles which cause erosion are 
present in the gas stream. 

Another type of surface calorimeter which has reportedly 
_ proved very successful under extreme conditions is the vari- 
_ able reluctance instrument (12e,37). This is used in are 
driven or “hot-shot’”’ tunnels where electrical interference 
is very high and is designed to measure heat transfer rates in 
the range of 10 to 1000 BTU/ft?/sec with flow durations up 
to 20 millisee. In this gage, a calorimeter sensing element 
‘0.18 in. in diameter and accurately flush with the surface 
varies the reluctance in a magnetic circuit. A voltage out- 
put proportional to the average temperature rise through 
- the element is obtained, and the slope of a voltage-time 
record gives heat transfer rate. The sensing element and 
_ inductance coils are housed in a cylinder of 0.25 in. diameter 
and 0.5 in. length. Calibration with a repeatability in the 
Bei _ range of 2 to 5 per cent is accomplished using alternating 

— hot and cool jets of water with an accurately known temper- 
: ature difference. The gage is reported to have a response 
time well under 1 millisec and an accurate reading duration 

of about 75 see (12e). The sensing surface is reported to 
withstand a small amount of melting and to normally endure 
20 to 50 tunnel runs with good repeatability (12e). It is 
ae noted that the low impedance of this calorimeter permits its 
use a few feet from a million amp are discharge. 
_ For some applications, a practical alternative to the surface 
_ calorimeter, where surface temperature rise is excessive for 
the thin film resistance thermometer, might be the Ben- 
 dersky-type nickel-steel thermocouple (6). This is a very 
rugged and fast instrument of sufficiently low thermal sensi- 
tivity that the (steel) surface temperature rise is only about 
one tenth that for glass. 


iv 


Conclusions 


The present review has covered some advances in transient 
surface temperature thermometry and calorimetry made 
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during the past five years in the field of shock tube and shock 
tunnel heat transfer research. Two basic techniques suc- 
cessfully evolved for shock tube heat transfer measurement 
are thin film surface thermometry and thick film surface 
calorimetry. The former method records instantaneous 
surface temperature from which instantaneous heat transfer 
rate is calculated. The latter method records total heat in- 
put at the surface, and instantaneous heat transfer rato is 
given by the corresponding time rate-of change. 

The most successful type of thin film instrument is pres- 
ently the thin film resistance thermometer, which is simple, 
sensitive, rugged, and which has microsecond response. In 
present applications, where the surface temperature does not 
exceed about 200 C, the film resistance thermometer is « ap- 
able of 5 per cent accuracy in measurement of surface tem- 
perature and heat transfer rate. Extension of its quantitative 
use to ionized flows is a matter of current research. Its de- 
sirable characteristics have led to its application in fi-lds 
other than shock tube research. 

The surface calorimeter has inherent low sensitivity :ind 
electrical impedance, and is advantageous under condit:ons 
of extreme surface temperature rise or high electrical inter- 
ference. Its increased thickness makes it more durable than 
the thin film, and its low impedance renders it less susceptible 
to flow ionization effects. Both resistance and variable re- 
luctance surface calorimeters have been successfully devel- 
oped for extreme conditions where use of the thin film resist- 
ance thermometer would not be feasible. 
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A program for the unmanned reconnaissance of the solar 
system at a moderate cost, in the near future, is suggested. 
{nonreturn, nuclear powered ion rocket propelled recon- 
naissance vehicle weighing 3300 lb is described. A vehicle 
weight statement, drawings and a description of a possible 
mission trajectory are included. 
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x = accelerating gap distance, em a 
= partial differential symbol 
£ = angle from final take-off point to Earth’s heliocentric 
velocity 
6 = angle from final take-off point to vehicle 
@ = angular acceleration of vehicle 
6 = angular velocity of vehicle 
6 = angle from hyperbola axis to cut-off point 
-@ = angle from hyperbola axis to asymptote 
g = angle from hyperbola axis to cut-off velocity 


fase 
HE primary aim of space travel in the near future is the 
acquisition of knowledge. Obviously, it is desirable to 
obtain this information at a minimum expenditure of both 
- money and manpower. The acceptance of this requirement 
~ leads logically to a system for utilizing robot, instrumented, 
- nonreturn space vehicles. Project Snooper is such a system. 
_ Snooper reconnaissance vehicles (Fig. 1) operating in satellite 
orbits about other planets of the solar system, would gather 
and transmit information about these planets to Earth. This 
can be accomplished using current electronic technology. The 
inclusion of humans in such vehicles, although possibly allow- 
ing the collection of more information, would greatly add to 
the complexity and cost of the system to meet the require- 
- ments for crew provision and return to Earth. 

The main task in the development of the Snooper recon- 
__-naissance system is then reduced to achieving a low weight, 
 interorbital space vehicle which can be placed in an Earth 
satellite orbit with a feasible conventional rocket booster; 
_ perhaps similar to the type necessary for an intercontinental 
ballistic missile. To satisfy this requirement, as well as that 
_ for carrying a reasonable size reconnaissance payload, a pro- 
_ pulsion system of specific impulse one hundred times greater 
than that available with conventionally powered rockets is 
required. Only the photon and ion rocket principles provide 
sufficiently high specific impulse. The photon rocket must 
be removed from consideration because its temperature and 
power requirements cannot be satisfied using our present 
_ technical knowledge. The ion rocket alone is suitable for 

_ propelling the robot Snooper vehicle. 

_ An ion rocket is a device for obtaining thrust through the 
7 high velocity expulsion of electrically accelerated ionized 

particles. In the case of Snooper, the source of electrical 
; np is a closed system nuclear powered turbogenerator. 
; The powerplant would provide most of its electrical power 
- output to the propulsion system during powered flight. All 
of its power would be available to the payload when the pro- 
- pulsion system is not operating, so that the powerplant could 
be considered part of the useful payload. 
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Fig. 2 Ion rocket power supply system 


Once developed, fleets of Snooper vehicles would have the 
capability of gathering the information necessary botl. to 
motivate further space exploration and to achieve practical 
manned space travel at minimum cost, maximum expediency. 


Nuclear Powerplant 


Power Source 


A nuclear powerplant was selected for the reconnaiss:unce 
vehicle after a comparison with a solar energy powerplant. 
The difficulties associated with packaging the solar mirrors of 
a solar energy powerplant for the boost period, the varying 
amount of solar energy received as a function of distance from 
the sun, and the inability to operate in shadow excluded fur- 
ther consideration of a solar energy powerplant. 

The low weight requirement of the nuclear powerplant for 
this application appears possible if radiation shielding is omit- 
ted. This is permissible in the Snooper vehicle because it is 
unmanned. On the basis of published data (1)* on current 
reactor design and performance, the estimated weight of the 
unshielded Snooper vehicle reactor appears to be reasonable. 

For the system assumed in this paper, a sodium coolant is 
used. There are, of course, other possible reactor coolants 
depending on the exact reactor type used (2,3). 

The Snooper propulsion system requires that the thermo- 
dynamic system have an output sufficient to produce 147 kw 
of electrical power. To do this the reactor must have a ther- 
mal output of approximately 1 megawatt. This requirement 
for reactor power, along with the total probable thrust dura- 
tion, is a fortunate circumstance, for it appears that there 
would be sufficiently low consumption of fissionable material 
in the reactor to allow long term operation without further 
addition of fuel. This should allow considerable simplification 
in the reactor design. 


Thermodynamic System 


In the system illustrated in Fig. 2, the sodium reactor cool- 
ant is passed through a heat exchanger to boil the turbine work- 
ing fluid, in a single-fluid Rankine cycle. Mercury was se 
lected as the turbine working fluid because it allows the system 
to operate at a moderate pressure at high temperature, i? 
comparison for example, with steam. 


Heat Rejection System 


The heat rejection system of this powerplant may be col- 
sidered one. of the major problem areas in the design of ion pre 
pulsion systems. Each heat generating component in the en- 
tire vehicle system must eventually have its heat rejected 
from the vehicle by radiation, the only means for rejecting 
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heat in the vacuum of space without mass ejection. Radia- 
tor sections would be constructed of thin gage aluminum sheet 
(0.010 in. or thinner if feasible) upon which are brazed, at ap- 
proximately 2-in. intervals, }-in. diam thin walled copper 
tubes for coolant flow. These tubes run parallel to the longi- 
ke tudinal axis of the vehicle, thereby requiring no bending of 
the tubes during radiator stowage (Fig. 3). The manifolds 
Re would be made of inflatable copper sections, which would be 
expinded (and remain expanded) when the radiator coolant 
fluid is initially pumped into the radiator, after it is un- 
' wrapped in the Earth satellite orbit. The entire exposed sur- 
face of the radiator would be suitably coated, perhaps with 
ey soot, to provide a surface having a high emissivity. 

Sodium was selected for the radiator coolant fluid because 
~ of its excellent heat transfer properties and low density. 
Upon comparison with mercury, the weight advantage of 
sodium was found to be quite pronounced, even though a 
mercury system can operate at approximately constant tem- 
periture by using the radiator as a condenser. 

The problem of the relative position of the radiator sections 
wit!) respect to incident rays from the sun was considered. 
re the It appears that no radiator hinging or gimballing mechanism 
tl: to is required because the Snooper vehicle would remain essen- 


tially in the plane of the ecliptic during its flight. 


Electrical Power System 7% 


The electrical power system would include an a.c. generator, 
a rectifier, and a control and distribution system. Both d.c. 
ssince — and a.c. generators were considered for this application. The 
plant. | dc. generator was rejected on the basis of the low power-to- 
rors of § weight ratio presently obtainable from such machinery. 
arying With an a.c. generator, it appears that a power-to-weight ratio 
efrom § of approximately 1 kw per lb can be attained (4). A judicious 
d fur- § arrangement of both the electrical and mechanical rotating 
machinery would be necessary to prevent precessional mo- 
int for § ments resulting from gyroscopic action. 


omit- 

of the 

nable Principles of Operation 

nak ® An ion rocket produces thrust by the reaction of a stream of 


olants F electrically accelerated ions. The acceleration of a mass to a 
high velocity requires power. The power in the ion exhaust 
1erm0- 7 beam is expressed by (5) 


47 kw 

me This beam power, although less than the power which must be 
“there § Provided by the ion rocket power source, is a measure of the 
aterial specific power requirements of an ion rocket propulsion system, 
varthet A general plot of the ratio of required ion beam power to ve- 


hicle weight (6) for a range of vehicle accelerations and ex- 
haust velocities is shown in Fig. 4. Applied to ion rockets, 
this curve provides a means of determining either the maxi- 
mum vehicle acceleration or the ion exhaust velocity for a 
given ratio of ion beam power to vehicle weight. 


ication 


yr 

> work- Propellant Selection 

system One of the primary requirements of an ion rocket propellant 

ure, in § 8 that it have favorable ionization characteristics as nonion- 


ized atoms contribute no thrust to the vehicle. 

Of all the elements, the alkali metals (lithium, sodium, 
potassium, rubidium and cesium) have the lowest ionization 
potentials. Because of this, the alkali atoms may be ionized 
be col § on incandescent metallic surfaces, such as platinum or tung- 
on pro § sten (7), two metals whose respective work functions (6.2 and 
the el § 45 ev) are higher than the ionization potentials of the alkali 


eject tlements. When an alkali atom contacts an incandescent 
je cting § platinum or tungsten surface it becomes ionized after a stay 
F ppapet time of only several microseconds; the probability of an al- 
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Fig. 3 Method of radiator stowage 
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Fig. 4 Specific power requirements for rocket propelled vehicles 


kali atom being ionized when contacting such a hot surface is 
essentially unity. 

Of the alkali metals, cesium and rubidium share the desir- 
able properties of low ionization energy, low heats of fusion 
and vaporization at relatively low temperatures, and high den- 
sity. 

Another criterion for propellant select‘on is the required 
area of the ion source. The accelerating voltage in an ion 
motor is expressed by (5) 


komc? 
[2] 
From Equations [1 and 2], the ion beam current is 
ks ma 
1c 


For a desired value of exhaust velocity (specific impulse) and 
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thrust the voltage and current can be related to the atomic 
weight as shown in Fig. 5. 
The ion current density, based on space charge considera- 
tions, is expressed by (7) 


THRUST = 0.5 LB 
EXHAUST VELOCITY = 657,000 FT/SEC Vie 


Using Equation [4] and Fig. 5, the required ion source area is 
given in Fig. 6 for several of the alkali metals. The ratio of 
the ion source areas of cesium and rubidium from Fig. 6 shows 
the ion source area for cesium to be approximately 40 per vent 
that of rubidium for gap distances of 1to3cm. Thus cesium 
appears to be the most desirable propellant for the Snooper 
design. 
Propellant Feed System 

ATOMIC WEIGHT vars pe ys 


_ _ Fig. 5 Typical ion rocket electrical requirements One possible propellant feed system suitable for the Snooper 
vehicle is shown schematically in Fig. 7. The cesium 
propellant is stored in the liquid state. Its temperature is 
maintained at approximately 100 F by an electrically heated 
blanket surrounding the propellant tank. (The melting point 
of cesium is 83 F.) 

The necessity for operation at exceedingly low vehicle ac- 
celeration (approximately 0.1 milligee) indicates the need of a 
positive expulsion propellant feed system for the Snooper 
vehicle. This is accomplished by an expandable metal bag, 
which is slowly inflated by pressurized gas as the propellant is 
consumed. The liquid cesium passes through a check valve 
and is then sprayed onto the wall of a vaporizing chamber by 
means of an atomizer. The vaporizing chamber wall is main- 
tained ata temperature of 1400 to 1500 F, providing a surface on 
which the atomized liquid cesium is vaporized. (The boiling § * 
point of cesium is 1301 F.) The cesium vapor then flows un- 
der vapor pressure through a sintered metal restrictor which 
functions as a precalibrated flow metering system. Twoon-off * 
control valves are in the system, immediately downstream of 
the propellant tank and at the vaporizer respectively. Dur- 
ing the startup of the propulsion system, the opening of the 
on-off valve downstream of the vaporizing chamber is delayed 
until the design operating vapor pressure is achieved in the 
vaporizing chamber by introduction of the propellant. 

One method of thrust level control which appears possible 
with this feed system is propellant flow control by varying the 
reference pressure setting on the pressurizing gas contro] valve 
and the wall temperature of the vaporizing chamber (and 
hence the cesium vapor pressure). The variation in the 
THRUST = 0.5L cesium vapor pressure will change the pressure drop across the TI 


EXHAUST VELOCITY = 657000 FT. . . 
restrictor and hence the cesium flow. 
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Two ion motors, each of $-lb thrust and an ion current 
rating of 2.67 amp are used in the Snooper design. These 
motors would provide an initial acceleration of 0.1 milligee to 
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Fig. 6 Typical ion source area requirements 


ge tari nitarhags The motor design is based on one proposed by Dr. E. Stuhl- the 
inger of the Army Ballistic Missile Agency, Huntsville, Ala. 
(7). Cesium propellant, delivered from the propellant 

ON-OFF wives ee, Se, : feed system, is fed into the motor as a vapor (Fig. 8). After rect 


LIQUID CESiuM AT 100 F 
/ passing through a distributor, the cesium impinges on incandes- ale 


cent tungsten surfaces, which are formed into a series of ion- traj 
izing grids. The cesium is ionized at these grids and thenac- § j, 
celerated through an acceleration gap to a velocity of 657,000 nlac 
ft per sec by a d.c. potential of 27,500 v. The accelerating or § 
_ cathode grid has a honeycomb cross-section to allow a nearly The 

uniform electrostatic field in the radial direction. lise 

If only electrically positive particles were ejected, the vehi- inst 
cle would soon attain a large negative charge which would the j 
hinder and even prevent the further expulsion of ions. This tient 
effect is cireumvented in the Snooper propulsion system in tive 
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Fig. 7 Schematic diagram of propellant feed system 
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much the same manner as that proposed by Dr. Stuhlinger— 
by the simultaneous emission of electrons which would mix 
with the ejected ions forming an electrically neutral plasma. 
As shown in Fig. 8, a thermionic electron emitting surface 
along with its associated heating element is mounted on the 
downstream side of the cathode grid. 


Payload 


The payload of the Snooper vehicle includes television, 
radar, communication equipment and auxiliary power supply 
systems. Devices for measuring cosmic ray and electromag- 
netic field intensity about the planets as well as solar effects 
at various distances from the sun could be included. This 
section would also include provisions in its design for environ- 
mental control where necessary, and would be equipped with 
hydraulie and electrical systems for the actuation and control 
of antennae, cameras and transducer devices. Much of this 
equipment would be specially designed to allow operation in a 
radiation field. 

The payload would be designed to accommodate a variety 
of tasks. Only one of these is discussed here to give an indi- 
cation of the information gathering potential of the Snooper 
recounaissance system. This task is the televising and radar 
scanning of the surfaces of planets and the transmission of 
this information to Earth. A general block diagram of the 
scanning, transmission and receiving system is presented in 
Fig. 9. Television and televised radar signals would be put 
on magnetic tape (continuously erased and revised). Follow- 
ing this taping process, the information is coded in order to 
transmit it in a narrow bandwidth signal. This signal would 
be received by monitoring stations on the Earth, where es- 
sentially the reverse process would take place. This system 
would reduce significantly the power required for the transmis- 
sion of information to Earth. Information theory can prob- 
ably be of great aid in further reducing power requirements 
and improving the transmission system. 

More advanced payload equipment might include small 
aerodynamic test vehicles (weighing perhaps 100 Ib each) 
which would be decelerated initially out of the observation 
orbit by a small rocket. The descent of such vehicles, in- 
strumented and carrying telemetering equipment, would allow 
the collection of considerable atmospheric data. The use of 
this telemetered information could conceivably lead to the 
placing of small instrumentation packages on the surfaces of 
the planets. This package could even carry a television cam- 
era-transmitter, possibly like the lightweight portable televi- 
sion camera developed for the U. S. Army, ‘“Creepy-Peepy.”’ 
The surface exploration package could even be self-propelled 

and thus deliver highly detailed information about a relatively 
large area of planetary surface. The possibility of such ad- 
vanced instrumentation systems indicates some of the recon- 
taissance potential of Project Snooper. 


> 
Trajectory Analysis 

In showing the feasibility of the Snooper propulsion system, 
the ability to move along a reasonable trajectory must be 
proved. To do this, a trajectory must be established and ex- 
amined to see if travel times are reasonable, if navigation cor- 
reetions can be made, ete. Snooper has been designed for 
wlar system reconnaissance. However, to set up a reference 
trajectory, Mars will be used as a sample objective (8). 

It shall be assumed here that Snooper has previously been 
jlaeed into a close satellite orbit, approximately 300 miles 
tbove the surface of Earth by means of chemical fuel rockets. 
The trajectory of the Snooper vehicle can be divided into four 
liserete parts, starting from this Earth satellite orbit. The 
irst part of the trajectory is powered flight, sustained until 
he ion rocket reaches some arbitrary orbit about Earth suffi- 
tiently far from the Earth’s center of gravity so that an exces- 
‘ve amount of energy need not be expended during the second 
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Fig. 9 Block diagram of TV-radar reconnaissance system 


part of the trajectory. The second part of the trajectory con- 
sists of resumed powered flight sustained until the necessary 
velocity is obtained to escape the Earth’s gravitational field 
and place the vehicle into the elliptical orbit connecting the 
solar orbits of Earth and Mars. The third part is the ellipti- 
cal orbit along which the rocket ‘“‘coasts.”” As the ship nears 
Mars, powered flight is resumed to maneuver the vehicle into 
a satellite orbit, about Mars, the fourth and last part of the 
trajectory. 

Only the theory describing the first three parts of the tra- 
jectory is presented in this report. Approximations have 
been made for the purpose of simplifying this analysis. Cir- 
cular, coplanar solar orbits are assumed for Earth and Mars. 
The perturbing effects of the moons of the two planets are not 
considered. During powered flight, solar effects are ignored 
and during the elliptical coasting orbit, the gravitational ef- 
fects of Earth and Mars are neglected. 

To derive equations of motion for the first sections of the 
trajectory the Lagrangian equation of energy is applied (9) 


d foT 
dt (37) ogi 


are the generalized coordinates, the Q’s the 
The equations of 


where the q’s 
generalized forces, and 7’ is the energy. 
motion are reducible (8) to (Fig. 10) 


The solution of the equations of motion must be done by nu- 
merical means involving a great deal of computational work. 
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The initial conditions are imposed by the orbit of departure 
and the angle of applied thrust. 

The vehicle is expected to spiral around the Earth, slowly 
increasing its energy, moving farther and farther from its orig- 
inal orbit. At some predetermined point the propulsion 
system is turned off, and the vehicle coasts into a new satellite 
orbit relatively close to the entry point of the circum-solar 
ellipse. While the ion rocket is in this new geocentric orbit, 
its exact position and velocity can be observed from Earth. 
These observations would then provide new initial conditions 
from which the interplanetary trajectory can be computed. 
Information concerning the position at take-off, the thrust 
duration (to determine cut-off point) and any other pertinent 
information can then be relayed to the vehicle. The equa- 
tions of motion are, as before, with the initial conditions a func- 
tion of the new orbit and desired angle of thrust application. 

When the observations and computations have been com- 
pleted, the vehicle resumes its spiral around the Earth until it 
builds up enough energy to escape the Earth’s gravitational 
field, and then proceeds into a circum-solar ellipse which has 
its aphelion in proximity to the solar orbit of Mars (Fig. 11). 
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The cut-off point will be timed so that as the vehicle follows 
the escape hyperbola, the escape asymptote of the hyperbola is 
parallel to the Earth’s circum-solar orbit and in the same di- 
rection. This escape hyperbola leads directly into the «le- 
sired circum-solar ellipse (Fig. 12). 

The point of take-off for the second part of powered flight is 
determined by the requirement that the escape asymptote be 
parallel to the orbital motion of the Earth 


> & = 6. + 90 deg — 6B + 90 deg — a - 
—a-+t 180 deg 


The above angles are defined in Fig. 12. 

The problem of errors occurring during the various portions 
of the trajectory is serious. These errors must be correcte: in 
order to ensure the successful completion of any mission. ‘lhe 
real problem then, lies in the ability of the vehicle to mike 
corrections for these errors. 

In the theory as discussed here only one attempt has ben 
made to compensate for errors. This is done by separating 
the Earth escape trajectory into two portions. Those errors 
incurred during spiraling from the first orbit about Earth: to 
the second can be compensated for by accurately determining 
the initial conditions for the second power flight and then pro- 
ceeding with the trajectory calculations and corrections. 


Guidance 


The navigation of the Snooper vehicle would be accom- 
plished with an automatic solar-planetary navigation system. 
Star triangulation does not appear to be usable because of the 
exceedingly small parallax of the stars. In turn, however, the 
stars would provide an extremely accurate reference system. 
A solar-planetary navigation system will enable the vehicle’s 
position to be determined at any instant by observing the 
angles between the planets and the sun. 

A highly accurate chronometer would be required based, 
perhaps, on the exact process of radioactive decay. A con- 
tinuously available frame of reference such as a gyrostabilized 
platform will probably be required. In the environment of 
total or near weightlessness, the drift of a gyro may be greatly 
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reduced because of the reduction of significant torque forces. 
The spatial orientation of the stable platform would be moni- 
tored by an automatic star tracking device to maintain the 
reference axes in the proper spatial direction. The perfect 
seeing available to any optical instrument in the vehicle’s 
guidance system, operating in the vacuum of space, is a defi- 
nite advantage in the development of such a system. 

Prior to any mission, an elaborate ephemeris of the voyage 
must be compiled. This information would be fed into the 
thrust level and vector control systems to achieve the desired 
traiectory in space. 


Vehicle Configuration and Operation 


The general configuration of the interorbital vehicle is 
shown in Fig. 1. The payload section of the vehicle, contain- 
ing the guidance and payload instrumentation and equipment, 
would be separated from the nuclear powerplant-ion propul- 
sion section during operation by a telescoping radiation isola- 
tion structure. The physical separation of the payload-guid- 
ance section from the powerplant nuclear reactor reduces the 
effects of harmful incident nuclear radiations on the payload 
and guidance instrumentation. The large, winglike radiators 
are extended on either side of the vehicle. 

The powerplant and propulsion systems of the Snooper 
form a compact package in the aft section of the vehicle. The 
nuclear reactor, turbogenerator, propellant storage tank and 
two ion rocket motors are located here. The ion rocket 
motors are gimbal mounted outboard of the reactor, and can 
be rotated to allow deceleration without changing the attitude 
of the vehicle in space. A preliminary weight estimate of 
this configuration is presented in Table 1. 


Table 1 Preliminary weight statement, Snooper 
reconnaissance vehicle, Ib 

Payload and guidance... 1500 
Reactor, heat exchanger, pumps and turbine 1000 
Radiator. 280 
Generator and electrical system............... 150 
Structure and miscellaneous.................. 150 


The Snooper vehicle would be accelerated into an Earth satel- 
lite orbit as the payload of a chemical propellant rocket 
booster similar to that necessary for a possible intercontinental 
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The problem of calculating the linearized motion of an 
incompressible inviscid fluid with a free surface is formu- 
lated as a variational principle. It is demonstrated that 
the variational principle is of the free boundary condition 
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ballistic missile. During the Earth satellite boost portion of 
the flight, the radiators are stowed overlay wrapped (Fig. 3), 
the ion motor mount arms folded, and the radiation isolation 
structure telescoped to form a compact assembly. A jetti- 
sonable bracing system and windshield are utilized to protect 
the Snooper during its acceleration into an Earth satellite 
orbit. Once in the orbit, the windshield bracing system is 
jettisoned along with all mass not a functioning part of the 
vehicle. Then the radiator sections are unwrapped; the 
radiation isolation structure is extended, and the ion rocket 
motor arms are unfolded. 

The nuclear power source is now activated and the low 
thrust ion rocket propulsion system accelerates the Snooper 
along its trajectory. When the desired orbit about its destina- 
tion planet has been achieved, the power of the reactor is uti- 
lized to operate the automatic data collection and transmission 
equipment. The apparatus would operate for approximately 
one year after it has reached its destination orbit, thus allow- 
ing sufficient time for the collection of useful data. 
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type. A procedure for calculating the natural frequencies | 


and modes of a shallow and a deep tank by applying the 


Rayleigh-Ritz procedure to solve the variational formula- 
tion is discussed. An approximate procedure for calculat- 
ing the sloshing modes and frequencies of a tank of inter- 
mediate depth in terms of the solutions for shallow and 


deep tanks is suggested. The highly satisfactory conver- 
gence and accuracy of the method are demonstrated by 
comparison with exact methods in the cases of flat bot- 
tomed tanks and 90 deg conical bottomed tanks. The 
results of numerical calculations for conical tanks are 
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— presented in the form of mode shape diagrams and fre- 


quency charts. 


Nomenclature 


= constants to be determined 
= defined by Equation [25] 
= y for shallow tank approximation mi 
= £, when the latter is function of r only ss a 
= gravitational or total vertical acceleration = 


= Bessel function of order s of the first kind 

= Lagrangian = T — V 
N = degree of thé polynomial 
n = normal to boundary surface 
r, 2,@ = cylindrical coordinates 

S = free surface of liquid 
number of mode in 6 
kinetic energy of the liquid 
= time 

velocity of fluid 

= potential energy of the liquid 
= rectangular coordinates 

= polynomial coefficients 

= dimensionless frequency defined in Equation [35] 
= phase angle 

= zon the fluid surface 

= frequency parameter 

= function of r and z; = &(r, z) cos s @ (or sin s 
mass density of fluid - 
= interface of liquid and tank . 

= circular frequency of time function, rad/sec 

= volume of the liquid 

= potential function dependent on r, z, 0(or x, y, z) and / 
= space part of ¢; @ = ¥(r, z, 0; or, y, 2) sin (ot + e) 
= potential energy per unit mass 
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Introduction 


bi HE flow field of the liquid contained in a tank with a 
if free surface when subjected to a transient acceleration 
field has been investigated since the middle of the eighteenth 
century. The basic partial differential equations and the 
classical treatment by linearization and separation of variables 
were essentially complete by 1895, as reported in Lamb’s 
first edition of ““Hydrodynamics.”’ Renewed interest in this 
subject has been generated recently by the engineering prob- 
lems arising in the design of liquid rocket powered missiles. 
In particular, the design of missile autopilots requires an 
accurate estimate of the natural frequencies of the sloshing 
modes which lie within the autopilot effective control fre- 
quency band during powered flight. In the development of 
propellant utilization systems which continually measure the 
residual propellants remaining in the tanks and adjust the 
engine mixture ratio accordingly, a knowledge of the dynamic 
behavior of the propellants is necessary. 

Linearization of the equations of motion and the imposition 
of the free surface boundary condition on a plane are employed 
in the classical theory. An order analysis to establish the 
range of validity of these assumptions is required, but will 
not be attempted here. The classical linearization will be 
employed with the understanding that the range of validity 
is unknown from the analytic viewpoint. In practice, ex- 
perimental data on mode shapes and frequencies seem to 
indicate that the linearization assumptions lead to satis- 
factory results for a wide class of problems. 

Additional assumptions used in the classical theory and 
adopted here are: The viscosity and compressibility of the 
fluid may be neglected, the fluid starts from rest and the forces 
acting on the fluid may be derived from a potential function. 
An important consequence of these assumptions is that the 
flow field of the fluid may be calculated in a relatively straight- 
forward manner if the natural frequencies and sloshing modes 
are known. 

Methods for « 
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calculating the sloshing modes have been de- 


= height of surface above bottom a positive function of r 


veloped in the classical theory for tanks of circular, elliptical 
and rectangular cross sections with flat bottom. These 
solutions make use of the technique of separation of variables 
in the differential equation. It can be shown that the appli- 
cation of this basie technique cannot lead to further useful 
solutions. 

If the tank is very shallow, the problem becomes two- 
dimensional regardless of the bottom shape. An exact solu- 
tion has been developed for a circular tank with a botton. in 
the form of a parabola (1).4 If the tank is infinitely deep, the 
bottom can always be considered as flat, and the equations 
also reduce to a two-dimensional problem. Such probl« ms 
may be solved for arbitrary cross sections by a variety of 
methods. The tank with a nearly flat bottom may be treated 
analytically by a perturbation analysis (1). Such cases, 
however, are of minor interest in technical applications. 

To develop approximate solutions to the sloshing moiles 
and frequencies for tanks of arbitrary depth and bottom 
shape requires the development of a technique which is busi- 
cally different from the classical method of separation of v..1i- 
ables. From an examination of similar mathematical prob- 
lems in other areas of physics, the use of integral equations 
or variational methods appears to offer promise. In this pa) er, 
the application of variational techniques is developed «nd 
exploited to produce a general procedure for approximating 
sloshing modes and frequencies. The method of intezral 
equations has not been studied but is believed to be equ:lly 
applicable and should be pursued. 

The variational method: developed here for the three- 
dimensional problem is solved approximately in terms of the 
solution of a pair of two-dimensional problems: The very 
shallow tank and the very deep tank. A variational pro- 
cedure for solving these special cases is presented. These re- 
sults are believed to be useful independent of their applica- 
tions to the three-dimensional problem. 

It will be noted that the variational procedures developed 
here are applicable to tanks with multiconnected free surfaces, 
such as an annulus as well as simple connected free surfaces. 

A more general approach to variational solution of sloshing 
has been initiated to include the nonlinear aspects of the prob- 
lem. Results of the study will be reported in a separate pape! 


Formulation of Problem as a 7 
Variational Principle 

The differential equation and boundary conditions for the 
linearized sloshing problem are developed in the literature (1). 
A formulation of the problem as a variational principle may 
be developed from these differential equations. However, 
the development of a variational principle directly from 
Hamilton’s variational formulation of mechanics is adopted 
here in order to clarify the assumptions employed in the deri- 
vation and to obtain a free boundary conditions statement 
of the problem. The subject of free boundary conditions 
in a variational formulation will be discussed after the basic 
mathematical apparatus has been presented. 

If 7 and V are defined as the kinetic and potential energy, 
respectively, of the fluid, Hamilton’s principle may be written 
as 


to te 
= f Lat (T — V)dt =0.. (1) 


The kinetic energy of the fluid under small waves can be 
expressed, to the first order 


T= Vol%dr... [2] 


where @ is the velocity potential, and 7 represents the volumé 


* Numbers in parenthesis indicate References at end of paper. 
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Fig. 1. 
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between the undisturbed surface z = 0 denoted by the sym- 
bol S and the wetted area of the tank & (see Fig. 1). 
The potential energy of the fluid can be represented by the 


following equation 
¢(pgtdS).. 


if the zero potential level is chosen to correspond to the fluid 
at rest. The symbol ¢ denotes the fluid surface height meas- 
wed from the undisturbed surface. Using the Bernoulli 
equation, it can be shown that for small waves the following 
equation holds at the free surface 


Axially symmetrical tank and symbols 


V = . [3a] 


Substituting [8b] into [3a], an expression of the potential 
energy is obtained in terms of the time derivative of the veloc- 


ity potential 
ra 
J = J, ot ds. 


Using Equations [2 and 3] for the expressions of the kinetic 
energy and the potential energy and neglecting the fluid 
lensity p, the Lagrangian may be written as 


Hamilton’s principle, Equation [1], will be applied to verify 
that Equation [4] is indeed the proper Langrangian and to 
larify the question of boundary conditions. Upon substitut- 
ig the expression for the Langrangian, Equation [4], in 
Equation [1] and carrying out the variational process there 
results the expression 
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In applying Hamilton’s principle the variation is assumed 
to vanish at the end points of the time interval. Thus the | 
last term in [5] vanishes. If the variation 6/ is to vanish for — 
all admissible variations in the potential function ¢ the follow- 
ing equations must hold 
= Oin volumerT.... 


0 on the wetted surface Z............ [7] 
n 

1 
+ = 0 on the free surface S......... [8] 
on 9g 


Equations [6, 7, 8] represent the classical formulation of — 
the linearized problem as a differential equation with bound- _ 
ary conditions (1). The conclusion to be drawn from the © 
preceding analysis is that the problem of calculating a poten-_ 
tial function ¢ which satisfies the differential equation and — 


boundary conditions of the linearized theory is equivalent 


to finding a potential function ¢ which makes the integral 


| l 9 
dt Vo|*dr — f, | [9] 


an extremum. The importance of this conclusion lies in the — 
fact that not only the differential equation but all necessary — 
boundary conditions are automatically satisfied. There 
exists an infinite number of variational principles which satisfy _ 
the differential equation. However, a variational principle | 
which satisfies the differential equation and the boundary 
conditions is unique. Such a variational principle is called — 
a free boundary condition variational principle in mathe- | 
matics literature (4). 

The well-known Rayleigh-Ritz procedure will be used to— 
obtain an approximate potential function from the statement _ 
that the action integral [9] is an extremum. 
is expanded in a series of known functions with undetermined — 
coefficients, and the action J is evaluated. The choice of — 
functions employed in the expansion has an important effect _ 
on the accuracy of the results. Since J is thereby represented 
as a function of the unknown coefficients, the usual methods 
of finding an extremum may then be used. The functions 
used in the expansion need not satisfy the boundary condi- 
tions in view of the fact that the variational principle is of 
the free boundary condition type. 

The sloshing modes and natural frequencies may be de- | 
fined as solutions of the extremum problem for the integral J 
of Equation [9] when the potential is a harmonic function of — 
time or symbolically 


y, 2, t) = y, z) sin (ot + €)....... [10] 


Here y represents a sloshing mode and o the natural fre- 
quency in radians per sec. Upon substituting expression 
[10] for ¢ in the action integral [9] and integrating over a 
single cycle in time, an expression for J is obtained in terms of 
the function y and its gradient. 


f — f y'dS......... [11] 
2g Js 

In proceeding to Equation [11], a constant factor rp/o has 

been replaced by unity. It is clear that this substitution will 

not affect the results, since the factor may be absorbed in the 

definition of the Lagrangian. 

The differential equations and boundary conditions for 
the modes and natural frequencies may be derived either by 
substitution of the expression [10] for ¢ into Equations 
(6, 7 and 8] or by carrying out the variation of the expression 


The potential 


a 


= r 

r 
|| 


for J in Equation [11]. 
are 


= Oin [6a] 
= on the wetted surface 2.......... [7a] 

Ww = 0on the free surface S........ [8a] 

g 


Upon multiplying Equations [6a, 7a and 8a] by y, integrat- 
ing over the regions of validity, and integrating by parts. 
the following result is obtained 


z+ 


The above equations may be solved for o?/g to give the result 


It can be shown by direct calculation that the frequency 
_ parameter \ is a minimum when the integral expression [11] 
for J is an extremum. 
The condition that the frequency parameter \ be an ex- 
- tremum where X is defined by Equation [12] thus represents 
an alternate statement of the variational principles with 
free boundary conditions. 
By using Equation [12], the natural frequency parameter 
d can be calculated to second order when the potential y 
is known to first order. This result follows from the fact 
_ that the first order error 6\ vanishes when the potential is 
calculated from a variational principle. 


‘ 


The Shallow Tank Approximation 


If the tank under consideration is shallow, i.e., the depth is 
sufficiently smaller than the distance between mode crests, 
the calculations can be simplified considerably. Equation 
_ [12] can under suitable circumstances be written in the form 


A=— = 


where A is the depth of the tank. It will be assumed that the 
potential y is independent of z. It is convenient to denote 
this special potential function by E(z, y) and the frequency 
parameter by Ae. Then the previous equation may be simpli- 
fied to the form 


dS( E,? + Ey?)h 


f _dSE? 
S 


The variational principle 6\, = 
boundary conditions. 


4 0 remains valid with free 


The Deep Tank Approximation 


If the tank is cylindrical in shape for depths sufficiently 
greater than the longest wave length of a sloshing mode, the 


In either case the resulting equations 


dS dz|(AE, + + (ABy + BES + \ 


= - [18 
by A? + 2bvAB + boB? 
where aj; and b;; are defined as 
ay, = JS dSWE, + 
en 


bottom contour may be ignored. 
on the surface S as expressed in Equation [8a] becomes 


ra 2 
14] 
Oz g 

when Equation [12] is employed. This boundary con- 


dition will be satisfied by any function of the form 

y = (x, y).. [15] 
where A, is the frequency parameter for a deep tank. Upon 
substituting expression [15] for the potential in Equation 
[12] and assuming the depth of the cylindrical tank to be 
infinite, the following result is obtained 


f, dS + Py? + 
2 = = 


f, dSF? 


f, dS(F, 2 + 


dSF? 
dSF? 


Equation [16] represents a free boundary variational principle 
for a deep tank which is cylindrical for a depth sufficiently 
greater than its maximum cross-sectional dimension. 


Tanks of Arbitrary Depth 

The potential functions Z and F for shallow and deep tanks 
‘an be calculated by the Rayleigh-Ritz method from the 
variational principles 6A. = 0 and 6A; = O where A¢ ani Ay 
are defined by Equations [13 and 16], respectively. The 
calculations necessary to determine these solutions are rela- 
tively short, since both cases involve only two dimensions. 
An approximate solution to the general three-dimensional 
case may be developed from the two-dimensional cases with 
little additional labor. 

The basic approach is simply to represent the potential for 
the general case W as a linear combination of the potential 
functions for the shallow and deep tanks. 

y, 2) = AE(z, y) + Be F(z, y)...... 17 
Upon substituting the preceding expression for y in the varia- 
tional principle [12], the following formula is obtained 


fi dS(A2E? + 2ABEF BYF?) 


Upon performing the indicated integration, X may be ex- 
pressed in the form 


aA? + 2a2,AB + 2B? 


ll 


gh 


2, 
= JS dSE? 
= Bn = J dSEF 


Be = J dSF? 
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It is evident from Equations [13] that the following identity 
holds 


The boundary condition 
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The frequency parameter \ as expressed in Equation [18] 
must be an extremum. Upon setting the derivatives with 
respect to A and B equal to zero, the following results are 
obtained 


(au — Aby)A + (ay — Aby)B = 
(ax — Aby)A + (a2 — = 


~ 


The condition that the above equations be consistent gives 
the following quadratic equation for X. 


ay 


Il 


Abu — Abe) — — — An) = O [21] 


Upon solving Equation [21] for A, the ratio of the coefficient 
B/A can be computed 


B (Gn — Abu) (aa — Aber) 


A (ay — Ab) (azz — Abe») 


The mode shape may then be computed by making use of 
expression [17] for y. 


Axially Symmetric Tank 


For most missile applications, the tanks are symmetrical 
about the missile axis. For these tanks, it is convenient to 
use -vlindrical coordinates. The potential function yY can 
then be written in the form 


=m t(r,z) cope @......... [23] 
a! 

= ¢(r, z) sins 6.. 


und the modes may be described by integral values of s. 
The dependence of & on s is not shown explicitly to simplify 
the notation. Upon substituting expression [23] for wy 
in Equation [11]. the action integral J becomes 


2 


r 2x £25 Pro 
[= ” der f dz aol + £2) cos?s 0+ sin? —X f drr f dE? cos? s 
—h 0 0 


l= drr + &2 + rf, drr&.. 


where a trivial constant factor + has been omitted. The 
expression [24] for J is unaltered if Equation [23a] is used in 
place of Equation [23] for y. 


The Shallow Tank 


For the shallow tank, & is assumed to be independent of 2, 
and Equation [24] is further simplified as 


dG\? at 2 
7) dr q 


vhere the potential G is a function of r only. 
To solve Equation [25], the unknown function G(r) is first 
represented by a polynomial, namely 
N 8s + 2n 


n=0 To 

This is the form of the exact solution for axially symmetric 
tanks with parabolie bottom (1). The polynomial form for 
“has been found to give rapid convergence in the treatment 
i shallow tanks. The presence of s in the exponent is also 
‘onsistent with the fact that when s ¥ 0 the singularity of 
Equation [24] at r = 0 requires that G(r) contain no constant 
term. The depth A will be represented by a polynomial 
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The polynomial coefficients a, in Equation [26] must be_ 
such that the integral J of Equation [25] has a stationary 
value. This condition requires that 


Oy, 


which produces a set of .V + 1 simultaneous homogeneous — 
algebraic equations as follows 


K 524 (n + m)s + 2nm B?/2 
[29] 
where 
= dimensionless frequency squared 


In order that the homogeneous equations be consistent. 
the determinant of the system must vanish. Thus 


where 
> (n + m)s + 2nm B?/2 
—, statmt+t/2 stntme+1 
and n and m (= 0, 1, 2.......1 VY) denote row and column | 


respectively. The solutions of Equation [30] will yield 
(V + 1) frequencies with the corresponding mode shapes 
given by Equations [26 and 29]. 


Shallow Tank Examples 


Parapo.ic Borrom In the case of a parabolic bottom where 


the exact solution (1) consists of terminating hypergeometric 
series. The variational procedure presented here reproduces 
the exact solution. 

Fiat Borrom The flat bottom case has been solved exactly 
in terms of Bessel functions. The exact value of the lowest 
radial frequency 6 for s = 0 is 3.8317. In comparison our 
approximation using a second degree polynomial is 3.8769. 
The agreement is considered satisfactory. 

ConicaL Borrom Since this case is of considerable practical 
interest, extensive numerical calculations have been _per- 
formed and are presented in this paper. The bottom has the 
simple analytical representation 


To 


Four values have been assigned to s, namely 0, 1, 2 and 3. 
A sixth degree polynomial, i.e., a polynomial with V = 3, 
has been employed in all the cases, yielding in each case 
four eigenvalues for the dimensionless frequency 8. These 
frequencies and the coefficients of the corresponding normal- 
ized polynomials are given in Table 1. The polynomials 
are plotted in Fig. 2 to give a graphical representation. 

In order to assess the effect of the number of terms of the 
polynomial on the accuracy of the results, polynomials of 
degrees from 1 to 3 have been used in the calculation of natu- 
ral frequencies. Results given in Table 2 indicate that poly- 
nomials with a very small number of terms are satisfactory, 


[22] 
: 

— 
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. Fig. 2 Frequency vs. depth for a conical bottom (s = 0) 


Fig. 4 Frequency vs. depth for a conical bottom (s = 2) 
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Fig. 3 Frequency vs. depth for a conical bottom (s = 1 
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Fig. 5 Frequency vs. depth for e conical bottom (s = 3) 


rt 
Table I Frequencies and normalized polynomial coefficients for shallow tanks with conical bottom : 
s Bsj n=0 n=1 n=2 n=3 
0 0 0 
—3.0296 2.1878 0.37387 
—11.672 0.18296 12.310 
47.902 113.74 73.286 
0.89134  —0.65843 0.27075 
—2.3593 6.1693 1.1262 
-17.930 = 13.350 27.457 
—121.99 214.19 —113.23 
2.7094 3.5354 1.8067 
3.4180 2.5563 8.7838 1.2027 
2 5.1098 26.217 29.610 44.987 
6. 8881 264.32 448.01 — 233.36 
1.8465 1.3202 —0.74504 0. 26920 
3.8511 —2.5725 —11.917 1.4863 
3 2 5.5776 2.4142 — 109.01 94.316 
3 6.9197 484.80 747.51 — 363.13 


circular mode number 
radial mode number 

coefficients of the polynomial representing the radial wave of the jth n 
dimensionless frequency 


10ode 
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6 Comparison of an intermediate wave shape with the 
asymptotic wave shapes 


| the calculated frequencies converge to the correct value 


rapidly when the number of terms increases. 


he polynomials given in Table 1 are normalized so that 


The orthogonality condition for these polynomials is 
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Table 2 Effect of degree of the polynomial on accuracy 
of the frequency of shallow tanks with conical bottom 
N j=0 j=1 j= 2 
0 0 
1 0 2.1899 
0 2 0 2.1732 3.8645 
3 0 2.1730 3.8078 5.5218 
0 1.1547 
1 1.1476 2.9370 
] 2 1.1466 2.9032 4.6177 
3 1.1459 2.9029 4.5382 6.2883 
0 1.5492 
l 1.5425 3.4483 
2 2 1.5424 3.4181 5.1843 
3 1.5424 3.4180 5.1098 6.8881 
0 1.8516 
1 1.8466 3.8771 
3 2 1.8465 3.8512 5.6742 
3 1.8465 3.8511 5.5776 6.9197 
s = circular mode 
N = degree of polynomial 
j= radial mode ded 
83; = dimensionless frequency 
The frequency sj for a given s and j shows rapid con- 
ergence with the degree of the polynomial used. 
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1 r r 
d ( ) [G5 ) = 
0 To To 


re = = 8’, and = Oif s 8’; Lifj =)’ 
6; = Oif 7] #7’, has been used as a check on the numeri- 
calculations. 


Circular Tank of Arbitrary Depth 


‘he procedure for obtaining the solution for the tank of 
trary depth is discussed in the section ‘“Tanks of Arbi- 
y Depth.” Fora circular tank, the shallow tank potential 

y) becomes G(r), a function of r only, and may be ob- 
ed by the method given in the preceding section. The 
» tank potential F has been solved (1) exactly in the form 
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Fig. 8 Normalized wave shapes for a deep tank of circular 
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of a Bessel function Js(uyr/ro), sy being the jth zero (2) 
J s'(usjr/ro), Where the prime denotes differentiation. 

For a circular tank of conical bottom with ho/ro = 1, the 
results obtained by this procedure are given in Figs. 2 through 


AMPLITUDE 


NORMAL 
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NORMAL IZED 
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Fig. 9 Normalized wave shapes for a conical tank of 90 deg 
vertex angle (ho/ro 1) 

+ 

> 


9. This happens to be a special case which has an exact solu- 
tion (3) for the lowest mode (i.e., fors = 1,7 = 0). The 
exact solution gives o%ro/g = 1, compared to 0.994 ob- 
tained by the present technique. Furthermore, the exact 
solution gives a linear function of r for the wave amplitude 
of this mode. The present technique gives a result, as shown 

in Fig. 4, in close agreement with this. —_ 

Conclusions. 

1 A variational principle for the approximate calculation 
of the normal modes of sloshing has been developed, including 
the specific formulas for the shallow tank, the deep ta ik 
and tanks of arbitrary depth. 

2 Application of this method to a shallow conical tank 
shows satisfactorily rapid convergence with the degree of thie 
polynomial used, as shown in Table 2. 

3 Application of the method to flat bottomed tanks of 
arbitrary depth shows satisfactorily accurate results. ‘I he 
error curve in Fig. 11 shows a maximum error of 0.8 per cc nt 
for the lowest natural frequency (Fig. 10) when the potential 
function of the shallow tank is approximated by a fou th 
degree polynomial. Similarly satisfactory results are obtain ed 
for 90 deg conical bottomed tanks. 

4 The variational procedure presented in this paper : p- 
pears to represent an accurate and rapid method for calcul:t- 
ing the sloshing modes and frequencies for a wide class of 
tank geometries. It must be emphasized, however, that the 
choice of functions used in applying the Rayleigh-Ritz pro- 
cedure has an important bearing on this conclusion. 
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using fourth degree polynomial approximation 
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Stagnation laminar heat transfer at hypersonic speeds 
depends on the rate of recombination of the dissociated air 
behind the detached shock wave. This paper is concerned 
with the case of a large recombination time compared to 
the time of diffusion across the boundary layer. 
ditions of existence of such a “‘frozen flow”’ and its coupling 
with the dissociation lag behind the shock are discussed. 
In order to account for finite catalytic recombination rates 


The con- 


at the wall, a nonsimilar boundary condition is introduced 
which can be reduced to similarity for stagnation flow 
In this latter case, Lees’ (1) and Fay and Ridell’s 
(2) heat transfer solutions are shown to correspond to the 
limiting case of an infinitely fast catalyst. The validity of 
their solutions is extended to the general case of a wall of 


oniy. 


finite catalytic efficiency, by introducing a correction fac- 
tor g. This factor is a simple function of the flight con- 
dition, nose geometry and the wall catalytic recombina- 
tion rate constant. For a given nose material, the per- 
centage of the heat transfer by catalysis is found to in- 
crease with the velocity, the nose diameter and the wall 
temperature and to decrease with altitude. Finally, the 
experimental values obtained for the catalytic recombina- 
tion rates of oxygen and nitrogen atoms on various sur- 
faces illustrate numerically the importance of the nature 
of the wall on the catalytic heat transfer to a missile nose. 
In particular, the superiority of pyrex over metallic sur- 
faces stresses the need for more experimental values for 
glassy and ceramic coatings. 


Nomenclature 
¢, C4 = atom mass fraction 
Cy = molecular mass fraction 2. 
Cp = specific heat of air at constant pressure, cal/gr K 
e = L~cicpi = frozen specific heat, cal/gr K 
d = nose diameter, cm oshe 
D = ordinary diffusion coefficient, em?/sec 
Dam = Damkohler number (Eq. [35]) : 
f = u/Ue 
= hs/hse 
h = enthalpy, cal/gr 
.. = frozen enthalpy, cal/gr 
he = heat of recombination (atoms — molecule), cal/gr 
(Eq. [32] ) 
j = net mass flux of atoms, gr/cm? sec 
ky = catalytic reaction rate constant (Eqs. [6 and 34]) 
l = 
Le = Lewis number = pDZép/d 
m = order of catalytic reaction (Eq. [6] ) 
M = average molecular weight of air, gr/mole a 
= atomic weight of air, gr/mole 
Vy, = molecular weight of air, gr/mole >a 
= molecules per unit volume, 
= pressure, dynes/cm? 


Presented at the ARS 12th Annual Meeting, New York, N. Y., 
Dec, 2-5, 1957 and revised, July 1958. This work has been sup- 
ported partly by the United States Army, through the Office of 
Ordnance Research, partly by the McDonnell Aircraft Corpora- 


On Catalytic Recombination Rates in 
Heat Transfer 


R. GOULARD? 


- Purdue University, West Lafayette, Ind. 


Z. compressibility factor, (» = r) 
B = du,/dzx stream velocity gradient (Eq. [23] ) 
1 = wall recombination efficiency (Eq. [34] ) 

€ = Lennard-Jones potential well a 
n = similarity coordinate 

d = heat conductivity, cal/em sec K 

= dynamic viscosity, gr/em sec 
p = density, gr/em’ 

o = Lennard-Jones collision diameter, A : 
= Prandtl number 

y = stream function 

¢ = correction factor for catalytic effects (Eq. [22] ) 
Subscripts 

A = atom 

M = molecule 

e = stream outside the boundary layer 

= stagnation 

w = wall 7 

© = ahead of shock 


the gas expands or flows past a cool body. 7 
the formulation of the boundary layer problem must account — 
for reacting gas mixtures processes (chemical rates, diffu-— 
sion) and their effect on thermal properties. 


formulation of the boundary layer equations for the case of a 
dissociating gas. 
transfer by ordinary molecular conduction and by diffusion 
of atoms, the latter releasing their energy by recombination 
in the stream or on the cool wall. 
possible: 
dynamic equilibrium) or very slow (frozen flow). 
latter case, no recombination takes place in the boundary 


tion, under the direction of Dr. H. DeGroff. 
* Assistant Professor, School of Aeronautical Engineering. 
Member ARS. 
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personic Stagnation 


Ie = heat transmission by conduction, cal/em? sec | 

qb = heat transmission by diffusion, cal/em? sec 

q = total heat transmission excluding radiation, g¢ + gp 

dr = heat transmission by radiation, cal/em? sec 

ro = distance from axis of revolution of body, em 

R = nose radius, em 

Ry = universal gas constant, 8.314 x 107 ergs/mole K 

8 = similarity abscissa 

S = Equation [30] 

Se = Schmidt number jhe 

t = time, sec = 

‘i = temperature, K 

u = velocity along the wall, cm/sec 

v = velocity perpendicular to the wall, em/sec 

w4 = rate of generation of atoms, gr/cm# sec 

x = abscissa, (distance along the wall from stagnation | 
point) 

7] = ordinate (perpendicular to the wall) 

= ¢c/Ce 


1 Introduction 


_ high stagnation temperatures generated by hyper-— 


sonic flights provoke a chemical process of dissociation — 
behind the shock, followed by a recombination process when | 
Consequently, 


Lees (1) and Fay and Riddell (2) have laid down the general 3 


They established the dual process of energy 


Two extreme regimes are 
The recombination rate may be very fast (thermo- 
In the 


layer, except possibly at the wall where catalytic reactions 
may occur. A third case allows for a finite recombination 
rate in the stream itself (mixed flow). The two extreme cases 
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present some simplifications: These permitted Lees, on the 
basis of physical arguments, to further simplify and uncouple 
the boundary layer equations and finally to present approxi- 
mate closed form solutions for these two cases. Fay and 
_ Riddell, on the other hand, reduced the equations into similar 
_ form and, in order to avoid oversimplification, programmed 
a numerical solution for all possible gas phase recombination 
rates. Although Lees’ closed results underestimate heat 
transfer, Probstein (3) has shown that they presented a good 
approximation of Fay and Riddell’s more accurate numerical 
results. This conclusion i is also verified on Figs. 2 and 3 of 
study. 

However, a study of the catalytic recombination process 
at the wall leads to the conclusion that both Lees’ and Fay’s 
- assumptions of wall recombination with zero atom concen- 
tration represents only the case of an extremely active cata- 
lyst. In general, chemical rules require that the rate of re- 
-combination at the wall be a power function of the atom 
concentration ((4) p. 140). As shown by Chambré and 
_ Acrivos (5) in the case of incompressible flow past a flat plate, 
_ the exact wall conditions generally destroy similarity solu- 
tions. A more exact representation of the catalytic process 
is of special importance, since heat transfer reduction by 
means of surfaces of low catalytic efficiency may be an at- 
_ tractive possibility for a safe satellite re-entry (6) or for long 
range hypervelocity vehicles. 

This paper is therefore concerned with the effects on stag- 
nation heat transfer of finite catalytic recombination rates 
frozen flow. 

Although catalytic recombination can be included in the 
“equilibrium” flow solutions, the high wall temperatures re- 
- quired for its effect to be appreciable are not practical, and 
this case is not considered here. 


2 The Laminar Boundary Layer Equations in 
Frozen Hypersonic Flow 


A frozen process is one in which the chemical reaction time 
is very long compared to some characteristic mechanical time 
of the stream. In such flows, the atoms generated by the 
detached shock wave diffuse towards the cool nose surface 
recombining (w, = 0). Therefore, the recombination 
energy release depends entirely on the catalytic properties 
of the wall. Hence the importance of such frozen regimes 
when catalytic processes are under study. 


The Conservation Equations 


Prandtl’s dimensional analysis reduces the equations of 
conservation of atom mass, mixture mass, momentum and 
energy to the classical form 


de de dc 
+v-—j]=- (op =) + mass 
ox oy oy 


| 
> (pur) > (por. = 0 mixture mass 
Ou + du _ dp, + °( ou 
oh re} or re) Setar 
Oy = oy ( =) oy (00 oy 
dp du \? 
Grtu > energy 


The basic assumption of frozen flow yields immediately: 
= 0. It was also shown by Lees that enthalpy variations 
ean be broken into a chemical energy variation due to the 

_ change in mixture composition and an energy variation due 
to the variation of temperature of the existing components. 


For a binary atom-molecule mixture 


dh = tp dT +hgde 


‘< 


It is further convenient in the case of the frozen flow where 
no reaction takes place in the boundary layer proper, to intro- 
duce, after Lees, a frozen stagnation enthalpy h, excluding the 
chemical energy that would be liberated by recombination 
of the atom fraction (hzdc) 


If we then multiply the momentum equation by wu and ald 


it to the energy expression, we obtain after simplification in 
expression similar to Von Karman’s total energy expression 


- 
~ dy dy oy oy\ 2 


The Similarity Transformation 


oh, dh, 
or oy 


For compressible viscous flow, a transformation origina ly 
devised by Doronitzin has made possible a set of similarity 


solutions for the compressible boundary layer equations. 

In Lees’ adaptation of Doronitzin’s transformation to the 
case of a blunt body of revolution, two new coordinates : re 
introduced 


u, 
= To pay 
(2s)'4 Jo 


4 
0 


In addition, the continuity equation is identically satisfied if 


we define a stream function yw such that 


= puro v, = — puro 


If the dimensionless stream function is defined, f = w/(2s 
it is then easily verified that 


A dimensionless frozen stagnation enthalpy is introduced: 
= h,/hse, and also a dimensionless atom concentration: 
= ¢/c,, where both stream values /,. and c, can be assumed 
constant if thermodynamic equilibrium is achieved in the 
stagnation area (see the second part of section 5 of this 
article). 

After substitution of the new variables in the equations, we 


obtain 
2s Ou, De 
u, Os \p 


xo 


(Uf")’ + Sf’ 


This set of equations is similar to Lees, and, provided we also 
assume in this paper a very cool wall (7, < T.;), the same 
_ boundary conditions exist for both momentum and energy 
equations 


f(0) =f’ (0) = 0 


However, the last equation (atom mass conservation) in- 
cludes a nonsimilar term in s. This dependence on s is «ue 
in general to the nonsimilar form of the boundary condition 
due to the catalytic process taking place at the wall. 
Strictly, the momentum equation also depends on z through 


f(»)= 1 g(o)=1 
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tions. 
the contribution of p to the momentum equation is neglected, 
and the terms in s have been left out for simplicity. 


Catalysis at the Wall 


l)xisting solutions (1, 2) assume a zero concentration to ey 
a suitable condition at the wall. This does not account for 
the chemical process of catalytic absorption (4,7) where the 
catalytic reaction rate is proportional to a power m of the 
concentration ¢c,p, at the wall.4 Since the net mass flux of 
atoms diffusing towards the wall, neglecting thermal diffu- 


si0}:, is 


it follows from the conservation of mass principle that in the 
ste: dy state 


Ju = Po De ( ‘) = k{c. po)” <2: 
OY] » 


‘+h becomes the boundary condition of Equation [3] at the 
wail. It is clear from this equation that no diffused atom 
ean reach the wall while the concentration c, is zero, unless 
the catalytic reaction rate constant &, becomes infinite. 
Lecs and Fay have thus taken for their treatment of diffusion 
as the rate controlling process, the limiting case of k, ~ @. 
The other extreme case of a noncatalytic surface (kh, ~ 0) 
implies, since c, is bounded, that the net flux 7. of atoms 
toward the wall is zero; it also follows from Equations [8 
and 16] that the trivial solution c(n) = c, applies to this case. 

In general, if we transpose the wall boundary condition 
[5] into the new coordinate system 7 and s, it becomes 


2s)2 
[8] 
on ToU, Pw D, 


(28)* ky 
—_— m [7 
role Dy 


This condition destroys in general the similarity of the 
equation. A nonsimilar solution to this problem was devel- 
oped by Chambré and Acrivos (5) in the case of an incom- 
pressible flow past a flat plate. In their solution, Chapman 
and Rubesin’s and Lighthill’s methods were transposed to 
this problem. In the case of stagnation flow, however, a 
simpler similar solution is possible. 

The Stagnation Flow With Ordinary Diffusion 

The solution of the system of three equations, simultaneous 
partial differential, nonlinear and coupled to each other is 
in general an impossible problem, short of numerical solu- 
In order to make possible an analytical solution, 
Lees (1) has introduced a major simplification to this prob- 
lm: For a highly cooled wall (g.<< 1), the third term of 
both momentum and energy equations can be dropped with 
good approximation. This procedure reduces the coupling 
of these two equations to the dependence of the thermal prop- 
erties | and on temperature. Lees further ad- 


tions.5 


‘This expression of surface recombination excludes the oppo- 
site surface dissociation reaction and implies therefore a high re- 
tombination activation energy and generally a low temperature 
wall. Scala (8) has shown that this wall condition was correct up 
‘o T,, = 2000 K. A more general approach is suggested by 
Rocner (9) in the case of small deviation from equilibrium at arbi- 
trary high wall temperatures. 

*Such numerical computation has been independently carried 
ut by Seala, including thermal diffusion effects. His results have 
ee) recently declassified and are briefly discussed in (8 and 10). 
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p, and the nonsimilar term should be shown in all three equa- — 
However, the solutions we will discuss are those where 
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Fig. 2. Stagnation heat transfer at altitude 200,000 ft (assuming — 


the outer edge of the boundary layer in chemical equilibrium) 


vanced that neither / nor & were too sensitive to temperature 
for most of the boundary layer and could be assumed constant 


= 0.715(Eucken) l=1 


Lees’ last assumption (1 = const) has been long successful in 
low temperature boundary layer theory where its adoption 
leads to a convenient uncoupling of momentum and energy 
equations. It was shown however in (11 and 3) that its 
validity decreases for higher temperatures, as can also be seen 
on Fig. 1 of (2). Fortunately, enthalpy wall gradients and 
concentration wall gradients, and consequently heat transfer, 
are weak functions of |. This point is illustrated on Figs. 2 
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Fig. 3 Stagnation heat transfer at altitude 250,000 ft (assuming 
the outer edge of the boundary layer in chemical equilibrium) 


and 3 of this paper, where in the case of a catalytic wall 
(kw ~ o) the exact numerical results of (2) lie within less 
than 20 per cent of Lees’ results, even for very high stream 
temperature, i.e., large variations of | across the boundary 
layer. Lees’ assumption / = 1 retains therefore a good ap- 
proximation of the stagnation heat transfer, where the in- 
fluence of the other physical factors are available in closed 
form, 

Bs The determination of the shear force and total frozen 
_ enthalpy reduces therefore in first approximation to the solu- 
a tion of the following set of equations : 


+ ff’ 
9” + afg’ = 0 


Both equations have been solved by Blasius and Polhausen. 
In particular, with the present definitions of f, g, and s, one 
. 

4 obtains 


0 


= (0.664) = [10] 
Wall Catalysis at Stagnation 


For stagnation flow conditions, we introduce the local 


approximations 
ro(z) = z u, 
where the parameter 6 depends on flight conditions and nose 


diameter only (see Eq. [23]). 
Substituting into the nonsimilar term of the boundary condi- 


4 tion [8], we read 

Ba? 


rou, 


Hence, for stagnation flow, the boundary condition [8] can 
be written 


on w 2B 
The only s-dependent variables left are the stream prop. r- 
ties p., we ande,. But the product u.p. can be assumed 
stant and equal to for the high temperatures and nea ‘ly 
constant pressure of the stagnation stream flow (see for n- 
stance Fig. 1 of (2)). On the other hand, c, has been assun °d 
constant in the stagnation area. Furthermore, the pr:c- 
tical necessity of a cool wall implies a catalytic reaction of 
the first order: m = 1.6 In this case, the term in ¢, (is- 


ky 


appears from Equation [13] altogether, and the bound: ry 
conditions of Equation [3] are 


The Continuity Equation Solution 


The continuity Equation [3] being entirely similar at st: g 
nation, the nonsimilar term in s can be dropped 


Scfz'’ +2” =0 


Integrating this equation in two steps gives 


fdn 


z(n) — 20) =2 (0) dn 


Se 
f"(0) 


The determination of 2(0) and z’(0) comes from the two 
boundary conditions [14 and 15]. When 7 ~ o, the numeri- 
cal value of the integral in Equation [16] was found by 
Polhausen’ to be the inverse of 0.47 Sc’4. Hence 


z’(0) = 0.47 Sc’4 [1 — 2(0)]............ [17] 


Substituting this value of z’(0) in Equation [14] gives the 
dimensionless concentration at the wall 


1 


(Bae) ky 1 
28 0.47 Sc’ p,, D,, 


The ‘eenoteiieilian profiles z(m) are then given by Equation 
[16] using the constant values [17 and 18] for 2(0) and z’(0). 
It is easily verified in expressions [17 and 18] that an infi- 
nitely fast catalyst (k. ~ ©) corresponds to 2(0) = 0 and 


2'(0) = 0.47 Sc’S (Lees and Fay’s values), whereas an in- 
finitely slow catalyst (k. = 0) gives z(n) = z (0) = 1, and 
2’(0) = 0. 


3 Heat Transfer to the Wall 


The heat transferred to the cool wall, excluding radiation, 
results in part from ordinary molecular conduction and in 
part from recombination at the wall of the atoms diffusing 
through the boundary layer. The rates of these two proc- 
esses, per unit time and area, will be called q, and qp, respec- 
tively. 

* In addition to this general statement provided by Laidler ( (7 
pp. 180-181) the first order reaction assumption is verified ex- 
perimentally for oxygen by Linnett and Marsden. (See section 4 


of this article.) 
7 The constant was the Prandtl number zg, and the vari:ble 
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Conducted Heat 


This fraction of the total heat is — 


oT (dh, 
= = h, e 
oy w Cp u oy w Cp w 


since at the wall h = h,. 
Substituting the expressions of 7 and g’(0) (Eqs. [10 and 12]) 


on 

g’(0).. . [19 

(5) 


A, 
Ge = — (2B)? (0.47 h 
de = 0.47 (2B 
[20] 
sine l, = = “= 1 in Lees’ approximation. 
MsPs MesPes 


Heat Released by Recombination at the Wall 


1 expression of the heat released by recombination is 
easy obtained by multiplication of the net atom mass flux _ 


to ‘he wall 7. (Eq. [4]) by the heat of recombination hp 
Ex. [82]) 


Substituting the expression [18] for 2(0) where ' 
iis 
Mw Sc 
dp = hrC Py = 
28 0.47 


0.1 
q = 0.763 (h,, — hy) ( 


sem se 


tearranging 

Qn = 0.47 (2Bu,, p,,)'2 Se 

(2B Bea) Se 0.47 Se 

Puky 

[21] 


The interest of this formulation lies in the fact that the effect 
of the catalytic rate k, is concentrated in the last fraction. 
When k, ~ , this fraction tends to unity, and qp(eo~ @) is 
equivalent to the previous solutions. This fraction is there- 
fore a correction factor ¢* 


1 
Pky 
The expression of gp becomes 
Gp = 0.47(2B [24] 


* Where 6 may be obtained, for instance, from Probstein (12) 


d Pes 


or ‘rom Li and Geiger’s (13) expression which uses the shock di- 
ameter rather than d. Probstein’s expression is verified by the 
experimental results of (14). 
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Total Heat Transferred to the Wall 


The total heat transferred to the wall is the sum of gq. and 
dp (Eqs. [20 and 24]) 


= + = 0.664 (Bue, Pes)? hye + peg). [25] 


Introducing now the Lewis number Le = @/Se and the 
stream stagnation enthalpy hy. = hee + Agee, We rearrange 
Equation [25] into 


q = 0.664 (Buse Pye)? hye + (Le? — 1) . [26] 


se 


where the correction factor ¢ is defined in expression [22]. 
This relation is easily converted into a Nusselt number rela- 
tion where Nuw = — hw) and Rew = Mw 
based on the wall conditions. One obtains after simpli- 
fication, still assuming = 1, he >he 


N 
= 0.6642," + (Les —1) 


Discussion of Results 


Expression [26] corresponds to Lees’ results when com- 
bining his Equations [19 and 12a] (1). The Nusselt oa. 
is also in agreement. Only the catalytic factor ¢ is added. 

Figs. 2 and 3 illustrate, for two typical flight conditions, the 
effect of flight velocity on the heat transfer q, for different — 
values of the catalytic recombination constant ky. It is | 
important to note that thermal equilibrium has been assumed — 
at the edge of the boundary layer behind the shock. The — 
validity of this assumption is discussed in section 5 of this 
article. 

The more exact numerical results of Fay and Riddell, as — 
correlated for the case of frozen flow by relation [8] of (9) 


Use 


is also illustrated on Figs. 2 and 3 for the two flight condi- | 

tions considered. We verify Probstein’s comment (3) that — 

for kw ~ ©, the assumption pwpw/MesPee = 1 reduces the © 

values of heat transfer by 20 per cent at most. It further — 

suggests that Fay and Riddell’s equation could be easily _ 
fitted with the correction factor ¢ to account for wall catalytic | 

effects. 

To investigate more closely the importance of the catalytic | 
correction factor ¢ on stagnation heat transfer, it is con- 
venient to analyze the ratio of the heat transferred in the case 
of any catalyst in general over the heat transferred in the case 
of an infinitely fast catalyst qz,~. (previous solutions) 
hpe, 


1 + —1) 
h 


= 
Use 


h,. 
[29] 
1 + (Les —1) 


se 


Up ww? 


However, the catalytic correction factor ¢g is not a func- 
tion of the catalytic reaction rate constant k, only, but also 
of the flight conditions (see Eq. [22]). It is therefore 
interesting, from Equations [22 and 29], to represent the 
variations of @ directly in terms of k,, for a given set of flight 
conditions.’ The results of such calculation are illustrated 


_ * The transport properties used in this calculation are discussed 
in section 4 of this article. 
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Fig. 4 Reduced heat transfer g = ¢/qkw ~ at 200,000 ft 


in Figs. 4 and 5 for two given altitudes, one nose temperature 
and diameter. Some typical values of k, are shown on this 
figure; they correspond to the experimental catalytic data 
discussed in section 4. 

These numerical values illustrate the importance of the 
wall catalytic effect in typical cases, since the proper choice 
of the wall material seems to be sufficient to reduce the stag- 
nation heat transfer by a large factor, assuming equilibrium 
o composition of the shock layer at the boundary layer edge. 


_ The Factors of Heat Transfer Reduction in Catalytic 
Frozen Flow 


In addition to the obviously determinant role of kw, these 
~ results and figures point out the effect of the following factors: 
An inerease in flight velocity u. increases both heat 
transfer q and the fraction ¢ of this heat due to dissociation. 
- A minor deviation from this trend is observed for qin the range 
of velocities where oxygen is nearly all dissociated but nitro- 
gen is still molecular (around 13,000 ft/sec). For these local 
conditions, a lesser fraction of the shock energy is converted 
into chemical energy, resulting in a lesser fraction of energy 
release by recombination. 


2 Ifwe rewrite Equation [22] in the form 


= 
where 


| — (? [2 |) [30] 
Pw d \Pes Pes 


it can be seen that the real catalytic variable is the ratio 
S/k., and that any increase in S will have the same effect as 
a decreased catalytic reaction rate constant k,. Therefore, 
a reduced nose diameter d will create a reduction of the cata- 
lytic effect and a diminution of the heat transfer. 

3 The two terms of S/k.~ which are affected by a wall 
temperature T,, variation are p, and k,. Across the boundary 
layer, the pressure is constant, so that a. 


hence p~ is proportional to 1/Z.w7’. for a given stagnation 
state. The variation of k, with temperature is discussed in 
section 4: ky is proportional to 7’,,°% for pyrex and a much 
stronger function of 7, for some oxides. Finally, it is easily 
verified by a short iteration process that Z.. (and z., from 
which it is calculated) is a very weak function of kw; Zew 
doubles at most when k,, increases a hundredfold; hence the 
variations of Z... with temperature can be neglected in com- 
parison with those of kw. 

The conclusion is that the ratio S/kw, being proportional to 
1/pwkw, increases when the temperature decreases. The in- 
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Fig. 5 Reduced heat transfer = 4/qkw ~ at 250,000 ft 


crease is mild for pyrex (proportional to 7',’2) and strong | or 
oxides. Correspondingly, Equation [30] shows that tie 
reduced heat transfer g decreases. Hence, wall cooling (: s- 
pecially for metal walls) has the double effect of evacuati ig 
the heat convected to the wall and preventing chemi :l 
energy from releasing more heat at the wall. Conversely. if 
the wall temperature is allowed to increase, heat transfer by ».- 
combination will increase; this inherently unstable situation 
is discussed by Rosner in (15). 

4 Finally, it is difficult to analyze in all cases the effect 
of a change of altitude on the catalytic heat transfer proce-s. 
This is because nearly all the properties in g are affected by a 
change in altitude. However, it can be stated in gross ap- 
proximation that for partially dissociated air, above 16,000 
ft/sec 


Pes Pes Pes 
and from Feldman’s charts (19) 
Po Po’ 


Hence, the simplified relation 


14 
~ 


R 
and after substitution into q and ¢ 
‘4 
and 


Relation [31a] has been already established by Kemp and 
Riddell (6). Relation [31b] illustrates the dependence of the 
catalytic parameter g on the ambient density p,, a strong 
function of altitude. The beneficial effect of altitude on heat 
transfer reduction by use of a low catalyst is then substan- 
tiated. It can be seen for instance that an increase of density 
by a factor of 4 due to a change of altitude from 250,000 ft to 
200,000 ft will have the same effect on ¢ as an increase of ky 
by a factor of 2 at constant altitude. This result is illustrated 
by the massive shift along the k,. abcissa of the g curves 
when going from Fig. 5 to Fig. 4. 

If the altitude is further reduced from 250,000 to 100,000 
ft, the density would be increased more than a hundredfold, 
and the q curves of Fig. 5 would be shifted by one order of 
magnitude to the left. In this case, all substances with k,, > | 0 
would behave as nearly perfect catalysts, as was concluded 
by Seala in (8). 


JET PROPULSION 


to 


Ni 


\ 
y j 
* = 
at 
: 
. 
Sey 
~ 
> . 
: 
rel 
2 


tie 
ating 
mi al 
ly. if 
ation 


effect 
CESS. 
by a 
6,000 


[Bla] 
_[31b] 


p and 
of the 
strong 
n heat 
bstan- 
lensity 
10 ft to 
e of kw 
strated 
curves 


00,000 
edfold, 
rder of 
10° 
cluded 


JLSION 


vatalytic character of the wall improves with altitude. 
Seala’s conclusion, illustrated only for an altitude of 100,000 
ft in (8 and 10), is therefore too restrictive for the possibilities 
of heat transfer reduction at higher altitude. This remark 
is especially important in view of the fact that frozen flow 
is nore likely to take place at higher altitudes. 


4 Properties of Dissociated Air 
Equilibrium Properties of Dissociated Air 


Che state and composition of the gas in equilibrium behind 
th shock is available from many sources. A very conven- 
ie: t reference is the set of charts published by Feldman (19). 

lilsenrath and Beckett’s (20) composition of argon free 
air has been adopted: 0.78847 moles of nitrogen and 0.21153 
moles of oxygen per mole of air. The standard heat of re- 
co nbination of nitrogen is 8089 cal/gr (9.76 ev) and that of 
ox gen is 3686 cal/gr. It has been assumed, after Hirsch- 
fel ler (21) and Lighthill (22), that the number of degrees of 
freedom of molecular air was double that of atomic air, re- 
su ting in the same specific heat per unit mass: c,y = ¢,4- 
Tle heat of recombination is therefore simplified to its stand- 
ar value 


hp = + —Cpa) (T — T°) = hop 

Finally, it is convenient to observe that practically all of the 
oxvgen is dissociated before nitrogen begins to dissociate. 
The atom concentration c, can then be used directly in the 
dissociation energy term hac, 


for c, <0.2346 hee, = hro, ¢, (oxygen only ) 


for c, >0.2346 hpe, hro, (0.2346) + hen, fe: —0.2346) 
(oxygen + nitrogen) 
[32] 
Transport Properties of Dissociated Air 
po OF; ss et 


Classic experimental techniques fail to give high temper- 
ature experimental values for transport properties of air. 
Although new approaches to this problem are being developed 
(i.e., shock tube techniques) the bulk of the recently published 
values rely entirely on extrapoletions of low temperature 
data based on molecular theories of gases (16,17). Since 
little is known about the atom-molecule interaction poten- 
tials of oxygen and nitrogen, these extrapolated values must 
be considered only as best available guesses. 

Among recent estimates, Lees’ values (1) were calculated 
by molecular theory methods on the basis of the following 
assumptions: 

| High temperature 
of atoms and molecules. 

2 The molecular models of air atoms and molecules are 
of the Lennard-Jones 6-12 type, with the following collision 
diameters and potential wells 


air is assimilated to a binary mixture 


= 0.80 A 


€4/k = 


Lees estimated these numerical values by extrapolation of 
atom-molecule hydrogen data. His molecular model has 
been adopted here. 

Because of the present uncertainty on the Prandtl, Schmidt 
and Lewis numbers, the following constant values, used by 
Lees, Fay and others, have also been used here 


= 3.34 oy = 3.69 A 


1.14 108 = €au/k = 84 °K 


0.485 Le = 1.4 ¢ = 0.71 


Sc = 


Vinally, the value of viscosity at high temperatures had 
to be calculated. In this paper, Wilke’s semi-empirical 
relation ((18) p. 533) was used. The lack of collision in- 
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It is however apparent from this discussion that the non- | 
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Fig. 6 Catalytic efficiencies for atomic oxygen and atomic 


nitrogen. The lines of constant k, (Eq. [34]) are practically 
the same for oxygen and nitrogen 


tegrals for the atom species, due to its very large potential 
well, made it necessary to extrapolate the available collision 
integrals to lower reduced temperatures. This approach, 
although approximate, is more consistent with Lees’ esti- 
mates than using a rigid sphere model for the atoms while 
retaining a Lennard-Jones model for the molecules. 


The Catalytic Recombination of Dissociated Air 


THE RECOMBINATION EFFicIENCy y The catalytic reaction 
rate constant k, defined by Equation [5] is not the constant | 
usually chosen to measure a catalytic process. Most chem- 
ists rather use the recombination efficiency y defined as the 
ratio of the number N, of atoms recombining on a surface per — 
unit area and time to the total number N of atoms striking — 
the surface per unit area and time 


N 


The number N, of atoms recombining upon hitting a unit — 
surface per unit time is easily derived from the mass flux at 
the wall. For a first order process, Equation [5] vields 


Ms ms 

The total number N of atoms hitting a unit surface per unit 

ys is easily calculated assuming a Maxwellian distribution — 
((23) Appendix B). In the form given by Langmuir (24) 


where k is Boltzmann’s constant. And since pa = 
Cap = nama the catalytic efficiency y can be written 


N, ky, (20 my kT Ty? 
N ngkT ki 


[34] 


nakT and 


RECOMBINATION OF ATomMIC OxyGEN Linnet and Marsden 
(25) have studied the mechanism of recombination of atomic — 
oxygen on various surfaces. They found that recombination — 
on both glass surfaces and oxides was a first order process. 

They further established some experimental values of the — 
‘atalytic efficiency y in a temperature range from 0 to 400 C. — 
These values are illustrated on Fig. 6, as well as the additional — 
values given by Linnett in Table 1 of (10) and the qualita-_ 
tive results of Reeves and Harteck (26). For convenient 
reading of heat transfer charts, lines of constant k,, have been 
‘alculated from Equation [34] and are also illustrated. 

The catalytic efficiency of several tested materials is low 


= 
ft 
: 
‘ 


enough, in the range of temperatures considered, to keep 
kw < 10. As can be seen from Figs. 4 and 5 such surfaces 
should allow a substantial reduction in heat transfer. This is 
especially the case for pyrex, for which Linnett’s measurements 
can be represented” between 300 and 700 °K by 


y™~4xX107T (T in °K) 


E: Other substances, such as the chlorides and oxides investi- 
‘(i gated by Linnett have a much stronger positive temperature 
dependence. Consequently a higher wall temperature, while 
reducing the temperature gradient and the conducted heat, 
_ will increase the recombination heat release at the wall. The 
possible consequences of these opposite effects on heat trans- 
fer instability have been recently discussed by Rosner in this 
Journal (15). 

Also note the possibility of “poisoning” an otherwise 
catalytic surface with some of the low efficiency elements 
shown on Fig. 6. As pointed out by Wahlen (27), the ad- 
sorption process suggested by Linnett and Marsden would 
be very sensitive to such poisoning. 

RECOMBINATION OF ATomic NirROGEN The recombination 
of atomic nitrogen, on various metals in a wide range of 
- temperatures, has been the object of a series of articles by 
~ Roginsky, Buben and Schlechter (28). Most metals tested 
— (Pt, Cu, Fe, Ag) seem to be highly catalytic for atomic nitro- 
gen in the whole range of temperature, while a temperature 
_ dependence is noted for nickel. Recently, the recombination 
efficiency of nitrogen on a clean pyrex surface has been 
: measured by Wentink et al. (29), and a value of y = 3 X 
- 10-* has been obtained. It is interesting to note that the 
same metals are highly catalytic for both atomic oxygen and 
: £ nitrogen while pyrex is nonecatalytic for both. 
7” It should also be noted that all these measurements were 
made for partial pressures of atomic nitrogen appreciably 
smaller (<10~ atm) than for most flight conditions of in- 
terest (even at 250,000 ft, u.. = 20,000 ft/sec, partial pressure 
of nitrogen is 2.8 X 107? atm). Catalytic experiments at 
higher partial pressures are needed. 
RECOMBINATION OF DissociaTED Air Although the catalytic 
- recombination values discussed above for pure oxygen and 
nitrogen show the overall trend of the catalytic phenomena, 
additional complexity can be expected in the study of the 
recombination of a mixture of oxygen and nitrogen. As in the 
case of the gas phase recombination discussed in section 5, fast 
intermediate reactions between oxygen and nitrogen might 
very well accelerate the recombination and increase the ‘“‘bulk”’ 
catalytic efficiency y. No experimental result has been Te- 
ased in this area of research. 


Ee 5 The Conditions of Existence of Frozen Flow 
Recombination in the Boundary Layer 
To determine whether atoms will recombine when diffusing 
_ across the boundary layer, it is useful to compare, after Fay 
and Riddell (2), the average time fairs for an atom to diffuse 
to the wall and the average lifetime tenem of this diffusing 
atom. A dimensionless Damkéhler number can then be 


For large values of this number (Dam > 1), chemical 


taitt 


Dam = 


chem 


; It is to be noted that the temperature 7’ used by Linnett is 
Br. actually the temperature of the particles hitting the surface and 
not the surface temperature itself, which can be different at the 

a Bon low pressure of his experiments. At the relatively higher 
i _ recovery pressures in hypersonic flight (i.e., partial pressure of 
oxygen = 7 X 10~* atm at 250,000 ft., 20 000 ft/sec) no tempera- 

- ture jump is expected at the wall and Te = T particles. We there- 

_ fore assume implicitly in Fig. 6 that the energy of the particles is 


the governing factor in their catalytic recombination and not the 
- energy of the wall fixed particles. 


reactions instantly adjust the diffusing particles to local ther- 
modynamic equilibrium. For small values (Dam < 1), no 
recombination takes place before the diffusing atoms hit the 
wall. By expressing the Damkéhler number in terms of flight 
conditions, it therefore becomes possible to determine the 
altitude and velocity range in which frozen flow and cataly tic 
phenomena will prevail for a given vehicle. 

To obtain this result, a simple formulation of taits and fexem 
is required. The diffusion time tais: is easily calculated from 
classical boundary layer theory. Unfortunately, no relia)le 
expression of the chemical time is available, essentially |e- 
cause the chemical process of air recombination itself is p:es- 
ently not well understood. 

As a first step, Fay and Riddell assumed a three-body re- 
combination model 


atom + atom + particle — molecule + particle 


with a recombination rate constant 


—1.5 
k, = 5x 104 (— 


This basic assumption led to the definition of a “recombinat ion 
rate parameter” C similar to the recombination Damk6l ler 
number.!! Frozen flow can then be shown to exist at hig'ier 
altitudes and lower velocities, as, for instance, in the case of 
a re-entry vehicle of nose diameter d = 100 cm, at 250,100 
ft altitude. 

Unfortunately, later experimental results indicated tiiat 
faster air recombination rates occur at high temperatures 
(30, 31), and it was soon suggested (32) that fast intermediate 
reactions between nitrogen and oxygen were responsible for 
the observed high recombination rates. Finally, theoretical 
work (33,34) indicates the possibility, at low temperature, of 
recombination rates lower than the present experimental 
high temperature values. Considering this rather fluid state 
of information, there is little to be done at the present time to 
predict the expected range of existence of frozen flow with 
much accuracy. 


cc? mole~? sec™!.. . 5) 


Dissociation in the Shock Layer 


Because of the interest of high altitudes for frozen flow, 
some of the problems of low density flight should be briefly 
examined. 

As Adams and Probstein have recently shown (35), the 
boundary layer thickness at stagnation increases strongly 
with altitude (~p,,~”), whereas the shock distance from 
the nose is a rather mild function of altitude (~p,~”). 
As a result, the inviscid shock layer will tend to disappear 
above an estimated ratio of the mean free path over the nose 
radius of 1/100, and in such case, the basic relations of this 
study should be reconsidered. At 250,000 ft, the mean free 
path is of the order of a few millimeters, and this problem may 
arise at this altitude in the case of small diameter noses. 

On the other hand, although the shock thickness is still 
very small in the classical sense of translational collisions, 
large dissociation relaxation times may prevent the particles 
from being at thermal equilibrium when reaching the outer 
edge of the boundary layer. For this problem, we again 
have to deal with the yet unsolved subject of air ehemical 
kinetics. In general, it can be expected, as shown by Camac 
et al. (36) in the case of pure oxygen, that low densities and 
velocities will cause dissociation lag (in Camac’s experi- 
ments, the relaxation time is inversely proportional to the 
density ahead of the shock p,). 

It is interesting to note, after Adams, that the same condi- 
tions which which are favorable for frozen recombination (low 


uc 1 (Eq. [ [57] of (2)) is different from the Damkohler number, 
since it is s baaed on flight and nose geometry conditions only and 
is not inclusive of wall conditions. For instance, a large increase 
in wall temperature and catalytic efficiency would accentuate the 
frozen character of the flow (lower Damkéhler number) without 
affecting C. 
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density and velocity) are also those that promote frozen 
dissociation. Depending on the chemical constants in- 
volved, the atom concentration of the shock heated air, 
reaching a frozen boundary layer at stagnation, may well be 
only a fraction of its equilibrium value at that temperature. 
Consequently, only a fraction of the heat transfer reduction 
illustrated on Figs. 2 to 5 may actually be achieved by using 
noncatalytic walls, the rest of this energy being conducted 
to the wall by molecules which did not have time to disso- 
ciate. However, such coupling is not necessarily rigid, since 
shock dissociation is a high temperature process, whereas 
boundary layer recombination is a low temperature one, 
possibly much slower (33, 34). 

\foreover, the stagnation area extends practically several 
shock distances from the geometric stagnation point. This 
situation multiplies the time available for complete stream 
dissociation by an appreciable factor. In some conditions, 
it will be possible, having a larger stream atom concentra- 
tion away from the stagnation point, to accordingly reduce 
he it transfer by use of a noncatalytic surface. Such a remark 
is naturally only qualitative, since to drop the assumption 
c, = const. introduces a nonsimilar boundary condition in the 
atom mass conservation equation; the similar solution pre- 
serited in this paper does not apply to this case. 
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Technical Notes 


Liquid Propellant Inertia and Damping 
Due to Airframe Roll 


HERBERT REISMANN' 


ar, _ The Martin Company, Denver, Colo. 


Nomenclature 


v = fluid speed 
r = radial coordinate 

a = radius of fluid cylinder Pip 
t = time 

w» = frequency 

v = fluid speed at surface of fluid cylinder (at r = a) 
6 = mass density of fluid 
= viscosity 


vy = kinematic viscosity 
p 


A = a, dimensionless radius of fluid cylinder 


= fluid shearing stress 
I = effective moment of inertia ae 
k = damping coefficient i 

Ant 
Introduction 


ALLISTIC missiles which obtain their propulsive thrust 

from a liquid propellant engine are characterized by a 
high mass ratio; i.e., the weight of the liquid propellant repre- 
sents a high percentage of the total take-off weight. If one 
considers the rolling characteristics of such a missile in con- 
junction with, say, the design of a roll stabilization and/or 
control system, the following question must be answered: 
What is the effective moment of inertia, and how much 
damping is provided by the fluid? It is intuitively obvious 
that not all of the fluid will participate in the rolling motion 
of the airframe, and that the fluid stresses set up by this mo- 
tion will result in a damping force acting on the airframe. 
Thus it is the purpose of this note to answer the following 
questions: To what degree does the liquid propellant in a 
rolling missile participate in its motion; what is its effective 
moment of inertia, and how much damping does it provide? 


Analysis 


We chose as our mathematical model a smooth walled tank 
of circular cross section and unit depth. Because of the 
manner of excitation (roll), it is reasonable to expect the fluid 
to respond in circular motion about the center of the tank. 
Thus the motion of a fluid particle is a function of its radial 
coordinate and time only. 

The motion of a viscous fluid is characterized by the 
Navier-Stokes equations (see, for instance, (1))?._ Neglecting 
compressibility, and assuming the motion to be independent 
of the polar angle, the Navier-Stokes equations in polar coor- 
dinates reduce to the following single partial differential 
equation (in dimensionless form) 

lov 1 ov v 


w Op? pop p 


Received March 26, 1958. 
! Systems Analysis Staff, Denver Division. Member ARS. 
? Numbers in parentheses indicate References at end of paper. 


where the fluid shear stress is given by 
&. 
y Op \p 


For harmonic excitation of the missile in roll, we assume a 
velocity distribution of the form 


v(p, t) = F(p) cos wt + G(p) sin wt........ 


Substituting Equation [3] into Equation [1] and equating 
coefficients of cos wt and sin wt, respectively, we obtain 


dp? p dp p* i 
dp? pdp p? 
Equations [4] can be combined into a single linear diffc r- 
ential equation by letting 
W=F 
We thus obtain _ 
2Ww 1dWw 1 
dp? p dp 


Squation [5] is the well-known Bessel differential equation 
of the first order. Its solution is given by 


Wp) = Al(pi~'/?) + 6 


where A, B are arbitrary constants of integration, and J; and 
K,; are known as the modified Bessel functions of the first and 
second kind, respectively, of the first order. Since it is awk- 
ward to manipulate Bessel functions with complex argu- 
ments, we note the following definition (see (2)) 


= ber, p — i bei p 
iK,(pi~'/?) = ker, p — i kei p 


The defining series expansions for the functions ber, bei, 
ker; and kei; are given in (2). 


A=a+ié 


B=y7+%6 


and substitute, together with the defining Equations [7] into 
Equation [6]. Thus we can separate the real and imaginary 
parts of the complex function W and obtain 


F(p) 
G(p) 


a ber: p + B bei; p + y ker; p + 6 kei: p 
B ben p — a bei: p + 6 ker: p — ¥ kei: p 


At the fluid-airframe interface there is no relative motion 
between the fluid and the airframe; i.e., the fluid sticks to the 
airframe, and the velocity of the fluid at that point is equal 
to the airframe velocity. In addition, since the airframe is in 
a state of pure rotation, the fluid velocity at the missile roll 
axis must vanish. Thus if the airframe is excited harmoni- 
cally, the fluid velocity function must satisfy the following 


conditions 
o(A, t) = vo sin wt 


~~ 
S 
~ 
Ss 
ll 
ay 


Eprtor’s Note: The Technical Notes and Technical Comments sections of Jer PROPULSION are open to short manuscripts describing new 
developments or offering comments on papers previously published. Such manuscripts are published without editorial review, usually 
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where J is the effective moment of inertia of the fluid, and k 
represents the damping coefficient. Substituting Equation 
[3] together with Equation [10] into Equation [11] we obtain 


I 2 
NX and performing the required differentiations and _ substi- 
\ | tutions in Equations [12] we obtain the results 
N 

(ber: A — bei A bers A) 
0.005 fat A(ber;2? A + bey? A) [13] 

0.004 
0.003 Aber; Abér; A + bei; A bei; 

For sufficiently large values of we can replace the Bessel 
V2.0 functions by their asymptotic approximations (see (2)). In 
Vy this case formulae [13 and 14] reduce to 
Fiz. 1 Dimensionless plot: Moment of inertia vs. frequency ; 
parameter 5 
/2 /9 
“7 a 2 2 
Conclusions 


4 
Fig. 2. Dimensionless plot: Damping vs. frequency parameter 


Consequently, if we substitute Equation [8] into Equation 
[3] and impose condition [9] upon the result, we shall find 
that 


—v, bei; A 


(ber;? A + bei,? A) 
Y= 6 = 0) Ps 
(ber;? AX + bei;? A) 
= (ber; ber p + bei; A bei; p) 
(ber;? + bei,? d) 


fu 
= = (:)| 
Op \p/ Jp=n 


torque = 

q Vp kup. 

I cos wt + sin wt 
a a 


The degree of propellant participation in the rolling motion 
of the airframe is frequency dependent. Plots of Equations 
[13 and 14] are shown in Figs. 1 and 2, respectively. It is 
seen that the higher the frequency, the lower the effective 
propellant moment of inertia, whereas damping increases 
with rolling frequency. These remarks apply to relatively 
smooth tank walls. Appropriate adjustments in the applica- 
tion of this analysis must be made for tanks with wall pro- 
trusions. 

The author wishes to thank J. C. Goodwyn, Chief, Systems 
Analysis Section, The Martin Company, Denver Division, for 
encouragement and permission to publish this note. 
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Optical Determination of Orientation 
and Position Near a Planet! 
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Introduction 


"WHE visible horizon is the primary reference for establish- 

ing the vertical in marine navigation. It is natural to 
consider the analogous operation of establishing the direction 
to the center of a planet from a satellite or space ship by ob- 
servation on the visible disk of the planet. This has, in fact, 
been proposed in astronautical literature (for example, in (1)*). 


1 Presented at the ARS Semi-Annual Meeting, Los Angeles, 
Calif., June 9-12, 1958. 


2 Chief, Astronautical Sciences. Member ARS. 
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Exactly the same kind of device required to find the center 
of the disk would be expected to give the extent of the disk— 
its subtended angle—as well. Such a measurement has long 
been proposed (2) as a means of determining position, either 
relative to a close planet or to the sun. 

The vertical, or direction to the planet’s center, may be re- 
quired for various reasons. In the case of a space vehicle ap- 
proaching a planet and desiring to orbit or land, it is an im- 
portant navigational orientation reference. It is important 
for the same reason in any ferry operation between the planet 
and an orbiting space station. Fora satellite, the navigational 
problem may vanish, but the problem of a primary orientation 
reference system for instrumentation and communication 
then arises. In each case, the usefulness of the extent of the 
disk as a measure of position is clear. 

The purpose of this note is to describe a rationale for de- 
termining orientation and position relative to a planet by 
observation of its disk. Specifically, the following measure- 
ment capability is assumed to exist aboard the vehicle. A 
physical device examines the intersection of the planet’s face 
with the family of planes containing a reference axis on the 
vehicle—the latter forming an axis of rotation for an 
azimuthal mode of examination. In each plane, the angle can 
be measured between the reference axis and the edge of the 
disk. 

Equations are developed which convert this measured 
angle into orientation and position information, at least in 
principle. They are derived purely on a geometric basis and 
are independent of the physical device used provided the sub- 
tended angle is available in each azimuth direction. 

The present note is not concerned with the details of such 
a device, and thereby ignores some formidable practical prob- 
lems. For example, if the edge of the disk is signaled by the 
albedo change between the planet and its space background, 
how is the case included where the planet has the appearance 
of a half-moon to the observer? It also leaves aside the prob- 
lem of choosing a working spectrum for the device, which 
may be conceived to operate in the thermal, visible or even 
other portions of the electromagnetic radiation spectrum. 
Nevertheless, it is expected that the basic results developed 
here can be modified suitably to apply to cases where Jess in- 
formation is present than we have assumed, but where ad- 
vantage can be taken of some of the operating details of 
actual physical devices. 

Nearly as important as the basic relationships used to infer 
orientation and position is the question of errors in the 
method. Some error sources are mentioned briefly. Most of 
these are stochastic; some depend on the particular mecha- 
nization used. One, however; is purely geometric and can be 
treated in detail, namely the planetary oblateness. It will 
be shown that this has a curious effect, leading to systematic 
errors in both orientation and position. at 

Basic Coordinate Systems = 

Choose axes x and y in the plane of the equator of the body 
to be observed and z along the polar axis so that xyz is a right- 
hand orthogonal system. Let the observer, i.e., the satellite or 
space ship, be situated at a radial distance p from the center 
of the body and at a geocentric latitude \.. Until the motion 
of the observer becomes a factor, there is no loss in generality 
in assuming that the observer is in the yz plane. Consider a 
set of reference axes x’y’z’ at the observer such that y’ is along 
the outward-pointing geocentric radius, and 2’ is parallel to z. 
The geometry of the solution is shown in Fig. 1. 

Two lines are of interest as seen from x’y’z’. One of these 
is the axis S, of the vehicle whose orientation, nominally ver- 
tical, it is desired to measure. The other is an arbitrary ray 
S: from the edge of the planet’s disk to the observer. Let the 
former be characterized by angles a, B and the latter by 0, 
as shown in Fig. 2. 

The angle @ between S, and Sz is the subtended angle seen 
relative to the vehicle’s nominal vertical at the true azimuth 
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S, (NOMINAL VEHICLE 
VERTICAL) 


y 
S2 


RAY FROM EDGE 
OF PLANET 


OBSERVER 


Fig. 1. The observer in relation to the observed planet 


Fig. 2. Basic angle variables relative to x’y'’z’ system 


¢. It is not difficult to see that it is given by 
cos @.= cos 6 cos a + sin 6 sin a cos (@ — B)... [1] 


If we suppose that the vehicle vertical has been established 
almost correctly so that a isa small angle, @ is approximatel 


@® = 6 — acos(@ — B)..... [2] 


It is not unreasonable to require that the axis of the sensing 
device be oriented approximately vertically—this is fairly 
easy to accomplish—so that attention can be confined here 
to the small a case. 


Orientation and Position by Averaging 


If the planet is a sphere with radius a, it is obvious that 
6 = aresin (a/p)..... 


should be used for 6 in Equation [2]. 

Now @® given by Equation [2] is the output of a physical 
device, and is a function of the azimuth angle g. Suppose 
that all values of g can be examined rapidly enough that 
changes in a, B and vehicle position (although the latter is 
not important until an oblate planet is considered) are in- 
significant. Then the physical output @ can be operated on 
in various ways; specifically, it can be averaged with respect 
to ¢ in various ways. Suppose that the following quantitics 
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Its importance depends on both the instrument design which 


can be formed by computation or filtering aboard the vehicle 
determines the g-frequency and on the attitude and transla- 


= , @.; = /2 ig <i aw tional vehicle velocities. Unfortunately, it may not be able 
to make the effects arbitrarily small by increasing ¢g-frequency, 
av av because the sensitivity and some instrument errors may 
Lo <¢s = <¢s 1/2 change adversely with increasing frequency. This source 
4 : may be troublesome if one is trying to infer rates of change of 

@u = Se (4) radius or attitude with the present method. 
Finally there is a class of errors which arise from the 
(The subscripts in the last four cases represent quadrants of ¢ geometry of the planet’s disk and the physics of its at- 
vues.) mosphere. It is clear that the edge of the disk will have a lo- 
Performing this averaging in Equation [2] and using Equa- cation error in the presence of clouds or mountains, at least in 
tion [3], it is easv to show that the visible spectrum. Atmospheric refraction and scattering 
oe also may be troublesome. One error of this class is particu- 
leah ee [5] larly interesting, that from any oblateness of the planet. It 


has systematic, pernicious components. These effects of 
asin = — @r) 4  oblateness are the concern of the remainder of the paper. 
€ Oblateness Effects 


In terms of the coordinate system xyz introduced previously, 
suppose that the body to be observed has the equation 


@ cos B — 


(i — + — [7] 
The distance p from the center of the planet to the observer is : : , 
g.ven directly by Equation [5], so @ is a measure of relative where a is the equatorial radius, and ¢& = 1 — (polar ra- 
position. For small a, the quantity a sin B is just the angle dius)?/a?._ That is, it isan oblate spheroid with polar axis Z 
0! rotation of S, about the —2’ axis, and a cos B about the and eccentricity €. (Note that this is not ellipticity.) Using 
z’ axis, so that the two quantities are a measure of the the transformation between xyz and x’y’2’ coordinates im- 
igular misorientation of S, relative to the y’ axis. plied by the previous definitions, the equation of the ellipsoid 
as seen from the x’y’z’ system has the equation 


(1 — &)(x’? — a?) + (1 — € cos? A)(y’ + p)? + 2e? sin A cos Az’(y’ + p) + (1 — esin? A)2” = 0........... [8] 

It is important to note in this connection that one really Now if r represents distance from the —y’ axis in any plane 

may desire these misorientation angles about body-fixed defined by ¢ = const (see Fig. 1), x’ = r cos g, and 2’ = —r 

horizontal axes rather than x’ and 2’, for if attitude control is sin g. Therefore, the cross section of the ellipse in this plane 

= — has the equation 

[1 — €%(1 — cos? AX sin? ¢) |r? — 2e? sin X cos A sin gr(y’ + p) + (1 — € cos? A)(y’ + p)? = (1 — €)a?.......... (9]’ 

At the same time, 6 is defined as the angle between the —y 

desired the corrective action would take place about body- axis and the ray r = —y’ tan @ which is just tangent to the 
fixed axes. It can be shown, though, that this distinction does ellipse given by Equation [9]. 

not vitiate the present results if the body-fixed horizontal A straightforward way to find @ is to substitute r = —y’ 

axes are rotating about S, slowly relative to the speed at tan @ into [9], and then to require that the resulting quadratic 


which the averaging over ¢ is accomplished. Although not in y’ have only a single root. In this way, @ is found to satisfy 


crucial in the simple picture presented here, the distinction is 
nevertheless of genuine interest and will be treated elsewhere. tan? + 26 tan 0 + ¢ 0 [10] 


where 
a2 a? 
A =]1—-——+ e&{— — cos?A 1 — €%(1 — cos?) sin’g) | —e* sin? cos? sin? g.............. . (11) 
p? p? 
: 
B= —e(1 — — sin AcosAsin¢g........ [12] 
Errors in the Method 
There are several obi ious sources of ¢ rror in the me thod de 
scribed above. One is the instrument noise, depending on the p? 


details of the physical device observing the planet. It may be 


stochastic, or it may have systematic components, such as We confine ourself to the case where ée? is small enough that 


e‘ and higher order terms can be discarded. In this case, the 
solution [10] takes the form 


= sie 


a ae’ a 


7 The upper algebraic sign must be chosen if @ is to remain a 


would be the case if it habitually triggered too early or too late positive vertex semi-angle, the physical equivalent of the op- 
n passing from regions of high to low albedo or vice versa. posite sign being realized for g + 7. Using Equation [3] for 
[his source would also encompass imperfections in the @, and continuing the same order of approximation in e’, 
iveraging process leading to Equation [4]. Equation [14] gives 

Another potential source of error is attitude or radial mo- a 
tion of the vehicle while g ranges over the values [0, 27]. ATE [15] 
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chanics and may be of interest because they show simple, 


cos \ sin da [16] exact solutions to apparently very difficult equations. 

se 2p? V p? — a? By the plane motion of a mass point the accelerations a, an‘ 
7 a, in the r- and ¢-directions, respectively, are given by 


Using Equation [16] for the oblateness error and perform- 


the averaging of Equation [4], we obtain Gy — MW 
eav/ p? — a? [1 2asinXcosh 49 


2 2a sin cos a’ 
2p? 2 rv p? -a 
b#2 = 6@u = 6@ If the origin of the coordinate system is placed in the ma:-s 
[18] center of the Earth, the satellite will move under the action 
se P, acting in the direction toward the origin, of tle 
From Equation [5], 6p/p = p? — a? 6@/a, whence the 
position error arising from oblateness is mag ¢ 
_ where m is the mass of the satellite. The equations of motion 


~ may now be written 


27,2 
2p 
‘ 
r dt 


These results are easily interpreted. Esuation [21] shows 
that there is no angular error about the z’ axis from oblateness, 
as would be expected from symmetry considerations. Equa- 


where F and G are the components of the frictional force in 
the r- and ¢g-directions, respectively. By neglecting friction 


tion [20] shows that there is, however, an error about the 2’ these equations represent the equations of motion of planets 
(east to west) axis. It is positive about x’ when the observer and give ellipses, parabolas and hyperbolas as only possible 
is in the northern hemisphere, negative in the southern orbits. 


It is at present very difficult to give any reliable expression 
for the frictional forces resulting from the motion of a satel- 
would vanish in the average. The technique would not neces- lite; one has to make more or less reliable assumptions. If 
sarily work for a nonorbiting vehicle, of course. On the other we suppose the friction to be caused by some sort of gas wit!) 
hand, the error in radial distance given by Equation [19] is an very low density, it seems justifiable to assume the frictional 


hemisphere. Thus if a sin 8 obtained from [6] were further 
averaged over an entire orbital revolution, this latter error 


essential one, and cannot be removed by orbital averaging. force to be proportional to some power of the satellite’s ve- 
To give some idea of the magnitudes of these errors, con- locity, and also to assume the density of the gas to be decreas- 

sider observations on the Earth from a distance p = 1/3a ~ ing with increasing distance from the origin. Two different 

11,000 km and a Jatitude 45 deg N. We have dp ~ 12kmand assumptions for the frictional force made in accordance with 

5(a sin B) ~ 0.547 mr = 2’.. This angular error is a maxi- this general point of view will now be treated, and exact solu- 

mum with respect to \ and decreases with increasing p/a. tions of Equations [4 and 5] will be established. 

One can ascertain from Equation [19] by properly grouping ee 


terms, that the relative error 6p/p is positive for all p and ); 
also, that 6p/p— (e cos *A/4) as p/a—> ©. It is assumed that the frictional force is proportional to the 
second power of the satellite’s velocity and decreases as 1/r as 
References rinereases. In this way the variation of the density is taken 


1 Stroud, W. G. and Nordberg, W., “Meteorological Measure- into account to some extent. Equations [4 and 5] may then 


ments from a Satellite Vehicle,’ in “Scientific Uses of Earth Satel- be written 
lites,”’ ed. J. A. van Allen, Univ. of Michigan Press, Ann Arbor, > j2 
2 Esnault, R. “Pelterie, L’Astronatique Complement,” Paris, 
1935, p. 87. id ry? 
r dt r 
where & is a dimensionless “coefficient of friction.’”’ In those 
Motion of a Satellite With Friction -ases where the orbit of the satellite is close to a circle, F may 


be assumed to be negligible in comparison with G, without 


introducing an error of any importance. The last term in 


Norwegian Institute of Technology, Trondheim, Norway Equation [6] may thus be neglected. Remembering that 


N THE following, two examples of a satellite’s motion 7 _ 7 
under the action of frictional forces will be presented. 
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Equation [7] may be solved by simple integration rendering 
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been calculated and is shown in Fig. 1. 
chosen unrealistically great to make the characteristic features 
of the curve come out more distinctly. 
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Fig. 1 Satellite orbit in case 1 with an unrealistically great value 


of k 


where C; is a constant of integration. The time ¢ is now re- 
pluced by ¢g as the independent variable of Equation [6] by 
means of the following expressions 


d? d? 4 d 1 


@ 


where the last expression is obtained from Equation [9]. In 
this way Equation [6] takes the form 


2 9 
dg? r dg Ci? 


This is a nonlinear differential equation with variable co- 
efficients, which, however, may easily be transformed into the 
following linear equation 


d? d 
+§ 
dg dg 
with the solution 
k ke 
r 4 4 


where C2. and C3 are constants of integration to be determined 
by the initial conditions. Assuming these to be 


Vo 
¥ ip 
and introducing the following numerical values 
16 2gk? 
yR? k= 0.02 


25 ro 


where R = radius of the Earth, the orbit of the satellite has 


ro 


The value of k is 


Case 2 
‘The assumptions made in the first case are not satisfactory, 


because F has been neglected, and furthermore the variation 


+ 2k?) 


1 1 
= cos — KB? sin g + K*B2{ — Ci(k — ¢) 


of the density seems to be too small. If the motion of the 
satellite is regarded as the motion of a body in a gas with 
very Jow density, a Reynolds number may be written 

Re = si 
which, in spite of a very high velocity v, may be assumed to be 
small because of the extremely low value of p. This makes it 
tempting to assume the frictional force to be directly propor- 
tional to the velocity; and if, in addition, the variation of 
density is taken into account by assuming F and G decreas- 
ing with the second power of r, Equations [4 and 5] may be 
written 


r 
r? r2 
ld ra 
2 


where c may be regarded as a dimensional coefficient of fric- 
tion. These equations may be solved in exactly the same 
way as those of the first case. By means of Equation [8] the 
solution of Equation [14] is found to be 


.. [15] 


from which the following relation may be derived 


lr 
= (cw + 2rw? 
dp 


With this expression for Equation [13] may be transformed 


into 
d*r 2 2 Ar? 0 (16), 
where the same procedure as that in the first case has been | 
used. Introduction of € = 1/r will further transform this 
equation into 
ad? 
+ = ——............... 
dg? (C, — cg)? 


the solution of which is easily found to be 


=, 
= sin + Cs cos — [cos + sin si(m)].. [18] 


Ne 2ke 


where Ci and si denote the cosine and sine integrals, respec- 
tively, and = (Ci/c — ¢). The constants of integration C), 
C’, and C; are now determined by the initial conditions of the 
first case. Furthermore the new nondimensional coefficient of - 
friction « defined by 


c= 
K 


and the abbreviation B? = 2gR?/vp?rp are introduced, where- 
upon the equation of the orbit may be written __ 


] cos (kK — + [si(x) — si(x — sin(k —g)}......... [19] 
Also, in this case the orbits turn out to be elliptical spirals. 
Examples for different values of 8 and « are given in Figs. 2 
and 3. 

It may be worth noticing that for the same initial value of — 
the frictional force, the spirals of the two cases (Fig. 3) are 
not distinguishable in the seale to which the figure is drawn. 
Furthermore, it may be mentioned that the second case per- 
mits an easy determination of the number of revolutions 
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Rg 3 Satellite’s orbit in case 2. 


ae 


Elliptic orbit 
without Friction 


i; — 


ig. 2 Satellite’s orbit in case 2. 
greatest in perigee 


a 500 


Elliptic orbit 
without friction 


Height loss in percentage 
greatest in apogee 


which the satellite will perform before reaching the origin. 
Once x is known, this number is x/27. 


_ It would have been interesting if these easy calculations had 


permitted a somewhat sensible determination of the lifetime 
of Sputnik I. From Equation [15] the following expression 
for the period is found 

2x 

0 (C, cy) 
where r = r(g) has to be inserted from Equation [19]. Re- 
placing c by c = c,voro will make 7 a function of c,, where c; is 
very small, and 7 may be expanded into a power series in ¢ 
Or 
oC; lo 1! 
where 7» denotes the period with no friction. 
the following expression for the time-lag 


T= 


This will give 


Ar 


Using the data for Sputnik I,? this time-lag is Ar/7 = 

0.146666/5760. The coefficient 8 may be calculated, giving 

B? = 1.89888 which further leads to 

or 

Oc; 

2 “Radio Observations of the Russian Satellite,” Nature, vol. 2» 
Nov. 1957, pp. 879-883. 


= 14.079 
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+ Satellite 


Height loss in percentage 


Physical Propulsion Systems! 


_ is the most serious problem in design, whereas the energy con- 
tent of the fuel, so to speak, is free. 


system in which the functions of supplying mass and energy 


4 partment of the Navy, Contract NOrd 7386. 


whereupon k = 1/c,; = 552,900 is obtained from Equation 
[21]. The values of x used for plotting the spirals in the dia- 
grams shown are extremely small compared to this value. 
The given value of x leads to 8800 revolutions before Sputnik I 
has reached the center of the Earth, but this is of no impor- 
tance in connection with the lifetime of the satellite. Because 
of the extremely high value of x, the Ci- and st-functions m:y 
be replaced by the asymptotic expansions, whereupon Equa- 
tion [19] takes the form 


1 1 
K 2 K 


(n = integer) 


1 1 \ 
= (1-5 (1-™*) [3] 
r 2 2 K 


this equation makes it possible to calculate r after n revo u- 
tions, or it may be used to obtain an approximate value or 
the number of revolutions after which r will have a given valie. 
When the satellite starts to dive into the atmosphere, the 
present calculations cease to be valid, but from then on it is 
only a question of a few days until a satellite without any «ie- 
celerating devices will burn out. Counting the lifetime of 
Sputnik I to the time when it starts diving into the atmos- 
phere, and putting the thickness of the atmosphere at 50 km, 
Equation [23] gives n = 1040 and a “‘lifetime”’ 7'59 = 70 davs. 
With an atmospheric thickness of 28 km, n becomes 1200, and 
the lifetime 73 = 80 days. If these calculations are realistic, 
one may expect that Sputnik I dived into the atmosphvre 
and burnt out shortly after the period Dec. 29, 1957 to Jan. 


On the Optimization of 


For ¢ = 


ROBERT R. NEWTON? 


The Johns Hopkins University, Silver Spring, Md. 


Introduction 


N A rocket with chemical propulsion, the fuel performs two 
quite distinct functions: It supplies the mass which car- 
ries the momentum of the jet, and at the same time it supplies 
the kinetic energy which is involved in giving this momentum 
to the jet. Fora single stage with a final mass of unity, initial 
mass V and an exhaust velocity c, the velocity in the absence 
of external forces is 


In all experience to date, it has been desirable to make c as 
large as possible, in order to make M as small as possible. 
This situation arises because obtaining the values of M needed 


By a physical propulsion system for a rocket, we mean a 


are separated, and are performed by different components of 
the system. In such systems, we are no longer necessarily 
limited to values of c of the order of 104 ft/sec, so that obtain- 
ing the values of M needed to achieve a high velocity may no 
longer be a serious problem. Instead, the limitation may be 
in either the power or the total energy which can be delivered 
by the source of energy. 

These considerations suggest the following problem: 
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. . . . . . . . 
regarded as a completely free parameter, is there any choice is discriminant of this quadratic gives the condition that 


of c which leads to minimum energy or power, for a fixed value 
of »? We shall analyze this problem briefly in two special 
cases, first with no external force and second with a constant 
gravitational field. The analysis also applies, of course, to 
chemical rockets, if we regard either energy content or de- 
livered power as more important than the mass of fuel. 


Optimization With No External Forces 


In this case, Equation [1] gives the final velocity, regard- 
less of the thrust or power level at which the motor operates. 
We are therefore concerned only with minimizing the total 
encrgy demanded of the energy source. The energy source, 
carried on the rigid framework of the rocket, must give to the 
exliaust, of mass M — 1, the velocity c relative to that frame- 
work. The energy demand W is therefore equal to the 
kinetie energy of the exhaust relative to the rocket, and is 


[2] 
Suostituting for c from Equation [1] 
W = — 1)/(In M)?].............. [3] 
W. wish to minimize W, for fixed v, which means minimizing 
th: quantity in brackets. 
sefore doing this, we remark that the useful accomplish- 
ment of the rocket is to give unit mass a velocity v in a “fixed” 
coordinate system. It is therefore proper to define the energy 
output as 4v?, and to define the efficiency as 


E = (30?)/W = (In M)?/((M 1)..... 
The efficiency is zero when M = 1 and ©; in the inter- 
val (0, ©), # has a single maximum, at 
[5] 
From this, we find 


This case approximates the launching of an escape vehicle 
from a satellite. For this application, v is about 104 ft/sec, so 
the optimum exhaust velocity is about 6300 ft/sec. 

In this analysis, the energy source is counted as part of the 
payload, and the kinetie energy given it as part of the useful 
output. In a more careful analysis, one would probably 
modify accordingly the definition of the output energy, but to 
do so would require knowing the relation between the mass 
and available energy of the energy source. Studying this re- 
lation is outside the scope of this note. 

Optimization in a Gravity Field 

If the rocket is required to achieve the velocity v while 
climbing vertically against gravity, the velocity is now given 
by 


Clearly, the maximum efficiency is achieved with t = 0; Equa- 
tion [7] is then the same as Equation [1], and the earlier 
analysis applies, if maximum efficiency is the criterion of 
optimization. 

However, achieving maximum efficiency in this case in- 
volves operation at an infinite power level. It is clearly of 
interest to investigate, instead, minimizing the power demand 
on the energy source. The power demand is 

[8] 
where 7 is the thrust, and m is the mass flow rate. Assuming 
that P, T and c are constant, the minimum allowable value of 
T is Mg, since the rocket must be able to lift its initial weight. 
Combining Equations [7 and 8] with 7 = Mg, and using the 
fact that ¢ = (M — 1)/m = (M — 1)c?/2P, we find that c 
must be a solution of the quadratic 


— 1)/2PJce? -clnM+v=0........ [9] 


1958 


P(In M)2/((M — 1) > 2gv..... 


Since 290 is fixed for a given application, the minimum value 0 
P occurs when the equality holds, and also when (In .1/)2/ 
(VW — 1) isa maximum. So, again, we have the optimum value 
of M to be 4.9214. Other quantities of interest are 


> = 1.2550 m = 3.9214(g/v) 
= (v/g) P = 3.0883gv 


™ Oo 


In units of horsepower per lb, using the sea level value of g, 
the minimum power is . 


P = 0.005615 


For v = 25,000 ft/sec, this is 140.4 hp per Ib of the final 
weight. 

The value of t shows that the gravity loss equals the final | 
velocity. Further, the value of c is just twice the value which 
gives maximum efficiency, so that the efficiency is only one 
fourth of the maximum or about 16.2 per cent. The distance — 
of powered flight is 


s = 0.24499v?/g 


For v as before, ¢ = 31,375 ft/sec; t = 776 sec, and s = 4.7552 
X 10° ft, about 900 miles. 


Summary and Comments 


Regarding the exhaust velocity ¢ as a free parameter, we 
have found the values of ¢ which give maximum efficiency | 
and those which require minimum power in climbing against — 
gravity, for single-stage motors. We have found that the _ 
mass ratio should be 4.9214 for all cases considered, and that | 
c should be proportional to the desired velocity. The value of 
c for the minimum power when climbing is twice the value | 
for maximum efficiency, and the efficiency is one fourth of the 
maximum. 

In launching a satellite, the minimum power is about 140 hp 
per lb, and the corresponding ¢ is about 31,400 ft/sec. In 
launching an escape vehicle from a satellite, the optimum c is 
about 6300 ft/sec, and the minimum power is quite low (zero 
to the accuracy of the present treatment). Thus, we have an 
ironic situation: The only energy sources which, to date, can 
deliver the necessary power for launching a satellite are— 
chemical, and they have ¢ values far from optimum. How- 
ever, for an escape vehicle, where power is not a limitation, 
chemical sources have about the optimum parameters. 

One cannot conclude from this that chemical rockets should 
be used for an escape vehicle, because there are important 
considerations other than maximum efficiency. However, 
since chemical rockets do operate at about the maximum 
possible efficiency for this application, they need to be 
considered. 


Measurement of Satellite Erosion Rates 
by the Backscattering of Beta-Rays 


R. C. GOETTELMAN! 


Stanford Research Institute, Menlo Park, Calif. 


A sensitive method of measuring thickness changes, 
using beta-ray backscattering phenomena, is described. 
Experimental data and theory are used to show that the 
backscattering method is capable of detection of thickness 


changes characteristic of anticipated satellite erosion rates. 
The method is especially attractive since it does not 
physically alter the ‘skin 


material surface, and therefore 
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should be representative of the actual satellite erosion rate. 
A null circuit is proposed which eliminates counting errors 
due to cosmic activity and provides a maximum sensitivity 


thickness. 


CCURATE measurement of surface erosion on orbiting 
vehicles is important for many reasons. In addition to 
its scientific value in adding knowledge about the universe, 

it has many practical applications. 

Absorptivity and emissivity are directly dependent on the 
surface condition of a metal. The control of temperature 
within an instrumented or inhabited satellite or space ship 
is completely dependent on radiation losses and gains. Thus, 
in the design of a controlled temperature vehicle, the amount 
of surface erosion becomes important if long times of flight 
are expected. 

At altitudes used for the employment of satellite systems 
we are concerned with free molecule flow, or superaero- 
dynamics. In this realm two drag mechanisms operate: 
Newtonian (inelastic) reflections and specular (elastic) re- 
flections of the gas molecules. Specular reflections are those 
in which the incident stream of particles is elastically reflected 
from the exposed body surface, and are directly dependent 
upon the “smoothness” of the reflecting surface, i.e., the 
roughness height (1).2- Although only 3 to 10 per cent of the 
total molecules are specularly reflected (2) such reflections 
are important over long time intervals and when accurate 
valculations for re-entry vehicles are required. 

Several erosion mechanisms (3) and erosion rates (3, 4) 
have been suggested by researchers. Singer suggests a maxi- 
mum erosion rate of 10u (microns) per 1600 sec, while 
Whipple favors a rate of 2 X gm/cm?/sec (or 107® 

per 1600 see for an aluminum target). This paper describes a 
oy  beta-ray scattering technique which can be used to measure 
surface erosion rates between these two values. 


Introduction 


Proposed Erosion Measuring Method 


Several methods of erosion rate measurement have been 
suggested in the literature (3,5). One of these methods (3) 
aA utilizes the number of beta particles that are backscattered 

from a thin target as a direct indication of the material thick- 
ness. By ing the change in backscattered 


should to the rate which a physical 
s process is removing material from the target surface, e.g., as 

in micrometeorite erosion. This paper describes the applica- 
___ bility of this process to satellite operation. 

_ The principle advantages of the use of a beta-ray backscat- 
allele technique are its sensitivity, its stability and the fact 
that it operates without phy sically altering the measured 
surface. 

a The purely theoretical treatment of backscattering is quite 
- complicated and not completely understood. Fortunately, 
thickness measurement by beta-ray backscattering has been 
a _ employed sufficiently in industry and medicine for evolution 
_ of several useful empirical formulas. It is hoped that, by 

treating the results of some backscattering experiments with 
these formulas, enough confidence will be developed to show 
_ that we may extrapolate this technique for use between the 
erosion limits cited above. 


Backscattering Experiments 


A technique was developed at Stanford Research Institute, 
Menlo Park, California (6) for measuring thicknesses of ice 
_ by using the backscattering of beta-rays. It was shown that 
for thin targets the number of bac -kscattered beta particles is 
proportional to the target thickness. 

The essential elements of this system were a beam of beta- 


* Numbers in parentheses indicate References at end of paper. 


which to a 6.009-mil change in aluminum 


Component arrangement for beta-ray scattering expe: i- 
(1) Collimator and source, (2) scattering material (tar- 
get), (3) detector, (4) scaling unit. 


Fig. 1 
ments: 


/ DETECTOR 


| 


LraRcet Lan 


/ LUCITE CAP 
+ IRIS 


| j 
/ BETA EMITTING 
77.25 
| LUCITE 
= 
| REAR PROTECTION 
PLATE 
| 
STAINLESS STEER 
|} source HANDLE 
Ls 
Fig. 2 Detail arrangement schematic for beta-ray scattering 


experiments 


rays, a target area and a radiation detector. These are shown 
photographi sally in Fig. 1 and schematically in Fig. 2 

The beta source was a Tracerlab RA-1A Medical Applica- 
tor. The applicator was a beta-emitting plaque which con- 
tained strontium 90 in equilibrium with yttrium 90. Stron- 
tium 90 disintegrates with a half-life of 20 years, emitting par- 
ticles of 0.65 mev to form yttrium 90. Yttrium 90 has a ha!f- 
life of 62 hr and emits beta particles with a maximum enerey 
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of 2.16 mev to form the stable isotope, zirconium 90. 

Fig. 3 shows that the number of backscattered beta par- 
ticles (counts/min) increases with the target thickness. 
The counts per min tend to increase until a state of ‘“satura- 
tion’ predominates. At this point a further increase of 
target thickness does not noticeably increase the number of 
backscattered particles. 

The two curves in Fig. 3 show the results of two different 
collimator “‘iris’”’ openings. The point of saturation at con- 
st::nt source angle and distance, constant detector angle and 
distance, is nearly identical (about 56 mil for polyethylene) 
for the 3- and 3-in. iris openings. The only significant dif- 
ference between the two curves is the improved sensitivity 
(change in counts/min per change in thickness) shown by 
th» larger iris opening. 

Che importance of the saturation thickness is that it sets 
a ‘thickness limit” for any satellite skin in which we wish 
to use backscattering techniques. It is also reasonable to 
as-ume that target material, target thickness and beta parti- 
ch energy will determine the amount of incident radiation 
sc.\ttered for detection. 


Scattering Theory 


The scattering of swift electrons occurs in an absorber both 
by nuclear deflections and by inelastic collisions with atomic 
el-ctrons” (7). At electron energies of 2 to 3 mev, the mini- 
mim distance of approach of the beta particles during colli- 
sion with the target atoms is significantly greater than the 
ra‘lius of the nucleus. Thus, the only force acting is the in- 
verse square coulomb force, so we may apply classical Ruther- 
ford scattering to our problem. 

[t can therefore be shown (7) that 


M,\?2 


E, M, 
where 
Z = atomic number of the target material ‘ 
E = energy of beta particles 
7, = mass of the beta particles a 
= mass of the target material 
o = Rutherford scattering cross section, em?/atom, and 
n = nmo(NAS).......... 
where 
n = number of particles scattered 
ny = number of particles incident on the target 
V = number of atoms, em? 


AS = target thickness, em 


We see that the assumption about the dependence of scat- 
tering on target material, target thickness and beta energy is 
supported by classical theory. on 


Empirical Formulas 7 


Theoretically (8) saturation is reached if the target thick- 
ness is equal to half the range (2) of the beta particles. Prac- | 
tically, the maximum backscattering occurs at a thickness 
d = 0.2R. Saturation values are found to be independent 
of the particle energy between 0.3 and 2.3 mev. This can be 
explained by the fact that, for increasing energy, the de- 
crease in scattering cross section, Equation [1], of the nuclei 
of the target material is compensated for by the increase in 
the number of nuclei encountered by the beta particles (be- 
cause of the increase in range). 

For energetic electrons the range is found empirically (9) 


to be 


?% = R times target density, gm/cm? 
E energy of beta particles, mev 
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Fig.3 Dependence of beta-ray backscattering on total thickness 
of polyethylene targets 


FO 


Fig.4 Dependence of beta-ray backscattering on target material 


Caleulations 


~ Comparing this theory with the results of the beta scatter- 
ing experiments for polyethylene targets (6) we find: If 
E avg = 1.3 mev, i.e., EF avg equals 0.6 Emax where Emax = 
2.16 mev, then Ry = 0.586 gm/cm?, and R = R,/density = 
247 mil. Fig. 3 shows that saturation occurs at about 56 
mil. Thus we see 56/247 = 0.22, or more properly, satura- 
tion thickness 


which is in rather good agreement with the value above, i.e., 
d = 0.2R. 

We may now with some confidence predict the saturation 
thickness (d) of aluminum and magnesium, two commonly 
used satellite materials. Using Equations [3 and 4], and E 
avg = 1.3 mev, we find the following results: 


Saturation 
Material Density Range R _ thickness d 
gm/cm$ mil mil 
3 polyethylene 0.93 247 56 
aluminum 86 19 
magnesium 1.8 128 28.5 


Thus, if the strontium 90, yttrium 90 source were used, the 
maximum skin thickness would be restricted to about 30 mil 
for magnesium. Skin thickness in excess of the saturation 
thickness would not be prohibited however, if “windows,” | 
i.e., areas which had been machined to the proper depth, | a 
were provided for the beta scattering targets. *] 

Fig. 4 shows the dependence of backscattering on the target 
density. The two curves indicate that aluminum is about 
twice as reactive (counts/min) as polyethylene at equal 
thicknesses. 


From Equation [2] we could expect magnesium | 
to fall somewhere between these two curves. The slope of | 


the aluminum-polyethylene curve suggests that, despite the 
high count rate, the counting tube was not overloaded in 
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Fig. 5 Proposed circuit arrangement for increasing measuring 
sensitivity and excluding cosmic background 


SATELLITE SURFACE 


this region, an important fact when considering that the ex- 
perimental saturation point might have been due to the tube 
and not the scattering mechanism being saturated. 

At the speculated upper limit erosion rate of 10 u per 1600 
see (3), or about 0.9 mil/hr, we see in Fig. 3 for the 4-in. 
iris curve, that the counting rate would have changed by 
about 5 X 10? counts/min for a polyethylene target whose 
thickness had changed by this amount. 

We consider only the upper portion of the curve since the 
satellite skin would naturally be very near the saturation 
thickness for the particular material and beta-ray energies 
at the start of our observations. Since the aluminum re- 
sponse is about twice that of equal thicknesses of polyethylene, 
this erosion rate would produce a change of about 1 x 10° 
counts/min, and for a magnesium surface possibly 8 X 10? 
counts/min. If we postulate that the erosion rate is 0.01 
to 0.001 of the upper limit of 0.9 mil/hr, we would expect a 
response change of only 1 to 10 counts/min for aluminum for 
each hour interval—a change too small to detect against 
natural cosmic background. 

The obvious solution is less target thickness and longer 
intervals between samplings. Thus, a 4-mil thick window in 
an aluminum satellite would give a response change of 1 X 
104 counts/min (Fig. 4) for every mil change in thickness. 
A 24 hr sample interval would show a response change of 
240 counts/min for an 0.009-mil change in surface thickness. 

At the lower erosion limit (4), which is equivalent to about 
2.5 X 10-6 mil/day for an aluminum vehicle, it would re- 
quire about 3600 days for the minimum response of a 0.009- 
mil depth change to become detectable. 

By making use of a bridge-type arrangement such as is 
shown schematically in Fig. 5, we can improve the detection 
of response change. The signal obtained with this new ar- 
rangement corresponds to thickness change rather than total 
thickness, and excludes the cosmic radiation background. 

An erosion rate of 0.009 mil/day (240 counts/min) ob- 
served by this method would require 12 years to completely 
erode away a 40-mil satellite. 


Conclusions 


Calculations based on experimental data have shown that 
beta-ray backscattering is a sensitive method for measuring 
the surface erosion of satellites by micrometeorites. Within 
the erosion rate limits envisioned by experimenters, the back- 
scattering method adequately covers this range, leaving some 
doubt as to its value only at the extreme low limit. The 
principle value of the method is its ability to measure directly 
the change in thickness of a vehicle skin without physically 
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Flashback and Blowoff Limits of 


Unpiloted Turbulent Flames’ 
JOSEPH GRUMER? 


Bureau of Mines, U. S. Department of the Interior, 
Pittsburgh, Pa. 


Introduction and Summary 


and blowoff limits of laminar flames have 
been explained previously by a concept based on bound- 
ary velocity gradients (1). This contribution has been most 
useful in treating gas appliance problems (2) but less so in ap- 
plications dealing with turbulent flames. It has been shown 
that the critical boundary velocity gradient for blowoff can 
be the same for laminar and for unpiloted pipe flow turbulent 
flames (2,3,4). On the other hand, flashback and blowoff 
gradients have been reported to be greater in turbulent than 
in laminar flow (5 to 9,11). Elucidation of the causes of 
these differences could contribute to the understanding of tur- 
bulent combustion and to the design of industrial equipment 
and engines. These problems have not yet been treated com- 
prehensively. In this note, the need for further study of 
stability limits of turbulent flames is indicated, and new data 
are presented on laminar tilted flame limits, and turbulent 
flashback and blowoff limits. Only unpiloted flames on un- 
obstructed ports are discussed. 

The mechanisms of flashback and blowoff of unpiloted turbu- 
lent flames have not yet been adequately elucidated. These 
instabilities stem from factors associated with the boundary 
laminar sublayer and with the core of the turbulent flow. 

The Flashback Problem 

1 Conceivably, in turbulent flow, one cannot distinguish 
between the flashback of tilted flames and true flashback. In 
laminar flow, the distinction is made readily (10). Conse- 
quently, a correlation should be attainable between tilted 
flame gradients in laminar flow and turbulent flashback gra- 
dients. 

Von Elbe and Mentser (10) studied tilted flames in laminar 
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flow only. Following the general lines of their theory,‘ the 
tilted flame gradient in laminar flow is obtained by means of 
Equation [1] 


R? 


TL. 
sag 
where K is the tilted flame gradient, which is a constant for the 
mixture, sec~!; gz, is the boundary gradient for the laminar 
flow prior to ignition, sec; R, cm, is the radius, and 8, is the 
burning velocity, cm/sec. In turbulent flow4 


K= 0.03328 71-7 [2] 


Here K is as defined above; S, is the turbulent burning veloc- 
itv; grr is equivalent in turbulent flow of gr, in laminar 
flow; o, em, is the major semi-axis of the tilted ellip- 
tical flame base (R is the minor semi-axis); p is the 
d-nsity, gm/cm, and 7 is the viscosity, poise. When R and 
o are large, Equation [3] holds for laminar and turbulent flow 


K OTL OFT. [3] 


The thesis that turbulent flashback is really the flashback 
o! tilted flames was tested by comparing the laminar tilted 
flame limits (10) with turbulent flashback limits (8) of hydro- 
gen-oxygen flames. As shown in Fig. 1, the laminar tilted 
flame limits (10) fall on the turbulent flashback curve (8), 
indicating that turbulent flashback can be caused by flame 
tilt. 

2 Another test of this concept was conducted employing 
slow-burning natural gas-air flames on much larger burners 
than those used in the hydrogen-oxygen tests. Measure- 
ments were made on three burners, 1.138, 1.435 and 1.914 em 
in diameter, to obtain laminar tilted flame limits. Turbulent 
flashback was measured using 300-cm long tubes of 5.08 and 
7.62 cm inside diameter. The data are given in Fig. 2.5 
There is no coincidence of laminar tilted flame limits and tur- 
bulent flashback limits, in contradiction to the results pre- 
sented in Fig. 1. A possible explanation is now offered. 

3 The turbulent flame stability gradients in Fig. 2 are 
numerically the same whether considered as tilted flame gra- 
dients or as turbulent flashback gradients. This is because 
Equation [3] applies for these large port diameters; actually 
no tilting of flames was observed before flashback. Next, 
consider that turbulent flames on the verge of flashback on 
wide diameter ports are short wide-angled flames. Photo- 
graphs of such a flame are shown in Fig. 3 (the port is 7.62 
cm, and the depth of focus of the optical system is about 3 
cm). Many flame surfaces within the body of the turbulent 
stream appear ready to flash back, whereas the flame surfaces 
near the boundary of the stream (laminar sublayer) are nearly 
vertical. According to Equation [4] of (5), the sublayer 
thickness / should be about 0.06 em 

66R 

The photograph in Fig. 3 strongly suggests the possibility 
that on reducing the flow, the flame surfaces leading the flash- 
hack will not all be contained in the laminar boundary layer. 


‘Certain modifications and add .ions were made to the treat- 
ment presented by von Elbe and Mentser (10). In arriving at 
l;quation [1], the sine of the angle of flame tilt was held to be 
equal to [1 — (R/o)?]'/: rather than 1.4[1 — (R/c)]'/. To 
arrive at Equation [2] the boundary velocity gradient for turbu- 
lent flow was equated to \U Re/16R, with X = 0.316/Re'/s. J is 
the average gas velocity, and Re is the Reynolds number. 

5 Efforts were made to induce flashback by flame tilt by fanning 
the base of the flame. These are all the points in Fig. 2 except 
those for the 5.08 em burner. This variation in experimental pro- 
cedure does not affect the conclusion that turbulent flashback 
gradients and laminar tilted flame gradients can be different. 
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Fig. 1 Comparison of curve for flashback gradients in turbulent 

flow with points for tilted flame gradients in laminar flow: 
Hydrogen-oxygen mixtures 
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Fig. 2 Comparison of laminar flashback curve and curves for 
tilted flame K gradients for turbulent and laminar flow: Natural 
gas-air flames 


Fig. 3 Direct short exposure photographs of a turbulent flame, 
stoichiometric natural gas-air mixture, Reynolds no. 10,000, tube 
diam 7.62 cm 


1 
or, 


Correspondingly, turbulent flashback need not always be 
caused by conditions in the laminar sublayer of wide diameter 
ports. Conditions in the core of the stream may lead to in- 
stability. This possibility may explain the disagreement 
noted in paragraph 2. 

The foregoing tilted flame theory involves the assumption 
that turbulent flashback is produced by events in the laminar 
sublayer. The sublayer is more likely to be significant on 
ports of small rather than large diameter. The natural gas- 
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Fig. 4 Comparison of laminar and turbulent blowoff curves: 


Natural gas-air flames 


Fig. 5 Photograph of water-cooled skirt stabilized flame of 


stoichiometric natural gas-air 


air turbulent flame tests were made with burners with diam- 
eters about 10 to 60 times those employed in the hydrogen- 
oxygen turbulent flame experiments. Equation [4] shows the 
sublayer thickness to be weakly dependent on port diameter 
and inversely proportional to roughly the average gas velocity. 


The Blowoff Problem 


1 Blowoff gradients of unpiloted pipe flow turbulent flamvs 
also show that stabilization of turbulent flames is not alwa)s 
in the laminar boundary sublayer. Ample evidence thit 
stabilization can occur in the laminar sublayer is recorded .n 
(2,3,4 and 6). This is based on the equivalence of blowc ff 
gradients in laminar and turbulent flow. But again thee 
exists contrary evidence in (10) and in some of the data n 
Fig. 4.6 

2 The volumetric flow at blowoff of pipe flow turbule:it 
unpiloted flames was raised considerably by adding a coll, 
short, slightly but sharply expanded cylindrical section >f 
tubing (a skirt stabilizer) atop a long pipe burner (12). Tie 
base of the flame was within the skirt stabilizer, general y 
about halfway up, the flame extending over virtually the e 
tire inside diameter of the stabilizer. Flows were stabiliz 
that were over 18 times the blowoff flow for a 300-cm log 
cylindrical pipe of the same diameter as the 15-cm long skit 
stabilizer (7.6 em ID). These skirt stabilizers were operated 
hot, water-cooled, and with inserts of varying geometry. |n 
all instances, large gains in stability were obtained. 

Fig. 5 is a photograph of the external portion of a 180,000 
Reynolds number flame of a stoichiometric natural gas-air 
flame, stabilized on a 15-cm long water-cooled skirt of 7.6 
em ID atop a 300-cm length of 5.08-cm tubing. The flame 
base shows an interesting abrupt expansion at the exit of the 
skirt, which suggests that the cross-sectional area of the flame 
within the skirt is reduced by a boundary layer of burned gas 
along the wall. This inward thrust on the flame due to the 
confinement of burned gas by the wall of the skirt may be 
highly important in achieving the great gain in flame sta- 
bilization that this geometry provides. 
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The Influence of the Launching 
Conditions on the Orbital Characteristics 


GEZA S. GEDEON! and RICHARD E. DAWLEY? 


Chance Vought Aircraft, Dallas, Texas 


Che relationship between the apogee and perigee alti- 
tudes and the burnout conditions are obtained from the 
equations of motion of a satellite orbiting in a vacuum 
around a spherical nonrotating Earth. This relationship 
is presented on a composite chart. It can be seen from the 
chart that for high reliability in achieving an orbit, high 


excess velocities are favorable. 


Nomenclature 
V = velocity 
R = distance measured from the center of the Earth and a 
point in space 
Rk, = radius of the Earth (here used as 3960 miles) = 
G = acceleration of gravity at sea level, 78,938 miles/hr? 
8 = angle between radius and reference plane 
@ = angle measured from the normal to the line connecting the 


point in space and the center of the 
local horizontal ) 


Earth (deviation from 


k = dR/dT 
R = d@R/dT* 

= d9/dT 

Subscripts 

@ = in the direction of 6 

B = burnout conditions 

A = apogee 

P = perigee 

S = conditions for circular orbit — 


_ equation of motion in polar coordinates is 


= R(6)? — {1] 
where 
R(6)? = centrifugal acceleration 
Ro?Go 


R? gravitational acceleration 


R = resultant radial acceleration 


These components of the accelerations are shown on Fig. 1. 
— 
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Ra.p = 


Vv 
qn 


Fig. 1 Polar coordinate system 


Since there are no tangential forces acting on the vehicle, 
the angular momentum is constant, i.e. 


where K isa constant. Hence 
,  K? 
The resultant radial acceleration R may be expressed as 
dR _ 1 d(k)? 
[4] 


dR 2 dR 


Using the above expression for R and 6 in Esuation [1], the 
equation of motion becomes 


d(R)? kK? Ry?G, [5] 
Integration of Equation [5] vields 
(R)? 1 
2 2 R? R 6) 


The angular momentum K and the constant of integration C 
‘an be evaluated at the burnout conditions. Thus 


K = R® = = cos [7] 


and 


. [8 
2 Reg 
Using the above values for K and C, Equation [6] becomes 
1 cos? op Ve? Ro*Go [9] 
R? R 2 Re 


The length of the radius vector to the apogee and perigee may 


be obtained from [9] since R is zero at these two orbital — 
points. Therefore 


cos? bg 


—— (10) 


2G.Re? 


Equation [10] may be simplified by introducing the circular 
satellite velocity at the burnout altitude 


R 

Vs 
¥ 
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to the apogee is obtained by 
using the positive sign on the 
radical, and the length of the 
radius vector to the perigee is 
obtained by using the nega- 
tive sign on the radical. 

Solutions of Equation [12] 
are illustrated on Fig. 2. To 
obtain the apogee and perigee 
altitudes from the initial condi- 
tions easily, Fig. 2 has two 
nomograms: One for convert- 
ing the burnout velocity and 
the burnout altitude into the 
velocity ratio (Ve/Vs), and 
one for converting the radius 
ratio (R4.p/Rz), and the burn- 
out altitude into the apogee 
and perigee altitudes. An ex- 
ample illustrating the use of 
the figure is given on Fig. 2. 

It can be seen from Fig. 2 
that high excess launching 
speed increases the probability 
of obtaining an orbit. It is ap- 
parent that at high excess 
launching speeds the angular 
errors affect, mainly, the apo- 
gee altitude, whereas at launch- 
ing speeds less than those re- 
quired for a circular orbit, the 
angular errors affect mainly 
the perigee altitude. 


Introduction 


APOGEE AND PERIGEE ALTITUDE — MILES 107 


MILTON FARBER! 


Thermodynamics of Al,O, 


MILES 


BURNOUT ALTITUDE — 


Hughes Tool Co., Aircraft Division, Culver City, Calif. 


HE utilization of aluminum as a propellant constituent 
is proceeding at an accelerated pace. 
fore arises for performance calculations of propellant systems 
containing elemental aluminum and various components of 


aluminum. 


The need there- 


Early performance calculations considered the 


formation of either solid, liquid or gaseous aluminum oxide 
as the only combustion product of aluminum or its com- 


pounds. 


However, since the reactions of aluminum with the 


other combustion species can be quite involved, it is the pur- 
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BURNOUT VELOCITY — FPS 107 
Fig. 2 Influence of launching conditions on the orbital characteristics 

per 


pose of this paper to discuss some of the significant reactions 


and to point out some of the thermodynamic relationships 


that are involved. 


Theoretical Considerations 


In most performance calculations of propellants containing 
- aluminum or its compounds assumptions have been made that 
Al,O; remains undissociated in the liquid stage at tempera- 


tures above 2500 K. Brewer and Searcy 
the AJ-Al.O3 system in an effusion cell. 
pressure measurements, they obtained a 
3770 + 200 K for Al,Os. 


(1)? investigated 
From the vapor 
boiling point of 


Their evidence for the dissociation 


of Al,O; to Al,O and AlO at high temperatures has been sub- 


stantiated by Hoch (2) and Cochran (3). 


From a combina- 


tion of their work and values listed by Kelly (4) and NBS 
(5) and Brewer’s Tables, some thermodynamic data have 
been deduced for the heats of formation of the dissociation 


products of Al,O3. 


A list of the heats of formation, entropies 


and specific heats (C,) at standard conditions of these proc- 


JET PROPULSION 


2? Numbers in parentheses indicate References at end of pape’. 


: 
16 
| 
18 12 1.10 | 
22 22 = 
16 
2 — 
= 
= 
- 
ip 
2 


hips 


ning 
that 
yera- 
ated 
apor 
t of 
ition 
sub- 
NBS 
have 
ition 
ypies 


ape”. 


SION 


ucts are presented in Table I. A value of 27 cal/deg/mole 


TableI Thermodynamic properties of aluminum and its 


oxides 

Compound 4H,° 298 So** C, 298 

(kcal/mole) (cal/deg/mole) (cal/deg/mole) 
Al(s) 0 6.77 5.82 
Al(l) 1.90 d 10.3 7.0 
Al(g) 39.315. 
AloO;(s) —399.1 18.88 
Al,03(1) 17.3 27 (estimate) 
AlO(g) 9.5 52.14 7.4 (estimate) 
AbO(g) —34 65 9.6 (estimate) 


* Measurements in an effusion cell by the author sub- 
stantiate this value. 

** Entropy and specific heat values are from Kelly 
Bureau Mines Bulletin 477) and NBS tables (circular 500 
NBS). Those for AlO(g) and Al.O0;(1) have been esti- 
mated as described in the text. 


for the specific heat of Al,O3 (1) at 298 K has been estimated 
from an examination of the literature for similar compounds 
and interpolation from a rule for estimating specific heats of 
lijuids given by Wenner (6). This rule gives a value of 34 
cal per mole for the C, of Al.Os; (1) while Kopps’ rule gives a 
rulue of 24.4 cal/deg/mole for C, of AlsO3(s). Since the ex- 
perimental value of the C, at 298 K of Al,O3(s) is 18.9 cal/ 
deg/mole the value obtained by Wenner’s rule is undoubtedly 
too high and should be corrected to a lower value. Herzberg 
(7) gives a value of 978 em~ for we of AIO from which a 
value of 7.4 cal/deg/mole is calculated for C, of AlO at 
298 K. The C, of Al,O was estimated by adding the vibra- 
tional contribution of four degrees of freedom (assuming that 
the molecule is bent) to that of rotation and translation. 
Employing the AlO frequency as a first approximation for 
the Al,O frequencies a value of 0.43 cal/deg/mole is obtained 
for the vibrational specific heat for each degree of freedom. 
The addition of 1.6 cal/deg/mole for the vibrational C, 
yields a total specific heat of 9.6 cal/deg/mole for Al,O(g) 
at 298 K. 

Besides the straight decomposition of Al,O3; to elemental 
aluminum and its suboxides of AlO and AIO, the reactions 
of Al,O; and the other combustion species must be considered. 
These include possible reactions with CO, CO., H,O and 
solid C. Example calculations of the free energies of these 
reactions are presented for the combustion temperature range 
of 3000 K and in the vicinity of the boiling point, 3770 K, 
as an indication of the possible decomposition of Al.Os. 


The Reactions of Al.O; With Propellant Combustion Gases 
The stability of Al,O; at 3000 K is tested in reactions with 
H., CO and solid C. 
1 The reaction of H, and Al,O3 may proceed as 


2H. + AlLO,(1) = 2H,0(g) + AlLO(g) 
AHo = +223.4 keal = 75.3 cal/K =—15.4 cal/K 
AF = 227,990 + 15.4 7In 7 — 178.47 


AF wo) = 62.8 keal AFs377 = 33.6 keal 


K3000 = 2.7 X 10-5 Kozo = 1.1 X 107? 


The partial pressure of Al,O becomes increasingly important 
above 3000 K. 

2 The reaction of solid C and He, with Al.O3. 

In some low temperature reactions or those in which large 
quantities of Al are added there is the possibility of the exist- 
ence of solid carbon. Therefore its reaction with Al,O; 
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will be investigated in the following equilibrium at 2500 and 


3000 K 
C + H, + Al,O,(1) = CO + + AlLO(g) | 
ASogg = 107.4 cal/K AC, = —11.6 cal/K 


= 254.8 keal 
AF = 258,260 + 11.6 7 In T — 185.1 7 
AF 37 = —79.4 keal 


AF 2500 = 22.4 keal AF 3000 = —18.4 keal 


P 1,0 
Pn, 


K 


From the above calculations it can be seen that the decom- 
position of Al,O;(1) with C and Hz becomes considerable at 
temperatures above 2500 K. Although this temperature is 
below the average chamber temperature for hot propellant 
systems, it is probably an upper limit for systems containing 
solid carbon, since the source of carbon would be a fuel that 
would be stripped of its hydrogen atoms, resulting in low 
chamber temperatures. 
3 The reaction of CO and Al.O; can be considered as 


2CO + Al.O,(1) = 2 CO, + Al,O(g) 
AA og = 203.8 keal ASoog = 55.3 eal /K AC» => 13.8 cal/K 
AF = 207,910 + 13.8 T In T — 147.7 T 


K P' coe 


= 1.01 K 107 Kaz = 2.3 X 1075 


This calculation also gives a negligible vapor pressure of 
Al.O; at 3000 K and at the boiling point. 

Dissociation of Al.O; 


Brewer and Searcy (1) indicate from their measurements | 
that Al,O; tends to dissociate to the suboxides at tempera- 
tures below its boiling point of 3800 K. A consideration of 
the extent of this dissociation is shown by the following 


reaction 
Al.O,(1) = 2 AlO(g) + O(g) 
= 456 keal 


AF = 443,000 + 12 7 In T — 212 T 7 
AF 3000 = 95 keal AF37 = 16.4 keal ee 
K = P%10 P “< 
= “Alt 
AlO > 
Kyo = 1.14 X 1077 Koz = 1.12 X 107 


This calculation shows that the pressure of AlO would be in — 
the vicinity of 10~4 to 10-3 atm at 3000 K depending on the — 
partial pressure of O in the gaseous species. However at_ 
temperatures near the boiling point the decomposition of 
Al.Os is considerable. 


Conclusions 


The conclusions derived from the experimental measure- 
ments and the above calculations show that Al,O; may be 
treated as an undissociated liquid resulting from the combus- 
tion of Al or Al compounds for the purpose of propellant per- 
formance calculations below 3000 K. At temperatures above 
3000 K it will be necessary to calculate equilibrium con- 
stants for the various reactions and dissociations in order to 
determine the extent of the Al,O; decomposition. = 
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Ge Characteristic Velocity as a 
Wrers Measure of Available Work 


S. L. BRAGG! and H. RATCLIFFE? 
Rolls-Royce Ltd., Derby, England 


HE purpose of this note is to demonstrate that the power 
available from a gas turbine is dependent primarily on the 
characteristic velocity of the entering gas, and is otherwise 
insensitive to gas properties, such as the specific heat ratio. 
The power available from an ideal turbine passing a flow 1/ 
of gas, at an initial total temperature 7;, and expanding the 
gas over a pressure ratio r (>1) is given by 


—(r) 7 [1] 


where C, is the specific heat of the gas at constant pressure, 
ecneeaed | in work units, and y is the ratio of the specific heats 
of the gas. 

In many cases, particularly where the turbine is driven by 
the products from a fuel rich gas generator, the exact composi- 
tion and properties of the exhaust gases are unknown. Some 
clue to these is, however, given by the calculation, from the re- 
sults of gas generator or turbine tests, of the characteristic 
velocity of the gases. This is defined by the relation 


where . 4 is the rare area choked by the flow J at total pres- 
sure P. 

From the normal] relations for isentropic choking flow of a 
perfect gas 


VRT 

APV/mg y+1 f(y) [3] 


Power function (550 HP g M c**) against pressure ratio 7 
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where R is the universal gas constant in work units, and m is 
the molecular weight of the gas. 

Also, for a perfect gas 
R 


m 7 


Combining these equations, one obtains an expression for the 
power function 
550HP g — 
Me 
By a happy chance the expression on the right hand side ‘s 
nearly independent of y over the whole range normally of in- 
terest. It varies by less than 12 per cent as y is increased 
from 1.0 to 1.2, at any particular value of r between 1 and 25, 
Thus, to a sufficient approximation, the idea] turbine pow:-r 
is given by 


550 HP g 
or 
HP = Mc*? fy 


Values for this function, based on a typical value of y = 
1.13 are given below. 


r rig + 1.0 1.5 2.0 3.0 
0. 162 0.268 
2 0 0.078 0.101 0.094 
oy 
5.0 9.0 17.0 25. a 
far) 0.592 0.782 0.974 
0.020 —0.149 —0.402 —0.620 


q- 
oy 
Since the ideal turbine horsepower can be obtained from 

Equation [6b] and compared with the horsepower actually pro- 

duced on test, it is possible to estimate turbine efficiency to 

within about + 2 per cent without any more detailed knowledge of 
gas properties than is afforded by the direct calculation of the 
characteristic velocity. Alternatively, the work available from 

a unit mass of gas expanding over a given pressure ratio is, to 

the same order of accuracy, directly proportional to the square 

of the characteristic velocity. 


Toxicity and Personal Decontamination 


of Boron Hydride Propellant Fuels 


ROTHBERG,' JOSEPH L. COLBOURN? and 
ROBERT SALVATORE? 


SIDNEY 


Directorate of Medical Research, U. S. Army Chemical 


Warfare Laboratories, Army Chemical Center, Md. 


The toxicities of the liquid boron hydride fuels (BHF) 
of the HEF-3 and HiCal class were determined for rats, 
guinea pigs, rabbits and cats by the intravenous, intra- 
peritoneal, oral, ocular, inhalation and skin routes. The 
high toxicity of these compounds when applied to skin 
prompted a search for an effective decontaminant. A 
number of candidate materials were tested and found to 


be unsatisfactory. Rabbits and cats survived, however, 


; when their skin was swabbed with a 2.8 per cent solution 


! Pharmacologist (Toxicology), Toxicology Division. 
2 Biologist, Toxicology Division. ; 
’ Enlisted Specialist (U.S. Army), Toxicology Division, 
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of ammonia water 15 minutes after application of four 
times the lethal dose of BHF. 


Introduction 
URRENT interest in high energy propellant fuels has 
A focused the attention of toxicologists on their potential 
hazard to health. The toxicities of the boron hydride fuels 
of the HEF-3 and HiCal class and procedures for their 
removal from the skin surface are the subjects of this report. 
sorax and boric acid have long been known to be moder- 
ately toxic. Accidental deaths have occurred when boric acid 
was inadvertently administered to infants with their feeding 
formulas (1)4 and when applied to extensive areas of skin in 
burn therapy (2,3,4). Krackow (5) has reported that deca- 
borane, diborane and pentaborane are highly toxic com- 
pounds presenting a serious handling hazard by inhalation, 
or: | ingestion and skin absorption. Other workers have re- 

ported on the toxicity of the boron hydrides (6,7,8,9). 


Toxicity Studies 

oxicity studies were conducted on animals to find the 
L150. This is the amount of material (in milligrams of fuel 
per kilogram of animal body weight) required to kill 50 per 
cet of the animals within a specified observation period, in 
this case, 24 to 48 hr. Intravenous and intraperitoneal tox- 
icity experiments were carried out to ascertain the inherent 
toxicity of the fuels. 
INrRAVENOUS Toxicity The LD50 of BHF by intravenous 
injection was 13 mg/kg for rats and ranged from 4 to 6 
mg/kg for rabbits. At these dosages, the animals were 
ataxic in 10 min, showed tonie and clonic convulsions in 15 

) 20 min and prostration and death in 25 to 30 min. 

INTRAPERITONEAL Toxiciry The LD50 of BHF by intra- 
peritoneal injection ranged from 20 to 71 mg/kg for rats and 
18 to 40 mg/kg for guinea pigs. In these animals, clonic 
convulsions were seen in 30 to 40 min and prostration fol- 
lowed by death in 45 to 120 min. 
OraLt Toxiciry Oral toxicity was not directly measured, 
but an estimate of the toxicity of BHF by ingestion was ob- 
tained by administering the material into the stomach of rats 
by means of a flexible plastic tube. The LD50 was 40 ml/kg. 
Time of onset of toxic signs was somewhat longer by this 
route, clonic convulsions appearing in 3 to 4 hr and death in 
3 to 5 hr. 
INHALATION Toxiciry® The vapors of purified BHF are 
highly toxic, the LD50 being approximately 23 ppm for 
rats and 6 ppm for mice for a 4 hr exposure period. 
Eye Damage A very small amount (0.02 ml) of BHF was 
instilled into the eyes of normal rabbits. Five minutes later, 
the eyes were washed with physiological saline solution and 
examined daily for several weeks. All animals showed 
severe eye damage as characterized by iritis, conjunctival in- 
flammation and edema, ulceration and corneal opacity. The 
opacity, which covered more than 75 per cent of the cornea, 
lasted longer than 16 days. 
Skin Toxiciry® BHF was applied to the skin of rats, 
guinea pigs, rabbits and cats. Local signs of poisoning ap- 
peared first. The skin immediately turned greyish-white 
at the site of application, where a vigorous reaction seemed to 
be taking place. Clonic convulsions appeared 3 to 6 hr 
later and death in 10 to 30 hr. LD50’s were 317 to 502 
mg/kg for rats, 158 to 251 mg/kg for guinea pigs, 57 to 105 
mg/kg for rabbits and 126 mg/kg for cats. 


Protection Offered by Rubber Gloves or Ointments 


Experiments were performed to evaluate the protection af- 


‘ Numbers in parentheses indicate References at end of paper. 
’ The inhalation toxicity expe riments were carried out by Leo 
Fcinsilver, Gassing Branch, Chemical Warfare Laboratories. 

’The high degree of skin toxicity might be better visualized if 
one were to convert mg to fl oz and extrapolate to humans. 
Thus, if BHF is as toxic to man as to rabbits, the LID50 for a 
150 Ib man would amount to + OZ. 
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forded by surgical rubber gloves and protective skin creams. 
No evaluation was made of the protection afforded by heavier 
rubber or synthetic rubber gloves. 

An application of 4 « LD50 (420 mg/kg) of BHF was 
made to the bare skin and to skin “protected” by a single 
layer of latex rubber from surgical gloves or by one of several 
silicone-containing ointments. In all cases, healthy clipped 
rabbits with intact skins were used. The unprotected animals 
died within 2 hr, the ointment-covered ones died in 3 to 48 
hr and the rubber-covered ones in 72 to 96 hr, indicating 
that the toxic material eventually penetrated the ointment 

the rubber. BHF 


caused the rubber to deteriorate to 


some extent, allowing penetration to occur during an over-— 


night period of contact. All animals showed local signs o 
skin damage, convulsions and prostration prior to death. 


control animals survived. 


Skin Decontamination 


Exploratory experiments were carried out on rabbits to — 
Applica-_ 
tion of 4 X LD50 (420 mg/kg) of BHF was made on the back — 
carried 


screen a number of candidate skin decontaminants. 


of each animal. Decontamination procedure’ was 
out 30 sec later. Results, shown in Table 1, indicate a lack 
of effectiveness of most of the materials tested. 


to death. 
Table 1 Effect of candidate decontaminants after ap- 


plication of 4 x LD50 (420 mg/kg) of BHF to the skin of 
rabbits, exposure time: 30 sec 


Approximate time | 
Candidate decontaminant of death(hrs) 

none (controls) 13-23 
water 13- 
acetone-alcohol (50-50) 13-2 
-ationic detergent (1 per cent triton 

X-400) 13-2 
anionic detergent (1 per cent triton 300) 3-45 
neutral detergent (1 per cent Tween) 2-5 
gasoline 2-27 
28 per cent tincture green soap—fol- 

lowed by water 48-116 
sodium bicarbonate—5 per cent 2-33 
glycerine followed by sodium _ bicar- 

bonate 11-273 
hard soap (castile)—followed by water 23-3 


Ammonia water was also tested as a decontaminant, based 


All. 


The animals 
in this experiment generally showed signs of poisoning prior — 


upon the fact that NH; reacts with decaborane to form hexa- _ 


ammoniates and diborane to form diammoniates. Results 
were encouraging; 29 out of 30 animals survived (Table 2). 
Even a weak solution (2.8 per cent) allowed 100 per cent re- 
coveries, and none of the surviving animals showed local or 
systemic signs of poisoning. A yellow color appears when 
ammonia is applied to a contaminated area. This, it is pre- 
sumed, is evidence of a chemical reaction. 

Additional experiments were carried out to estimate the 
longest period of exposure which could be successfully decon- 
taminated (Table 3). A 2.8 per cent solution (1:10 dilution 
of stock ammonium hydroxide) was selected as an effective, 
yet relatively nonirritating decontaminant. The results indi- 
cated that 4 * LD50 (420 mg/kg) can be decontaminated 
by a 2.8 per cent solution of ammonia water even after a 15 
minute contact of BHF with the skin. Longer exposures prior 
to application of ammonia resulted in decreasing numbers of 
recoveries. All animals decontaminated with water after a 
30 sec exposure died within 2 hr. Doubling the dose of 
BHF to 8 & LD50 (840 mg/kg) decreased the number of re- 


7 The proc edure consisted in flushing with approximately 1300 
ml of decontaminant while swabbing vigorously with cotton. 
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Pox 
Agent animals Decontaminant 
BHF 10 14 per cent solution ammonia water 
BHF 10 7 per cent solution ammonia water 
BHF 10 2.8 per cent solution ammonia water 
BHF 2 none 
none 6 none 


Table 2 Decontaminating effect of ammonia water after application of 4 x LD50 (420 mg/kg) of BHF to the skin of rabbits, 
exposure time: 30 sec 


No. of 

Approximate time of death (hr) recoyv- 

24 48 72 96 120 144 168 eries 
0 0 0 0 1 0 0 9 
0 0 0 0 0 0 0 10 
0 0 0 0 0 0 0 10 
0 0 0 0 0 0 0 6 


the number substantially. 


coveries somewhat, and 16 X LD50 (1680 mg/kg) reduced 


7 Table 3 Decontaminating effect of 2.8 per cent ammonia 
water after application of 4 x LD50 (420 mg/kg) of BHF 
to the skin of rabbits 
Exposure No. of 
No. of time recov- 
Agent animals (min) Decontaminant eries 
BHF 6 0.5 2.8 per cent 6 
solution 
NH,OH 
BHF 6 1.0 6 
BHF 6 6 
: BHF 6 2 6 
BHF 6 6 
BHF 6 5 6 
BHF 6 10 6 
BHF 6 6 
BHF 6 | 4 
BHF 6 45 2 
- BHF 6 0 5 water 0 
: none 6 2.8 per cent 6 
solution 
NH,OH 


| 


Studies with cats indicate that they can also be effectively 
decontaminated with ammonia. 


Studies of the Decontaminating Action of Ammonia 


_ BHF was mixed with an excess of ammonia water in test 
tubes and allowed to react for about 30 min. These mixtures 
were injected into the peritoneal cavity of guinea pigs or 
_ applied to the skin of rabbits. The resulting LD50’s showed 
no decrease in lethality by these routes. It appears then 
that the addition of ammonia to the fuel may produce a reac- 
‘ tion product, which although toxic and capable of penetrating 
skin, can be easily washed away. 


Discussion 


Boron hydrides of the HEF-3 and HiCal class are com- 
pounds which have a relatively high degree of toxicity by all 
the practical routes of entry into the body. Hence, it is 
necessary to take suitable measures to prevent accidental 
exposures, as for example, protective clothing. The area of 
operation should be well ventilated and protective masks 

_ readily available in the event of a spill. If, however, the 
skin does get splashed, immediate action is required. A 
simple application of ammonia water to the contaminated 
- area is not sufficient, since the reaction product is toxic when 
_ left on skin. Copious flushing, along with swabbing is essen- 
; tial to remove the toxic material. Goggles should be worn 
to prevent splashing of the eyes with the ammonia water. 


First Aid Decontamination Procedures 
Immediate action is required in treating victims of accidental 


skin contamination with BHF. 
suggested: 

1 Assure that both the victim and the person ministering 
to him are wearing goggles. 

2 Flush contaminated skin with large amounts of 2.8 },er 
cent ammonia water solution (preferably from a spray attach- 
ment) while swabbing vigorously with a ball of cotton o: a 
sponge attached to a long holder. 

3 Continue washing affected and adjacent areas until 
all of the BHF has been removed (at least 5 min). 

4 Rinse thoroughly with a plain water shower to remove 
ammonia solution. 

5 Follow up with soap and water, if available. 

6 Report to plant physician for examination. 


The following procedure is 


Conclusions 


On the basis of the work described herein, it is concluded 
that: 

1 Boron hydride fuels of the HEF-3 and HiCal class are 
highly toxic materials which present a serious health hazard 
by inhalation, oral ingestion, ocular contact and skin ab- 
sorption. 

2 A dilute solution (2.8 per cent) of ammonia water as a 
flushing agent is effective in removing these fuels from the 
skin of animals used in experiments; hence it is presumed to 
be effective for man. 

3 Surgical glove material and the silicone-containing 
ointments used in this study provide only a limited degree of 
protection against BHF, the former being more effective than 
the ointments. 
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Technical { 
of 
v- 
. Possibility of Aerodynamic Descent on material of the sail must be suitable for prolonged exposure to 
these conditions. 
’ _ the Moon - 3 On occasions when the shadow side is emitting radi- 
“s ee ee ae > ation, there is a force exerted on the sail opposite in direction 
LANCE W. SMALL! 
to that produced by solarradiation, 
ECENT investigations of Elsmore? by use of radio Reference 
—— astronomy have noticed an effect at certain wave lengths 1 Garwin, Richard L., “Solar Sailing—A Practical Method 
whch implies refraction at the lunar disk. This refraction of Propulsion Within the Solar System,” Jer Proputsion, vol. 
may be explained by assuming an atmospheric pressure of 28, March 1958, p. 188. 
lure is 10 '* atm on the moon. In this connection, Borst and _ 
Edwards* have suggested that such an atmosphere may be 7 
tering cor posed of the rare gases krypton and xenon. 
f it were possible to glide into the moon without the ex- A Remark = ee Iteratiy . Method of 
8 ner cessive use of rocket braking, the vehicle used would either Determining Equilibrium Compositions 
ttach- be smaller than those now planned or could utilize the power oo of Reacting Gases | 
nora led for braking for the trip home. : 
low does the lunar atmosphere compare with ours? Ob- J. L. SACKMAN! 
until servations of Sputnik I have revealed that the ARDC atmos- ~ } re 
phere is too small by a factor of 40. This means that a Columbia University, New York, N. Y. 
emove satellite with a perigee altitude of 400 km initially will have 
=) its lifetime reduced by a factor of 40. Note that this applies N A recent article in this Journat (1),? an iteration pro-- 
to 400 km. The atmospheric pressure at 400 km would be cedure for obtaining the roots of a system of simultaneous 
1.44 & 10-8 & 40, that is, 5.8 X 10-9 millibars. On the moon equations governing the equilibrium compositions of reacting 
the pressure would amount to 10~° millibars. Therefore, the gases was presented. It should be noted that the numerical — 
lunar atmosphere provides a greater drag by a factor of 5. procedure there presented is in reality an application of a 
ded \lso, if the mass of the lunar rocket is 10‘ greater than that method which has been well known in the field of numerical | 


ss are where Cd is the drag coefficient) will be greater by a factor o 


nazard 10°. Thus, we have the possibility of aerodynamic descent tion” (2). This method, as presented in (2), is briefly re-_ 
is. ‘obs on the moon. Of course, a winged vehicle would be better viewed here for the special case of two simultaneous eek 
able to accomplish this. The roots of the simultaneous equations F(x, y) = 0 and 
Tr as a - G(x, y) 0 are desired. Starting with an approximation — 
m the Received March 24, 1958. y™), the next approximation (2+), is deter-_ 
dt Student, Roosevelt High School, Yonkers. Student Mem- mined in the following manner 
nec ber ARS. 
Elsmore, B. and Whitfield, G., Nature, vol. 176, 1955, p. 457. zint) = 4 = yl) 4 
aining Borst, L. and Edwards, W., Science, vol. 127, 1958, p. 325. a 
vree of where 6x2 and éy are obtained from the two simultaneous _ 
linear algebraic equations 
Solar Sailing—Practical Problems F, yi) + (2, y™) = — [2] 
S. W. GREENWOOD! y™) + SyMG, (2, y) = —G(2™, y™) [3] 
. Acid College of Aeronautics, Cranfield, England The comma symbol denotes partial differentiation with re- 
n, vol. spect to the indicated variable. 
, N A recent Technical Note, Garwin (1)? considers the solar ne now the system of Equations [1 to 6] of (1). 
radiation forces acting on sails in space, and accelerations 
theoretically possible with sails of g given area and vehicles of [4] 
oat a given mass. As are sult of his analysis he cone ludes that solar i=l 
sailing is a practical method of propulsion within the solar = + — + (1 — + = Be 
‘ec Acid | SYstem In his examples he selects a thin sail made of plastic (5) 
1 Labs § 2nd aluminized on one side. It seems to me that before the 2 F 
method can be assessed as practical it is necessary to consider polpo,)~"* = Ko.............. [6] 
Boron § the following problems: 
decupa- | If all the incident solar radiation is reflected to produce 
maximum thrust, the temperature of the sail would on many = Ku,o 
[ydride § occasions fall to near absolute zero owing to the radiation of 
lustrial energy from the shadow side. The problems arising here are = Ku,o (9) 


‘irbely, J tures, and effects produced by temperature cycling as the sail 
oral’, is furled and unfurled. 

Health, 2 The sail is exposed to cosmic radiation, meteorite bom- 
me bardment and the full spectrum of solar radiation. The 
Received May 5, 1958. 

oar iak ‘ Lecturer, Department of Aircraft Propulsion. Member ARS. 
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of Sputnik I, the ballistic drag coefficient (mass/Cd xX area 


the suitability of the material for operation at low tempera- 


computations for a considerable period of time. It is com- 
monly referred to as “the Newton method of approxima-— 


f 


It is clearly seen that by using [6 to 9] to express four of the 
partial pressures in terms of the remaining two, this system | 
of equations can immediately be reduced to one involving — 

two simultaneous equations in two unknowns. To parallel — 


Received April 30, 1958. 
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the choice made in (1), pz,o and po, are selected as the pri- 
mary unknowns, and then [6 to 9] yield 


' i is now desired to locate the point of intersection of the two 
plane curves F 
method, as prev discussed, is employed. 
_ Equations [2 and 3] it is required to obtain F,po,, F,px,o, 


= 0 and G = 0, and to this end Newton’s 
To utilize 


G,po,, G,pu.o. These quantities are readily obtained from 
[14 and 15] and [10 to 13] by simple partial differentiation. 
For example 


Opu:0 , por 
Op 420 

1 
=14+0+4P% 4 — Pow 
Puo 2 
In like manner 
F Po: = B G,p =C G,Por =D...... [17] 


It is also noted that 


where A, B, C, D, E are quantities defined in (1). (The ex- 
pressions for C, D, E given in (1) appear to be in error. As- 
suming Equation [2] of (1) is correct, the multiplier in paren- 
theses before pon should be (1 — 8) instead of (1 — 28) as is 
shown in (1) for C, D, E.) 

These results are utilized in [1 to 3] to yield 


Pio = Pino + Por = Por + Spor -- - - [19] 
6 
A + Spor = p — [20] 


It is seen that [20 and 21] are exactly the Equations [22 and 
23] of (1), thereby demonstrating that the method proposed 
in (1) is in actuality a specific application of Newton’s method 
to the particular system of equations considered in (1). A 
discussion concerning convergence of Newton’s method may 
be found in (2). 
References 
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On a Generalized Optimization 
Procedure for N-Staged Missiles 


Israel Institute of Technology, Haifa, Israel 
PH = Ku(Ku,0)'? (Po) ot [12] 
pow = (pu (po,)"4...... EFERENCE (1)? established the general relation betwe: n 
weights of consecutive stages of a missile in drag-frie 
Inv view of [10 to 13], Equations [4 and 5] may be written as space so that the launching weight be a minimum. It mig!it 
6 be worth noting that the solution is implicitly contained in 
F(pH.0, Po.) = Pi | [14] 
Po.) = (1 28)pu.0 + 2po0, — 2Bpu, + (l — B)pou + po — Bpu [15] 


(2) where the more difficult problem of acceleration-limited 
stages is treated. 

Using the notation of (1), the minimal solution must satis. y 
Equation [15] of (2), which with no restrictions on the perm ;- 
sible acceleration takes the form 


Mz, 
(7) 
Mz, 
Mo, 


N 
> Mins [N+1] 


n=1 Mon 
1-—k,+K, 
+ Mons 


= 


By solving for the Lagrangian multiplier 7 one obtains 


Ka (1 — 


Ky-1) ry Cn-2) + K,_(1 = Mz, 1 


1 Weisbord, Leon, ‘A Generalized Optimization Procedure for 
N-Staged Missiles,” Jer Proputston, vol. 28, 1958, pp. 164-167. 

2 Schurmann, Ernest E. H., “Optimum Staging Technique 
for Multistaged Rocket Vehicles,’ JET vol. 27 
1957, pp. 863-865. 


which is identical with Equation [13] of (1). 
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Comment on “A Practical Approach to 


Grain Design” 
8 
J. A. VANDENKERCKHOVE! 
U nidennael de Bruxelles, Brussels, Belgium 


N HIS interesting paper, Stone (1)? has proposed a method 
for the design of internal burning stars. This method is 


Received June 24, 1958. 
1 Assistant. Consultant of the Poudreries Reunies de Belgique. 
? Numbers in parentheses indicate References at end of paper. 


JET PROPULSION 


» 
x 
a 
= 
i 
/ 
or 
Mo, 
M 
Mz, 
- 
~ 
> 
1 
| 
] 
Vey. 
] 
4 
- 


ween 
nig! it 
ed in 


nited 


itis. v 
rmis- 


[2] 


M 


re for 
-167. 
Lique 


neer- 


aper. 


to 


thod 


od is 


ique. 
aper. 


Fig. 1 Definitions of the design parameters of an internal 


burning star 


limited to well-defined values of certain parameters (the fillet 
radius for instance) and does not take advantage of a few 
basie properties: 

1 The fillet radius must not necessarily be imposed. In- 
deed, in view of the evolving nature of the combustion, the 
star with a fillet is identical to the star without a fillet which 
would have burnt on a distance equal to the fillet radius. 

2 The star points are often cut at the outer diameters to 
increase the web or the port area and, overall, to reduce 
the sliver fraction. The amount cut can be expressed by the 
parameter € defined on Fig. 1 (Stone is only concerned with 
the case ¢ = 1.) This solution is of particular interest when 
tapered port areas are used. 

3 The problem is complicated by taking the character- 
istic length equal to the grain outer diameter. 

Indeed, for the same internal perforation, the web can be 
increased without changing the beginning of the thrust-time 
curve. By doing so in a region where the star points have 
disappeared, one generally introduces an increasing pro- 
gressivity which fortunately can be easily neutralized by 
tapering the grain at the outer diameters. 

Therefore a whole family of curves given by Stone can be 
replaced by one single curve independent of the fillet radius 
and of the web thickness, if the maximum radius of the basic 
performation is taken as the characteristic length. Such a 
curve readily gives the variation of the burning area during 
combustion. In view of these remarks it seems possible to 
propose a much more general solution which is summarized 
in Figs. 1 to 4. 
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Fig. 3 Dimensionless port area during the first phase as a — 
function of the reduced burnt distance. 
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Fig. 2 Dimensionless burning perimeter of a seven-pointed 
star as a function of the reduced burnt distance 
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Table 1 Internal burning star 


n f 
sin 


2 


Fig. 4. Dimensionless sliver area of a seven-pointed star as a 
function of the reduced burnt distance 
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_ Fig. 1 represents a sector of an internal burning star, with 
the symbols defined. Table I gives the resulting equations 
w hich can be easily solved with a desk computer. The port 
Breer equation which is of particular interest for the grain 
designer is given rather than the loading fraction equation 
which isa byproduct. Figs. 2,3 and 4 represent vs. the burnt 
distance y sg the fillet radius f, the burning perimeter S, 
the port area A, during the first phase (during which the 
star points Press not disappeared) and the final sliver area 
Ay during the second phase. These curves are given in a 
dimensionless form for different values of ¢, in the case of 
= seven-pointed star (n = 7). Fig. 3 has been calculated 
for neutral burning during the first phase. Fig. 2 also repre- 
sents in the case € = 1 a regressive and a progressive solu- 
tion (in dotted lines) obtained by taking 6 ¥ 6 during the first 
phase. The burning perimeter S varies linearly vs. y during 
the first phase and tangentially joins the curve describing 
the second phase (this curve is fully drawn in the case € = 
1). Finally, the curve A represents the limit between the first 
and second phases for neutral burning during the first phase. 
We believe that the method summarized above gives a 
simple, clear and useful picture which provides a maximum 
amount of information with the minimum number of graphs. 
The method can be readily extended in the case of the so- 
called wagon wheel design. A complete set of charts for 
ranging from 3 to 12 has been prepared under an ARDC 
contract (2). 
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A Comment on Combustion Instability 
in Solid Propellant Rocket Motors 


L. GREEN Jr.! 
Aerojet-General Corp., Azusa, Calif. 
2. 
= 


and 


W. NACHBAR? = 


Lockheed Aircraft Corp., Sunnyvale, Calif. | 


N THEIR recent note describing experiments with a closed- 

end, tubular research burner, Price and Sofferis (1)? re- 
ported a “contention” of one of the writers to the effect that: 
“Changes in burning rate during unstable performance are 
favored by conditions of high quasi-steady stream velocity 
in the conduit of the burner.” This report requires correction. 
On one hand, experience has been cited (2) that (other factors, 
such as mean chamber pressure, being constant) gross reac- 
tion irregularities (i.e., secondary peaks in the low response 
pressure vs. time history manifesting an increased propellant 
burning rate) generally tended to increase in severity and 
frequency of occurrence as the gas velocity leaving the aft 

1 of side-burning charges was increased; i.e., as the throat 
to port area ratio was increased. It was surmised (2) that 
this effect might be attributed to the excitation and/or partial 
sustenance of combustion oscillations by fluid dynamically 
excited sound resulting from high velocity flow across a 
channel discontinuity, a possibility also noted by Price and 
Sofferis. 

On the other hand, analytical results (2,3) indicate that a 
change in mean burning rate is favored by a high magnitude 


of both the mean value and the amplitude of the periodically 


Applied to Metal-Based Propellants 


fluctuating component of a parameter group involving the 
ratio piw/l (where p; is the chamber gas density, v; is the 
parallel stream velocity and / is a length parameter identified 
with the distance from the stream stagnation point to the 
station in question). This implies that changes in mean 
burning rate are consequently favored by a high mean value 
and by a high amplitude, in-phase fluctuation of gas pressure 
and of stream-wise velocity gradient as distinguished from 
a high stream velocity. These conditions would appear in 
the region of a stagnation point. This result is consistent 
with observations that the areas of locally excessive pro- 
pellant burning rate, as revealed by interrupted burning 
experiments with conventional motors, are found near the 
stagnation point of the flow system. The analytical results 
do not conflict with the experimental observations of Price 
and Sofferis, which was the interpretation given (1), but are 
in agreement with their experiments, which showed that the 
observed instances of pronounced oscillations and changes in 
burning rate corresponded to instances of gas phase oscilla- 
tion in those modes where the velocity nodes (and hence 
anti-nodes of pressure and axial velocity gradient) were 
located at or adjacent to burning propellant surfaces. 

Recent observations by Price and Sofferis and others indi- 
cate that, under some conditions, gas phase oscillations can 
effect a decrease in mean propellant burning rate. This possi- 
bility has been considered in a re-examination of the behavior 
of the simplified analytical model (4). A principal assump- 
tion involved in the model is that the burning rate is an 
Arrhenius function of the surface temperature: r = B exp 
(—E/RT.). Ther vs. T, curve is thus S-shaped, with zero 
slope at 7, = Oand 7, = ©, and with an inflection point at 
T, = E/2R. For T, < E/2R the curvature is concave up- 
ward, for 7, > E/2R it is concave downward. In the numer- 
ical example previously examined (3), the assumed values 
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(taken as characteristic of a high energy, fast burning pro- 
pellant) of B = 100 cm/sec and E = 11,900 cal/mole yielded 
a steady-state value of 7,, = 1150 K. In this case E/RT., = 
5.2, and any fluctuations in 7’, around this value of T,, 
(located in the region of upward concavity of the nonlinear 
Arrhenius curve) would effect an increase in mean burning 
rate, as averaged over a complete cycle. 

However, if a different set of parameters (such as might 
characterize a slow burning but easily ‘“activated’’ pro- 
pellant of high flame temperature) were chosen, for example 
E = 6000 cal/mole and B small enough that the steady- 
state surface temperature is high enough to fall in the region 
of downward concavity (e.g., T.. > 1500 K), then oscillations 
in 7’, about this value would effect a decrease in mean burn- 
ing rate. The writers wish to emphasize this possibility 
here, since it was not noted in the original analysis (3) of the 
simplified model behavior. 
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High Temperature Chemistry as 


JOHN S. GORDON! 


Thiokol Chemical Corp., Denville, N. J. 


“last frontiers’ of rocket propellant chemistry are 
high energy oxidizers and novel fuel structures containing 
light metals. These materials may find application in both 
liquid and solid rocketry. High energy oxidizers will con- 
tain oxygen and fluorine and as few other elements as pos- 
sible; their combustion products are already well-known. 
Such is not the case (in varying degrees) for the combustion 
products of the light metals lithium, beryllium, boron, mag- 


nesium and aluminum. The status of these fuel elements 


may be summarized briefly. 

Lithium: Excellent performance with fluorine-based oxi- 
dizers. Not outstanding with oxygen oxidizers. Density 
unattractive. 

Beryllium: Performance probably excellent. 
temperatures. Toxicity unattractive. 

Boron: Extremely desirable in liquid fuels and solid fuel, 
binders. Complex combustion products due to trivalence. 

Magnesium: Probably mediocre due to high atomic weight 
but no conclusive data yet. 

Aluminum: Very attractive additive when oxidized by 
hydrogen peroxide or ammonium perchlorate, due to high 
heat of combustion to liquid Al,O;. Attractive density. 

Silicon and higher: Poor energetics. Useful as central 
atoms in building novel molecular structures. Attractive 
densities occasionally occur. 

The reduction-to-practice of these propellant families is 
being hindered by poor knowledge of the thermodynamic 
properties, heats of formation and molecular constants of the 
products of combustion of the metals in question. Improved 


High flame 
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BeO, BeF:, BeCle, BeFCl, 
compounds 


B.O., HBOs, HeBeO,, 
FBCI, B,O;, H2B,O;7 


cule, 


H;BO;3, OBF, 


data would be useful on the following gaseous species: 
Li,O, LiOH, LiO, LiCl, LisCle, Lie(OH). 
BeCl and corresponding Mg 


AlxOz, AIO, AlCle, AICl;, AleCl., ALF3;, 
AIFCh, OAICI, OAIF, AIC, AIN, LiAIF;, LiAICl, 

There is no present evidence for the gaseous HOAIO mole- 

but this should be checked further. 


OBCI, 


BF:, BCh, 


Work has been 


done in the NasO and NaOH liquid-vapor systems, and 
this should be extended to Li systems. 
underway on determining unequivocally the structure and 
molecular constants of BsO3; and gases. 
type of work is properly being done 
fundamental research institutions, 
techniques and interpretative thinking required. 

Improved fundamental data along the lines suggested here 
will improve the reliability of the theoretical performan e 
analyses being made by various agencies. 


Studies are also 


Most of this 
at universities and 
because of the advanevd 


New Patents 


Control of reheat in turbojet engines 
(2,844,936). H. S. Fowler and D. A. 
J. Millar, Cyrville, Ontario, Canada, 
assignors to National Research Council. 

Automatic variation of the nozzle area 
opening by an amount theoretically 


necessary to maintain the velocity of 
gases passing through the nozzle at a 
constant value. 

Jet engine control (2,844,937). S. J. 


Mikina, Pittsburgh, Pa., assignor to the 
U.S. Navy. 

Device for operating a fuel supply and a 

potentiometer for temperature control in 
response to the force of pressure from the 
compressor in excess of a predetermined 
pressure. 
Fuel pressuring system for supersonic 
ramjet (2,844,938). J. P. Longwell, 
Scotch Plains, N. J., assignor to the 
U.S. Navy. 

Projectile with outer and inner walls 
forming a tubular duct containing a 
diffuser. A compressible wall is located 
between the inner and outer walls, and a 
fuel chamber is defined by the diffuser 
wall and compressible wall. Vapor pres- 
sure compresses the compressible wall and 
forces fuel into the duct. 

Tube-bundle combustion chamber 
(2,844,939). F. Schultz, Scotia, 
N. Y., assignor to General Electric Co. 

Tubes secured together to form a 

hollow combustion chamber surrounded 
by a rubber coating forming a fluid 
impervious layer for preventing leakage of 
combustion gases between tubes. A 
binding means surrounds the coating for 
keeping it in contact with the tubes and 
for maintaining them in a predetermined 
shape. 
Noise suppressors for jet engines (2,- 
845,775). J. M. Tyler and R. E. Meyer, 
Glastonbury, Conn., assignors to United 
Aircraft Corp. 

Hollow vanes, having large apertures 
on one side and small apertures in another. 
Vanes are rotatably mounted from move- 
ment in a clustered position downstream 
of an inner body within a jet engine tail- 
pipe, outwardly to present the large 
apertures to the gas passage. 

Protective device for jet engines (2,846, 
023). V. Milliman, San Diego, Calif., 
assignor to General Dynamics Corp. 

Inlet for filtering foreign material from 
air entering a jet engine, consisting of a 


Epritor’s Note: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 


chamber with a top wall formed of coarser 
mesh than the other walls. 

Vertical take-off and landing aircraft 
(2,846,164). E. Haberkorn, Dammarie- 
les-Lys, France, assignor to SNECMA 
Co. 

Annular wing with flaps along the 
trailing edge pivotable to a retracted 
or extended position from both the outer 
and inner wing surfaces. Some of the 
flaps may be displaced longitudinally 
toward the trailing edge to form supports 
for the aircraft on the ground. 
Impregnated ferrite (2,846,655). A. H. 
Iversen, Santa Monica, Calif., assignor 
to Hughes Aircraft Co. 

Microwave ferrite element impregnated 

with a vitreous material along a _pre- 
determined direction so that the element 
is bonded and hermetically sealed at its 
surface. 
Variable area convergent-divergent ex- 
haust nozzle (2,846,843). D. B. Clark, 
G. 8. Kelley and W. W. Mahnken, Union 
City, N. J., assignors to Curtiss-Wright 
Corp. 

Regulation of the nozzle area so that 
increase in total pressure of the exhaust 
gases and decrease in static pressure of 
the surrounding atmosphere both effect 
an opening adjustment of the nozzle exit. 


Aircraft reaction propulsion units (2,849,- 
861). J. Gardiner and J. 8S. Wallett, 
Allestree, England, assignors to Rolls- 
Royce, Ltd. 

Propulsion nozzle at the outlet of a jet 

pipe porting upstream to auxiliary ducting 
inclined forward to the flight direction to 
an outlet to the atmosphere. Gas flowing 
through the auxiliary duct produces a 
braking effect on the aircraft. 
Variable area thrust deflector-augmenter 
for jet engines (2,846,844). N. W. 
O’Rourke, LaJolla, Calif., assignor to 
The Ryan Aeronautical Co. 

Nozzle having complementary elements 
shiftakle to vary the area of the primary 


orifice. Direction controlling deflectors, 
pivotally mounted, define a secondary 
orifice. The deflectors are connected to 


move the primary and secondary orifices 
in cross-sectional area 


George F. McLaughlin, Contributor 


Fuel system for gas turbine engines 
(2,846,846). F. C. Mock, South Bend, 
Ind., assignor to Bendix Aviation Corp. 

Control of the rate of fuel feed to a g 
turbine having means for avoiding coia- 
pressor instability. Fuel is metered at a 
rate proportional to compressor inlct 
pressure and temperature. Effects of 
temperature on the rate of fuel feed are 
modified in the intermediate speed range 
Zero-length launcher (2,843,020). V. |3 
Bertagna and J. J. Pfarr, Inglewood, 
Calif., assignors to Northrop Aircrait, 
Inc. 

Device for launching a missile. A 
frame structure spaced on a base and per- 
mitting a vehicle supporting a missile to 
advance within a passageway so the frame 
can be lowered to position for launching. 
Angle mounted tip engine for aircraft 
sustaining rotor (2,843,210). R. K. 
Grove, Grosse Pointe Farms, Mich., 
assignor to Studebaker-Packard Corp. 

Tip mounted aerodynamic pod secured 
to a propelling blade and enc losing a gas 
turbine. The exhaust discharge in a jet 
opening in the bottom of the pod is set in 
a direction rearward and slightly down- 
ward so as to produce a thrust having a 
positive blade twisting moment. 
Composite aircraft system and method of 
flight (2,843,337). J. A. J. Bennett, 
San Mateo, Calif., assignor to Hiller 
Helicopters. 

Rotary wing aircraft tug with jet en- 
gines at the rotor tips, and hollow legs 
and mechanism for detachable connection 
to a fixed wing aircraft. Landing of the 
composite aircraft is accomplished by the 
lift of the aircraft tug. 

Circular wing aircraft having pressure 
induced radial airflow (2,843,338). H. 
F. Streib, Chula Vista, Calif. 

An impeller directing a flow of air 
radially inwardly over the upper surface 
of a circular wing to create lift. Air is 
driven downwardly through a central 
circular opening to supplement the lift. 
Jet aircraft (2,843,340). M. H. M. J. 
Wibault, Paris, France. 

Propulsion system comprising a turbo- 
compressor unit obliquely inclined with 
respect to the line of flight. Adjustable 
means permit deflecting the jet discharge 
from normal to a vertical downward 
direction so the jet and air stream pro- 
duce a combined lifting force on the air- 
craft. 

Mobile hydraulic catapulting apparatus 
(2,843,342). EK. F. Ward, Pasadena, 
Calif., assignor to Task Corp. 

Motorized vehicle with portable guide 
track sections adapted for assembly rear- 
ward of a gun for jetting liquid to guide a 
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Rocket engine performance is recorded by oscillograph at this test 


stand high in the Santa Susana Mountains. Test facilities are main- 
tained by Rocketdyne, a division of North American Aviation, Inc. 


Readout: you go 
Why wait for the answers? 


You can immediately evaluate your test results if they 
are recorded on new Kodak Linagraph Direct Print Paper used 
in appropriate moving-mirror galvanometer oscillographs. __ 
There is no need for development: Just run 
your test and read the records as they 
are being made. You'll get sharp, legible traces 
over a wide frequency range—0 to 3000 cps. 
You can make pen or pencil notations on these records. 
And you can make the records permanent by fixation. 
New Kodak Linagraph Direct Print Paper is available in 
5” x 200’, 6” x 100’, 7” x 200’, and 12” x 200’ rolls. Other 
sizes on request. For complete details, write to: 


EASTMAN KODAK COMPANY 


Graphic Reproduction Division 


Rochester 4, N. Y. 
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Reviews 


- Gasdynamics, by A. B. Cambel and B. H. 
Jennings. 


jet impelled carriage for catapulting an 
aircraft. 

Launching device (2,844,073). C. Ré, G. 
_F. Rice, ei F. Birkeland, V. T. Koozin, 
G. C. Duncan and D. Kinda, Torrance, 
] — assignors to Royal Industries, 


“Ripple firing multiple rocket launcher. 
all rockets. Separate series-connected 


pair of adjacent connections with a fused 


Air fuel control for jet 


- (2,849,862). C. H. Jorgensen, J. M. 
_ Barr, H. H. Dietrich, W. T. Nickel and 
M. R. Rowe, Rochester, N. Y., assignors 


to General Motors Corp. 


Fuel control unit ae the control by 


_ the pilot’s lever and also under control of 
i automatic devices which sense certain 
_ parameters affecting the operation of the 
engine. The pilot may change the status 
of the control unit so control by some or 
ail automatic devices is removed. 
Fast rising balloon (2,844,336). H. J. 
ao, Falls Church, Va., assignor 
to the U. S. Navy 
High altitude Ae balloon with an 
inverted conical skirt attached at an angle 
to provide a couple to maintain stability 
during ascent but unrestrained in their 
natural ability to extend. 


Aircraft control arrangement incorporating 
deflectable surface and boundary layer 
control jets (2,844,337). R. B. Mac- 
Arthur and D. F. Swearingen, Hawthorne, 
Calif., assignors to North American 
Aviation, Inc. 

Primary airfoil with a rearward opening 


for discharging air over the upper portion 
of the elevator or aileron when deflected. 


within its trailing edge with nozzles 


A series resistance circuit is provided for 
resistances are provided between each 


circuit between each connection and the 

ground, including a fuse of lesser resistance 

adapted to burn through at a predeter- 
mined time interval. 


Vertical velocity thrust control for hover- 
ing VTOL aircraft (2,844,338). I. W. 
Keith, San Diego, Calif., assignor to 
Ryan "Aeronautical Co. 

Sensitive vertical acceleration sensing 

means and signal and a velocity signal 
connected to an automatic throttle posi- 
tioner with a signal delay means in- 
terposed. 
Improved method and arrangement for 
starting reaction A propulsion units 
(2,840,987). D. Bloomberg and N. 
Burgess, Ohio, 
General Electric Co. 

Auxiliary reaction type propulsion unit 
having low energy starting character- 
istics and adapted to initiate the operation 
of several jet engines. 

Fuel control apparatus for supersonic 
ramjet (2,840,988 ). J.P: 4° West- 
field, N. J., assignor to the U. S. Navy. 

Aerial missile with a body in the form of 
an elongated tubular duct, and containing 
a diffuser and combustor. Static and 
impact pressures in the duct control a 
flow control valve in the fuel feed line. 


Multistage fuel injector for ramjet com- 
bustor (2,840,988). J. P. Longwell, West- 
field, N. J., assignor to the U.S. Navy. 

Set of fuel injection nozzles in an aerial 
missile arranged to inject fuel into air- 
streams passing through venturi throats. 
Nozzles are separated from each other by 
a distance of one-half the wave length of 
the periodic pressure fluctuation in the 
combustor. 


Boundary layer control 2 841,344). M. 
_ Stroukoff, Lawrenceville, N. J. 

Depression in the surface of a wing, 
between the leading and trailing edge 
portions. A series of high velocity jet 
streams of fluid are directed into a pas- 
sageway associated with the depression, 
inducing boundary layer fluid into 
entrance openings fixed to the surface. 

Jet propelled aircraft (2,841,346). W. 
E. Petter, Curdridge, England, assignor 
to Folland Aircraft, Ltd. 

Fairing on the exterior of an engine 
mounted in a fuselage, extending aft 
from the air intake opening, providing a 
compartment for a landing gear unit or 
armament forward of the engine. 


Fastening device (2,841,850). A. T. 
Zohorsky, Altadena, Calif., assignor to 
Aerophysics Development Corp. 

Device for fastening together the ad- 

jacent edges of a material comprising a 
pair of strip members, the inner face of one 
strip in surface engagement with the outer 
face of the other. One strip has a row of 
projections and the other strip a row of 
recesses spaced along its length. 
Radar tracking demonstrating and training 
instrument (2,841,884). H. P. Birming- 
ham, Washington, D. C., assignor to the 
U.S. Navy. 


assignors to 


Circuit for measuring the time duration 
that the deflection of a cathode ray beam 
is maintained within an octogonal area on 
the screen of a cathode ray tube 

Rear inlet annular diffuser (2,841,952). 
L. S. Billman, Lancaster, Calif., assignor 
to the U. 8. Na avy. 

Ramjet aerodyne with a divergeut- 
convergent fuselage. Air entering in- 
lets in the path of travel flows through 
ducts, struts and complementary openi:igs 
to the atmosphere. Air entering air i: let 
ducts may follow an unrestricted path to 
passageways to the fuel ignition mean: in 
the combustion chamber. 


2,860,076 


(2,850,976). 
Calif., as- 


Thrust cancellation device 
H. S. Seifert, Playa del Rey, 
signor to the U. S. Army. 
Device to control the trajectory o! a 
solid propellant rocket without affecting a 
change in the characteristics of the motor 
design or in the propellant powder to 
initially propel the rocket in free flight. A 
neutralizing force overcomes the thrust of | 
the main nozzle, operating with the least 
disturbance or transient force. 
System for pressurizing the fluid propel- 
lant tank of a self-propelled missile 
953). We. Teague, Al- 
pine, N. J., assignor to the U. Army. 
Diaphragm normally sealing a choueee 
way between two compartments of a 
closed pressurizing chamber. A combus- 
tible charge in the first compartment 
moves a piston, ruptures the diaphragm 
and pressurizes the tank. 
Directional control for jets (2,841,955). 
G. H. McLafferty, Manchester, Conn., 
assignor to United Aircraft Corp. 
Converging-diverging diffuser down- 
stream of the air intake, with partitions 
dividing the cooling duct into annular 
segments. Servos selectively close one 
or more valves in each segment to vary 
the flow direction from the main duct. 


Guide vane arrangement for rocket 
missiles (2,842,058). N. E. G. Kuller, 
and K. J. T. Thorildsson, Bofors, Swe- 


den, assignors to Aktiebolaget Bofors Corp. 

Sheet metal vanes attached to a tubular 
rocket by protruding members and slots 
in releasably locking engagement. 


McGraw-Hill Book Co., 
New York, 1958, 353 + xiii pp. text, 


_ pp. tables and 11 pp. index. $11. 


Reviewed by S. KLINE 
Stanford University 


_ This book covers the elements of one- 
_ dimensional gasdynamics, combustion, 
aerothermochemistry and also an intro- 
duction to two- and _ three-dimensional 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


flow problems. It is written at the ad- 
vanced undergraduate level. The cover- 
age and choice of material are excellent, 
and the text is very readable. For 
these reasons, and since no other book 
with such a coverage is available for the 
many courses in this area now developing 
within Mechanical and Aeronautical Engi- 
neering curricula, this book should be of 
considerable assistance to undergraduate 
instructors, and will also undoubtedly 


find a useful place as a reference boc 
some industrial situations. 

A more detailed picture of the mat 
included can be most readily obtained 
from the chapter headings; these are: 
1 Basic Concepts of Gas Dynamics and 
Gas Properties; 2 Fundamental Equa- 
tions of Steady Flew; 3  Isentropi¢ 
Flow; 4 Diabatic Flow; 5 Flow With 
Friction; 6 Wave Phenomena; 7 Vari- 
able-Area Flow; 8 Thermochemistry; 


PrRoputsion 


> 
3 
: 
> 
a 
) e) = 
> 
E 
- 
» 
By 
= 
. 
* 


ration 
beam 
rea on 


1,952). 


signor 


rgeut- 
ig in- 
1rough 
eniigs 
ir inlet 
ath to 
‘ans in 


0,976). 
if., as- 


y ora 
cting a 
motor 
der to 
ht. A 
rust of 
e least 


yropel- 
nissile 
Al 
Army. 
assage- 

of a 
ym bus- 
rtment 
hragm 


1,955). 
Conn., 


down- 
‘titions 
nnular 
one 
O vary 
ct. 

rocket 
Kuller, 
Swe- 
Corp. 
ubular 
d slots 


Editor 


ook in 


aterial 


tained 
e are: 
cs and 
Equa- 
itropic 
y With 
Vari- 
nist ry; 


LSION 


WHAT IT’S LIKE TO BE A CREATIVE ENGINEER AT 


“Test equipment engineering,’’ says John Lloyd, ‘‘is particu- 


TEST EQUIPMENT ENGINEER John W. Lloyd tells why his 
work on the B-70 Weapon System at IBM Owego affords him — 
the creative engineering career he always wanted. es 


larly stimulating when it’s part of a project as new and 
important as the B-70, sometimes described as ‘a huge flying 
computer.’ Right now I’m coordinating the design of engineer- 
ing support equipment for an advanced digital airborne com- 
puter, part of the B-70’s bombing-navigational and missile 
guidance system. There’s a minimum of routine. In order to 
design test equipment you must know — or learn — about the 
equipment to be tested; among these are radar, servo sys- 
tems, digital and analog computers, inertial guidance. | see 
my professional growth assured as IBM continues to develop 
computers for airborne applications.”’ 


Challenging assignments now open 


* Airborne digital & analog computers 
Inertial guidance & missile systems 
Information & network theory 
Magnetic engineering 
Maintainability engineering 

Optics 
Radar electronics & systems 
Servomechanism design & analysis 
Theoretical physics 

Transistor circuits 

® Units & systems test equipment 


There are other openings in related fields to broaden your 
skills and knowledge. 


Qualifications: 


B.S., M.S. or Ph.D. in Electrical or Mechanical Engineering, 
Physics, or Mathematics, and proven ability to assume a high 
degree of technical responsibility in your sphere of interest. 


FOR DETAILS, just write, outlining background and interests, 


to: Mr. P. E. Strohm, Dept. 572Y 
International Business Machines Corp. ° 
Owego, New York 


IBM is a recognized leader in the rapidly expanding electronic 
computer field. Its products are used for both commercial and 
military applications. Continuous growth affords excellent ad- 
vancement opportunities. The ‘‘small-group”’ approach assures 
recognition of individual merit. Company benefits set stand- 
ards for industry, and salaries are commensurate with your 


abilities and experience. 


Plants and laboratories: Endicott, Kingston, Owego, Poughkeepsie, 
Yorktown, N. Y.; Lexington, Ky.; Rochester, Minn.; San Jose, Calif. 
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9 Introduction of Flames and Combus- 
tion; 10 Introduction to Multidimen- 
_ sional Flow; 11 Dimensional Analysis 
and Similitude; 12 Experimental Tech- 
i niques and Measurements; 13 Aero- 
thermochemistry. It should perhaps be 
noted that the authors use the term 
isentropic flow to mean what others have 
called introductory concepts of com- 

_ pressibility, and include such matters as 
the flow in nozzles with variable area in 
_ chapter 7. Chapter 7 also includes a 
__-very brief introduction to the now ca- 
- nonical methods of Shapiro and Hawthorne 
for one-dimensional flow with more than 
one effect occurring simultaneously at 
appreciable Mach numbers. In each 
ease after the first few chapters, the 
- authors carry the work considerably 
forward toward the best available modern 
results. For example, the chapter on 
wave phenomena includes a discussion 
of oblique shock waves and shock wave 
build-up; the chapter on thermochem- 
istry includes discussions reaction 
disassociation, transport phe- 
nomena, ete., and the chapter on dimen- 
sional analysis includes not only the usual 
treatments of Buckingham and Bridgman, 
but the more powerful methods of Ray- 
leigh, Kayser and Huntley, and a brief 
discussion on use of differential equa- 
tions for such problems. Similar ex- 
amples occur in all chapters beyond the 
first three. This reviewer feels that the 
inclusion of the three chapters on com- 
bustion are particularly useful and ap- 
propriate, and these will be doubly wel- 
since the senior author is well 
a known for his many contributions and his 
grasp in this area. The inclusion of the 
material on combustion and on experimen- 
tal methods with the more usual gas- 

_ dynamics should provide a much more 

compact and usable text than has hereto- 
_ fore been available. 
The knowledgeable reader will by now 
undoubtedly have asked himself the 
‘ | es, “How can so much material be 
r placed in a single volume without some 
loss of rigor or completeness?’’ The 
answer, of course, is that it cannot. The 
present book in many instances provides 
only brief introductory treatments par- 
ticularly of the more advanced topics, 
such as multidimensional flow and _ re- 
action chemistry. However, this is 
deemed appropriate for an undergraduate 
textbook since it provides the beginning 
student with at least a picture of the 
scope of present day knowledge. Also, 
by including in most instances the best 
modern references at the end of each 
chapter, the authors are able to make 
reference to appropriate treatments which 
are more complete and rigorous at many 
points, and thus give the student a good 
start on further work in any topics of 
particular interest or concern. 

Also to some extent the concision of 
the present volume is achieved at the 
expense of thorough derivations and the 
sacrifice of a unified fundamental point 
of view. Different methodologies are 
adopted at different points in the text, 
and little is done to associate or rational- 
ize them. The authors have a tendency 
to begin each section or example from the 
view that will lead them to the desired 
answer by the simplest and shortest means 
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for the problem in question. While this 
is permissible and understandable for the 
advanced worker, it may cause some 
difficulty for beginning students who are 
attempting to find for themselves some 
unity and fundamental structure in this 
topic. This reviewer admits he is highly 
biased toward a uniform fundamental 
treatment of the expressions of continuity, 
momentum, energy and the Second Law 
of thermodynamics as the best rational 
means for construction of a uniform 
fundamental basis for gasdynamics, and 
many teachers may find the absence of 
such material in an introductory book far 
less vexatious than the writer. De- 
spite this, the work would be improved, 
in the reviewer’s opinion, by a more com- 
plete derivation of all the fundamental 
equations including more complete de- 
lineation of the conditions and restrictions 
surrounding their use, and by a more sys- 
tematic methodology, both physical and 
mathematical, for introducing these equa- 
tions in the examples given. Examples of 
this problem are the introduction of the 
momentum principle in various forms at 
various times with only a very short and 
incomplete discussion of the relation be- 
tween momentum and Newton’s Law of 
motion as applied to the system of fixed 
mass. In the most general momentum 
equation given (in chapter 10) the 
authors fail to note that the equation 
presented assumes the observer stands on 
the control volume and hence that the 
control volume must be such that its 
coordinates are inertial. Similarly, in 
the primary presentation of the energy 
equation, the discussion is restricted to 
steady, one-dimensional flow of a single 
stream, and the restrictions relating to the 
need for reversibility in the flow work 
term and in the use of an internal energy 
requiring only two independent intensive 
properties are not discussed. Admit- 
tedly, these are not matters that bother 
the average undergraduate, but despite 
this, both the precision of work and the 
information made available to the reader 
would be improved by their inclusion. 
It is also clear that the omission of dis- 
cussion of these matters means that the 
book cannot be considered comparable 
in rigor to the recent advanced discussions 
of this topic, as for example those of 
Shapiro or of Liepmann and Roshko. 

It is not to be inferred from these com- 
ments that the book contains technical 
errors; on the contrary, the results stated 
seem to be remarkably free from techni- 
cal errors and misleading implications. 
Thus the comments offered in the pre- 
vious paragraph concern methodology 
rather than fact and accordingly lie in the 
realm of opinion. The present version of 
the book already represents considerable 
improvement on these matters over an 
earlier draft seen by this reviewer, and the 
comments are offered in the hope that 
possible later editions of this useful work 
will be still further improved in the direc- 
tion of combining maximum rigor and 
logical unity with the very readable, con- 
cise and extensive treatment given. - The 
work would also be improved by inclusion 
of a few examples not entirely dependent 
on perfect gas behavior. 

In summary then, the present volume 
is a contribution and should be of real aid 


in improving the effectiveness of teaching 
in advanced undergraduate courses on 
gasdynamics and combustion. It is not 
as rigorous as the best available mono- 
graphs on these topics, but is nevertheless 
accurate and readable, and it provides a 
surprisingly good coverage of the material 
and introduction to the best modern re- 
sults and references in a single compict 
volume. It is ‘well worth the serious 
attention of anyone interested in a com- 
prehensive introductory treatment in 
these areas for either educational or in- 
dustrial purposes. 


The Russian Literature of Satellites, 
Part I. Edited by A. V. Shpol’skii, 
International Physical Index, Inc., 
New York, 1958, x + 181 pp. $10. 

Reviewed by F. J. KR1tEGER 

The Rand Corporation 


This book is a collection of translations 
of the first six papers (pp. 1-122) of a 
special series published in Uspekhi Fi-i- 
cheskikh Nauk (Progress in the Physical 
Sciences), vol. 63. no. 1a, Sept., 1957, pp. 
1-144, commemorating the centennial of 
the birth of K. E. Tsiolkovskii, founder of 
the science of astronautics. 

In the first paper, ‘Certain Variational 
Problems Associated with the Launching 
of an Artificial Earth Satellite,’ D. 1. 
Okhotsimskii and T. M. Eneev examine 
the problems of: 1 Selecting the opti- 
mum mode of propellant consumption and 
the optimum program for thrust direction; 
2 motion when the mode of propellant 
expenditure is specified; 3 guiding a 
satellite into its orbit while taking into 
account the variation of the gravitational 
field and the rotation of the Earth. 

Using the methods of the calculus of 
variations (assuming a uniform gravita- 
tional field and constant acceleration), 
the authors establish that the most ad- 
vantageous program for rocket thrust 
direction obtains when the tangent of the 
thrust angle is a linear function of time, 
and that the most advantageous mode of 
propellant expenditure is the instantane- 
ous burning of all the propellant at some 
definite moment of flight. 

The motion of a multistage rocket is 
considered in a uniform gravitational 
field (planar Earth), assuming that the 
rocket thrust is constant. The solution 
of the problem is given in analytical form. 
Graphs are included which make it pos- 
sible to determine the various flight char- 
acteristics. 

In an analysis of the more complex 
problem of rocket motion subject to the 
variation of the gravitational field and the 
rotation of the Earth, Okhotsimskii 
Eneev show that the tangent of 
thrust angle, giving the optimum prog! 
for rocket thrust direction, is expressec 
terms of hyperbolic functions of time. 

In the second paper, ‘Determining 
Lifetime of an Artificial Earth Satel 
and Investigating the Secular Pertur 
tions of Its Orbit’? by D. E. Okhots 
skii, T. M. Eneev and G. P. Taratyno 
the variation of atmospheric density w 
altitude is approximated by exponent 
formulas fitted to the data presented in 
S. K. Mitra, “The Upper Atmosphere” 
(2d ed., 1952, p. 582). By assuming that 
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AERO-THERMODYNAMICISTS 
EXPLORE HIGH-SPEED RE-ENTRY 


Advanced weapon system technology has brought to the forefront problem areas 
requiring attention to interaction between aerodynamic and thermodynamic “a 
_ phenomena. Typical of these is the problem of high-speed atmospheric re-entry. 


A report to Engineers Expanding research and development activities have coincided with acceleration 
ees _ on top priority programs like our Polaris IRBM. At the same time, positions 
and Scientists from for qualified engineers and scientists have opened up that are unequalled in 
Lockheed Missile Systems — _ responsibility or in opportunities for moving ahead. 


Positions in aero-thermodynamics include such areas as: aerodynamic 

characteristics of missiles at high Mach numbers; missile and weapon system 

programs insure more design analysis; boundary layer and heat transfer analyses in hypersonic flow 
promising careers fields; and calculation of transient structural and equipment temperatures 


where expanding missile 


resulting from aerodynamic heating and radiation. 


In addition, openings exist at all levels in Gas Dynamics, Structures, Propulsion, 
Test Planning and Analysis, Test Operations, Information Processing, Electronics, 
and Systems Integration. For these and other positions, qualified engineers 

and scientists are invited to write Research and Development Staff, 

Dept. 2511, 962 W. El Camino Real, Sunnyvale, California. —_ 


Lochheed 


SUNNYVALE, PALO ALTO, VAN NUYS, SANTA CRUZ, VANDENBERG AFB, CALIFORNIA 
CAPE CANAVERAL, FLORIDA * ALAMOGORDO, NEW MEXICO 


MISSILE SYSTEMS DIVISION 


Maurice Tucker, Aero-Thermodynamics Department Manager, 
discusses combined aero-thermodynamic re-entry body tests being 
conducted in Division's new “hot-shot” wind tunnel. Others are 

Dr. Jerome L. Fox, Assistant Department Manager, Thermodynamics, le 
and Robert L. Nelson, Assistant Department Manager, Aerodynamics. 
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the Earth’s gravitational field is central 
(i.e., the Earth is spherical) and neglecting 
the diurnal rotation of the atmosphere, 
the authors reduce the problem of de- 
termining the lifetime of a satellite to the 
integration of a system of two differential 
equations, the right-hand sides of which 
are expressed in terms of definite integrals 
whose integrand functions are known 
functions of the semilatus rectum, the 
eccentricity and the argument of latitude 
(i.e., the sum of the true anomaly and 
the argument of perigee). These func- 
tions are also independent of the aero- 
dynamic drag coefficient c; and of the 
satellite’s design parameters (i.e., weight 
G and maximum cross-sectional area F). 
_ The computational results, presented 
in both graphical and tabular form, give 
values for a quantity » and for the initial 
perigee velocity as functions of the peri- 
gee and apogee altitudes. Satellite life- 
time—in terms of the number of revolu- 
tions N—is obtained from the formula 
N = vG/Fgcr, where g is the acceleration 
of gravity. 

An analysis of the secular perturbations 
of orbital elements caused by other per- 
turbing forces (e.g., the oblateness of the 
Earth and the diurnal rotation of the 
atmosphere) indicates that the pertur- 
bations are small and “that the method 
used in this paper to compute the lifetime 
is a sufficiently correct one.” 

_It is interesting to compare the com- 
‘ puted lifetime of a satellite with its actual 


lifetime. Sputnik I data published by 
TASS—perigee altitude 228 km, apogee 
altitude 947 km, diameter 58 cm, weight 
83.6 kg—and a drag coefficient equal to 
2 show that » = 264 and N = 4260 
revolutions. This value of N is consider- 
ably higher than Sputnik I’s actual 
lifetime of about 1400 revolutions or 92 
_ days, indicating that the upper atmos- 
_ phere is somewhat denser than heretofore 
supposed. 
‘The third paper, “The Motion of an 
_ Artificial Earth Satellite in the Non- 
- Gravitational Field of the Earth When 
_ Atmospheric Resistance Is Present’’ by 
G. P. Taratynova, and the fourth paper, 
_ “The Effect of Geophysical Factors on 
_ the Motion of a Satellite’. by I. M. 
& Iatsunskii, are ancillary, and in a sense 
oe complementary, to the second paper. 
_ They both take into consideration the 
perturbing forces caused by the 
_ oblateness of the Earth and by its atmos- 
: pheric drag, including the diurnal rotation 
_ of the atmosphere with the Earth. Iat- 
_ sunskii’s paper discusses the possibility 
of determining the density of the air, the 
~ constants of the Earth ellipsoid and gravi- 
- tational anomalies by means of satellite 
observations. 
- Although not concerned with satellites, 
— the piece de résistance of the collection is 
A. Egorov’s paper, “Certain Problems 
This paper 
is a result of a systematic investigation 
undertaken at the Steklov Mathematics 
Institute of the USSR Academy of 


4 


satisfactory solutions for the funda- 
mental problems of the theory of flight 
dynamics, including the shape and classi- 
fication of trajectories in the passive 
_ (unpowered) phase; possible trajectories 
Se for cireumflight of the moon and return 


to the Earth; the minimum initial veloc- 
ities required for reaching the moon; the 
problem of hitting the moon, and the very 
important problem of the effect of dis- 
persion in the initial conditions on the 
characteristics of various moon-flight 
trajectories. 

At the outset, Egorov reduces the prob- 
lem of rocket flight to the moon to the 
circular restricted three-body problem by 
neglecting the oblateness of the Earth and 
the perturbations of the sun and the large 
planets. He also assumes that the rock- 
et’s flight trajectory begins outside the 
Earth’s atmosphere at a height of several 
hundred kilometers. 

The problem of minimum velocities 
Egorov solves by means of Jacobi’s in- 
tegral. He finds that the minimum 
velocity required to reach the moon (from 
an altitude of 200 km above the Earth’s 
surface) is 10.84890 km/sec, and that 
required for a rocket to proceed to infinity 
is 10.84968 km/sec. 

To determine minimum-velocity trajec- 
tories Egorov uses the method of nu- 
merical integration. These trajectories, 
although interesting academically, are of 
no practical interest because they remain 
approximate ellipses with a focus at the 
Earth’s center for a long time, so that the 
rocket must make a large number of 
revolutions around the Earth (hundreds 
or more) before reaching the moon. Of 
more interest is the author’s demon- 
stration that the moon cannot possibly 
capture a projectile launched from the 
Earth on the first circuit of the trajectory, 
no matter what values of the initial con- 
ditions are specified. 

Egorov describes a simple but accurate 
approximate method of analyzing ap- 
proach trajectories based on the concept 
of separate “spheres of action’’ for the 
Earth and for the moon. Thus, it can 
be assumed the motion outside the moon’s 
sphere of action occurs along a geocentric 
conic section, and the motion inside its 
sphere of action occurs along a seleno- 
centric conic section. This method shows 
that motion relative to the moon along 
approach trajectories always takes place 
at velocities greatly exceeding the para- 
bolic selenocentric velocities. Egorov also 
analyzes the evolution of a complete set of 
approach trajectories when the magnitude 
and direction of the initial velocity are 
varied. 

In solving the problem of striking the 
moon, Egorov develops a method for 
determining the initial conditions that 
correspond to hitting the moon with a 
predetermined accuracy. He establishes 
that the deviation of the trajectory from 
the center of the moon is a quadratic, 
not a linear, function of small errors in 
the initial conditions. Similarly, in the 
problem of circumflight of the moon, he 
develops a method for determining the 
initial conditions for a return to Earth. 
Here the deviation of the return-flight 
trajectory from the center of the Earth is 
a quadratic function of the errors in the 
initial conditions. It is found that at a 
distance greater than the sphere of action 
of the moon (66,000 km), the accuracies 
required for a return to Earth would be no 
greater than in the problem of hitting the 
moon for the same initial velocity. Thus 
the problems of: 1 Hitting the moon 


and, 2 circumflight of the moon at a 
sufficiently great distance and returning 
to Earth without any correction during 
the passive (unpowered) phase of the 
motion are evidently — technologically 
possible. 

Egorov also examines three other 
problems: 1 The special problem of 
lunar circumflight, i.e., the determination 
of circumflight -trajectories by which the 
rocket enters the Earth’s atmosphere «it a 
shallow angle; 2 the problem of whether 
or not it is possible to launch a projectile 
into the plane of the lunar orbit in such a 
way that it will periodically circumfly the 
moon, return to the vicinity of the Ea:th, 
and regularly give information about the 
moon (e.g., about its opposite side), this 
problem has not been previously treated in 
the literature; 3 the problem of accel- 
erating or decelerating a cosmic rocket by 
means of the moon without making use of 
engine thrust. 

In his paper, Egorov draws the jol- 
lowing conclusion: For close circumflight 
of the moon (and especially for a shaliow 
re-entry into the Earth’s atmosphere) 
the accuracy requirements in the initial 
conditions are so stringent that a cor- 
rection of the passive phase of the trajec- 
tory definitely seems necessary. For tra- 
jectories that hit the moon and for those 
that achieve distant circumflight of the 
moon, the accuracy requirements are 
moderate. Thus, it will be possible to 
achieve flights along these trajectories, 
without introducing any correction during 
the passive phase, as soon as rockets can 
achieve velocities of the order of the 
parabolic velocity. 

The sixth paper, “The Use of Artificial 
Earth Satellites for Verifying the General 
Relativity Theory” by V. L. Ginzburg, 
is simply a résumé of a more comprehen- 
sive article published in Uspekhi Fizi- 
cheskikh Nauk (vol. 59, no. 1, 1956, pp. 
11-49). The following four effects, which 
emerge from the general theory of rela- 
tivity, are considered: 1 The rotation of 
a planet’s perihelion around the sun; 2 
the bending of light rays that pass close 
to the sun; 3 the gravitational shift of 
spectral lines; 4 the displacement of a 
planet’s perihelion caused by the rotation 
of the sun. Ginzburg discusses the 
observational data obtained by astro- 
nomical methods to give a quantitative 
verification of the theory. He indicates 
how the use of artificial satellites will make 
it possible to verify the general theory of 
relativity (i.e., measure relativistic ef- 
fects) more rapidly and with greater 
accuracy than heretofore. 

In general, the quality of translation of 
the papers comprising the book is good. 
The editors have corrected numerous 
typographical errors in the original text 
but have made a few slips of their own, 
such as the persistent use of the word 
“perihelion” instead of “perigee” and of 
the redundant expression “conic cross 
section.”” In spite of its price (due to the 
enormous effort and expense of issuing 
such a translation), the book will un- 
doubtedly be welcomed by the growing 
number of technical people who are se- 
riously interested in problems involving 
the motion of a satellite (or lunar vehicle) 
and the scientific research that can be 
accomplished with it. 
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New extreme-high-temperature lubricants for missiles 


and supersonic aircraft SHELL ETR GREASES 
One of the serious lubricating problems _ tures. Co-operation with representatives special thickener—an organic vat dye— 
faced by designers of missiles and super- of bearing manufacturers and military _ which has exceptional heat stability and 
sonic aircraft has been solved by scientists personnel resulted in a completely new _jelling efficiency. 
at Shell Research Laboratories. class of greases—SHELL ETR GREASES. If you are presently in the market for 
The problem: to find a grease which These greases can easily withstand tem- _—an_ultra-high-temperature-range grease, 

would permit components to operate with peratures up to 600°F. They give superior _ we will be glad to provide more informa- 
certainty under extreme high tempera- _lubricating performance because of a tion on Shell ETR Greases. 


 §HELL OIL COMPANY 
‘was, 50 West SOth Street, New York 20,N.Y. 
ieee 100 Bush Street, San Francisco 6, Calif. ae 
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ENGINEERS 
& SCIENTISTS 


Here is your opportunity to grow 
with a young, expanding subsidiary 
of the Ford Motor Company. Out- 
standing career opportunities are 
open in Aeronutronic’s new RE- 
SEARCH CENTER, overlooking the 
Pacific at Newport Beach, and the 
facility in Glendale, California. 
You will have all the advantages of a 
stimulating mental environment, 
working with advanced equipment 
in a new facility, located where you 
can enjoy California living at its 
finest. 


PhD and MS RESEARCH SPECIAL- 
ISTS with 5 to 7 years’ experience in heat 
transfer, fluid mechanics, thermodynamics, 
combustion and chemical kinetics, an 

thermoelasticity. To work on theoretical 
and experimental programs related to re- 
entry technology and advanced rocket 
propulsion. Specific assignments are open 
in re-entry body design, high temperature 
materials studies, boundary layer heat trans- 
fer with chemical reaction, thermal stress 
analysis, and high temperature thermo- 


dynamics. 

PROPULSION ENGINEERS with 5 years’ 
experience in liquid and solid rocket design 
and test. Familiarity with heat transfer 
problems in engines desirable. To work on 
program of wide scope in R & D of advanced 
concepts in rocket _— components, and 
for missile project w 

ADVANCED AERODYNAMIC FACIL- 
ITY DESIGNER. Advanced degree de- 
sired. To supervise work in design and in 
instrumentation of advanced aerodynamic 
test facilities such as shock tubes, shock 
tunnels, plasma-jets, and hyper-velocity 


guns. 
STRUCTURAL ANALYSIS SECTION 
SUPERVISOR with 8 to 10 years’ experi- 
ence, including supervision, in the missile 
field. Graduate degree for design and 
required. Will be to 
apply k ledge of high t ure ma- 
terials and methods, + le stress, dynam- 
etc. to advanced vehicles, 
e-entry bodies, and space vehicl 
FLIGHT TEST & IN: STRUMENTATION 
ENGINEERS with 5 to 10 years’ experi- 
ence in laboratory and flight test instru- 
mentation techniques. Will develop tech- 
niques utilizing advanced instrumentation 
associated with space vehicles. 
THEORETICAL AERODYNAMICIST. 
Advanced degree and at least 5 years’ ex- 
perience in high-speed aerodynamics. 
Knowledge of viscid and inviscid gas flows 
required. To work on program leading to 
advanced missile configurations. ork in- 
volves analysis of the re-entry of hypersonic 
missiles and space craft for determining 
optimum configuration. 
DYNAMICIST. Advanced degree, ap- 
plied mathematics background, and ex- 
perience in missile stability analysis de- 
sirable. Work involves re-entry dynamics 
of advanced vehicles and dynamic analysis 
of ice craft. 
ENGINEER or PHYSICIST. With ex- 
perience in the use of scientific instru- 
ments for making physical measurement. 
Work related to flight test and facility 
instrumentation. Advanced degree de- 
sired with minimum of 3 years of related 
experience. 


Qualified applicants are invited to send 
resumes and inquiries to Mr. L. R. Stapel. 


AERONUTRONIC SYSTEMS, 
INC. 


A subsidiary of Ford Motor Company 


1234 Air Way, Bidg. 19, Glendale, Calif. 
CHapman 5-6651 
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The James Forrestal Research Center, Princeton University 


Jet and Rocket 


Propulsion Engines 


Replacement of Piston Engines by Gas 
Turbines in Air Liners, by H. Sammons, J. 
Roy. Aeron. Soc., vol. 62, no. 556, Feb. 
1958, pp. 94-104. 

Effect of Fuel Density and Heating 
Value on Ram-jet Airplane Range, by 
Hugh M. Henneberry, NACA Res. Mem. 
E51121, Feb. 1952, 56 pp., diagrs., 3 
tabs. (Declassified from Confidential by 
authority of NACA Res. Abstr. 123, p. 8, 
1/7/58.) 

Experimental Investigation of Turbojet- 
engine Multiple-loop Controls for Non- 
afterburning and Afterburning Modes of 
Engine Operation, by Donald B. Dirsch, 
Leon M. Wenzel and Clint E. Hart, NACA 
TN 4159, Jan. 1958, 61 pp., diagrs., tab. 

Fundamental Rocket Research. Sum- 
mary Report Covering the Period July 1, 
1954, to July 1, 1956, by American Ma- 
chine & Foundry Co., Turbo Div., Rep. GA 
104 (AFOSR TR-56-55; ASTIA AD 
110373), Aug. 1956, 64 pp., 26 figs. 

Vapor Fuel Distribution Effects on Com- 
bustion Performance of a Single Tubular 
Combustor, by Richard J. McCafferty, 
NACA Res. Mem. E50J03, Dec. 1950, 27 
pp. (Declassified from Confidential by 
authority of NACA Res. Abstr. 123, p. 7, 
1/7/58.) 

Preliminary Investigation of an Annular 
Turbojet Combustor Having a Catalytic- 
coated Linear, by Carl T. Norgren and J. 
Howard Childs, NACA Res. Mem. E- 
53L07, Jan. 1954, 11 pp., diagrs. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 123, p. 10, 1/7/58.) 

Solid Propellant Rocket Motors, by 
E. T. B. Smith, B.Se., D.C.Ae., A.M.I. 
Mech. E., A.F.R. Ae.S., J. Brit. Inter- 
planet. Soc., vol. 16, Oct.-Dec. 1957, pp. 
198-211. 

The Calculation of Fuel Distribution in 
Step-rockets, by M. L. Williams, JJ. 
Brit. Interplanet. Soc., vol. 16, Oct.—Dec. 
1957, pp. 211-215. 

Evaluation of an Automatic Inlet-pres- 
sure Control Valve for Study of Transient 
Engine Performance Characteristics, by 
Lewis E. Wellner, Robert J. Lubick and 
Harry E. Bloomer, NACA Res. Mem. 
E55L13, April 1956, 25 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 124, p. 44, 2/10/58.) 

Altitude Wind Tunnel Investigation of 
Tail Pipe Burner with Converging Conical 
Burner Section of J35-A-5 Turbojet En- 
gine, by H. Carl Thorman and Carl E. 
Campbell, NACA Res Mem. E9116, Feb. 
1950, 60 pp. (Declassified from Confiden- 
tial by authority of NACA Res. Abstr. 
124, p. 40, 2/10/58.) 

Altitude Performance and Operational 
Characteristics of 29 Inch Diameter Tail 
Pipe Burner with Several Fuel Systems 
and Fuel Cooled Stage Type Flame Hold- 
ers on J35-A-5 Turbojet Engine, by Rich- 
ard L. Golladay and Harry E. Bloomer, 
NACA Res. Mem. E50A19, April 1950, 57 
pp. (Declassified from Confidential by 


Ep1tor’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the Heat 
Transfer Laboratory, University of Minne- 
sota, are gratefully acknowledged. 


VACA Res. Abstr. 124, p 
, 2/10/58.) 

Effect of Diffuser Design, Diffuser-exit 
Velocity Profile, and Fuel Distribution of 
Altitude Performance of Several After- 
burner Configurations, by E. William Con- 
rad, Frederick W. Schulze and Karl H. 
Usow, NACA Res. Mem. E53A30, J:ly 
1953, 63 pp. (Declassified from Confid«n- 
tial by authority of NACA Res. Abstr. 
124, p. 42, 2/10/58.) 

Investigation of a Supersonic-compres- 
sor Rotor with Turning to Axial Direction, 
I—Rotor Design and Performance, by 
Edward R. Tysl, John F. Klapproth «nd 
Melvin J. Hartmann, NACA Res. Mem. 
E53F23, Aug. 1953, 36 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 124, p. 42, 2/10/58.) 

Investigation of a Supersonic-compres- 
sor Rotor with Turning to Axial Direction, 
II—Rotor Component Off-design and 
Stage Performance, by Melvin J. Hart- 
mann and Edward R. Tysl, NACA Res. 
Mem. E53L24, March 1954, 31 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 124, p. 42, 2/10/5s.) 

Design and Test of Mixed-flow Im- 
pellers, I1I—Design and Experimental 
Results for Impeller Model MFI-2A and 
Comparison with Impeller Model MFI- 
1A, by Joseph T. Hamrick, Walter \M. 
Osborn and William L. Beede, NACA 
Res. Mem. E52L22a, March 1953, 34 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 124, p. 41, 
2/10/58.) 

Design and Test of Mixed-flow Im- 
pellers, IV—Experimental Results for 
Impeller Models MFI-1 and MFI-2 with 
Changes in Blade Height, by Joseph T. 
Hamrick, William L. Beede and Joseph 
R. Withee Jr., NACA Res. Mem. E53L02, 
Feb. 1954, 32 pp. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 124, p. 42, 2/10/58.) 

Design and Test of Mixed-flow Im- 
pellers, V—Design Procedure and Per- 
formance Results for Two Vaned Dif- 
fusers Tested with Impeller Model MFI- 
1B, by Joseph T. Hamrick and Walter M. 
Osborn, NACA Res. Mem. E55E13, July 
1955, 32 pp. (Declassified from Conf- 
dential by authority of NACA Res. Abstr. 
124, p. 43, 2/10/58.) 

Design and Test of Mixed-flow Impeli- 
lers, VI—Performance of Parabolic-bladed 
Impeller with Shroud Redesigned by 
Rapid Approximate Method, by Kenneth 
J. Smith and Walter M. Osborn, NACA 
Res. Mem. E55F23, Sept. 1955, 22 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 124, p. 44, 
2/10/58.) 

Design and Test of Mixed-flow Impel- 
lers, VII—Experimental Results for Para- 
bolic-bladed Impeller with Alternate 
Blades Cut Back to Form Splitter Vanes, 


by Walter M. Osborn, NACA Res. Mem. 
K55L15, March 1956, 13 pp. (Declass i- 
fied from Confidential by authority of 


NACA Res. Abstr. 124, p. 44, 2/10/58. ) 
Design and Test of Mixed-flow Impel- 
lers, VIII—Comparison of Experimental 
Results for Three Impellers with Shroud 
Redesigned by Rapid Approximate 
Method, by Walter M. Osborn, Kenneth 
J. Smith and Joseph T. Hamrick, NACA 
Res. Mem. E56L07, Feb. 1957, 32 pp. (De- 
classified from Confidential by author- 
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Application of Radial-equilibrium Con- 
dition to Axial-flow Turbomachine Design 
Including Consideration of Change of En- 
tropy with Radius Downstream of Blade 
Row, by James E. Hatch, Charles C. Gia- 
mati and Robert J. Jackson, NACA Res. 
Mem. E54A20, April 1954, 42 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 124, p. 43, 2/10/58. ) 

Performance of a 16-inch Ram-jet En- 
gine with a Can-type Combustor at Mach 
Numbers of 1.50 to 2.16, by Donald P. 


124, p. 44, 


Hearth and Eugene Perchonok, NACA 
Res. Mem. E54G13, Aug. 1954, 30 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 124, p. 43, 


2/10/58.) 

Effects of Inlet-air-flow Distortion on 
Steady-state Altitude Performance of an 
Axial-flow Turbojet Engine, by E(arl) 
William Conrad, Morgan P. Hanson and 
John E. MecAulay, NACA Res. Mem. 
E55A04, Sept. 1955, 47 pp. ( Declassified 
from Confidential by authority of NACA 
Res. Abstr. 124, p. 48, 2/10/58.) 

Power Packages for Medium-sized 
Missiles, by A. E. Maine and E. Wall, 
Missiles and Rockets, vol. 3, no. 2, Feb. 
1958, pp. 115-120. 


Aerodynamics of Jet 
Propelled Vehicles 


On the Numerical Calculation of De- 
tached Bow Shock Waves in Hypersonic 
Flow, by P. R. Garabedian and H. M. 
Lieberstein, J. Aeron. Sci., vol. 25, no. 2, 
Feb. 1958, pp. 109-116. 

Estimating Aerodynamic Characteristic 
Times in Hypersonic Flow, by Sinclaire 
M. Seala, J. Aeron. Sci., vol. 25, no. 2, 
Feb. 1958, pp. 131-132. 

A New Method for Computing Drag 
Coefficients from Ballistic Range Data, by 
Alvin Seiff, J. Aeron. Sci., vol. 25, no. 2, 
Feb. 1958, p. 133. 

Stability of Bodies of Revolution Having 
Fineness Ratios Smaller Than 1.0 and 
Having Rounded Fronts and Blunt Bases, 
by Stanley H. Scher and James 8. Bow- 
man Jr., NACA Res. Mem. L52L08, Jan. 
1953, 23 pp., diagrs., 5 tabs. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abstr. 123, p. 14, 1/7/58.) 

A Flight Technique for the Measure- 
ment of Thrust Boundaries and of Drag 
Due to Lift, by H. D. Rylands, North 
Atlantic Treaty Org., Advisory Group for 
Aeron. Res. and Dev., Rep. 123, May 
1957, 14 pp. 

Aerodynamic Characteristics of a Ro- 
tating Model of the 5.0 Inch Spin Sta- 
bilized Rocket Model 32, by R. H. Heald, 
J.G. Logan Jr., H. Spivak and W. Squire, 
Bur. Ordn., TN 37, June 1957, 18 pp. 

Summary of Flutter Experiences as a 
Guide to the Preliminary Design of Lifting 
Surfaces on Missiles, by Dennis J. Mar- 
tin, NACA TN 4197, Feb. 1958, 21 pp. 
(Supersedes Res. Mem. L51J30. ) 


Heat Transfer and 
Fluid Flow 


Heat Transfer in the Condensation of 
Metal Vapors: Mercury and Sodium up to 
Atmospheric Pressure, by Balabhadra 
Misra and Charles F. Bonilla, Chem. 
Engng. Progr., vol. 53, 1957. 

Effect of Agitation on the Critical Tem- 
perature Difference for a Boiling Liquid, 
by F. S. Pramuk and J. W. Westwater, 


NOVEMBER 1958 


TYPICAL APPLICATIONS of 
thread. weave it through or around the propellant 


string, 


New Du Pont Pyrocore’” 


core’ Ipnite L d 
Pyrocor Ipniter. You can 


TRADEMARK 


... provides faster, more uniform ignition , 
of rocket propellants 


“Pyrocore” is a small diameter flexible 
metal tubing, available in practically 
any length you desire, containing a 
detonating-ignition composition as its 
core. 

Propellant ignition may be effected 
along its length at velocities ranging 
from 12,000 to 21,000 feet per second, 
depending upon the type and amount 
of core composition you specify. 

In laboratory tests, “Pyrocore” cut 
ignition time 99.5% when compared 
with a standard primer. Total ignition 
of a 2114” length of cannon primer was 
achieved in less than 4 millisecond, as 
against 50 milliseconds required with 
a squib primer. 


Greater Design Freedom 


The stringlike form of ‘“‘Pyrocore” gives 
you new freedom in designing rocket 
motors. For example, you can thread 
“Pyrocore” through jelly roll or basket 
type igniter assemblies, or string it 
around the propellant grains. Or you 
can weave coils or loops of ‘“‘Pyrocore” 
directly into the propellant grain either 
during or after grain formation. 
Further, “Pyrocore” has very low 


PONT “PYROCORE” IGNITER 


TRADEMARK 


PRODUCT OF DU PONT RESEARCH 7 


REG. 5. PAT. OFF 


Better Things for Better Living .. 


brisance (shattering effect), and can 


convey a non-electric stimulus safely _ 
past expensive equipment and bulk- | 


heads. 


Initiator assemblies may be located — 
outside the propellant chamber, elimi- | 


nating a major source of residue in the 
reaction zone. The “Pyrocore’’ itself 
leaves very little residue. The portion 
of “Pyrocore” not actually used for 
ignition may be easily insulated with 
any one of a variety of plastic com- 
pounds. 


Sample Kit Available 


Sample kits, con- 
taining several 
variants of the 
“Pyrocore” Ignit- 
er and various 
primers and ac- 
cessories, are 
available from 
DuPont. Kit 
prices and _tech- 
nical bulletins can 
be obtained by writing to E. I. du Pont 
de Nemours & Co. (Inc.), Explosives 
Department, Wilmington 98, Del. 
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Correlation of Maximum-heat-flux Data 
for Boiling of Saturated Liquids, by W. 
. and P. Griffith, Chem. Engng. 
Progr., vol. 53, 1957. 
Some Experiments Related to the 
Noise from Boundary Layers, by Harvey 
—H. Hubbard, J. Acoust. Soc. America, 
29, 1957. 
he: Investigation of Burnout Heat Flux in 
Rectangular Channels at 2000 Psia, by 
H. S. Jacket, J. D. Roarty and J. E. Zerbe, 
Wi estinghouse Electric Corp., Atomic Power 
Div., (U.S. Atomic Energy Comm.), 
WAPD-A2W (IM)-3, Dec. 1955. 

Review of Boiling Heat Transfer with 
Particular Reference to Unstable Flow, by 
J. MeNeillis, Engng., vol. 183, 1957, p. 

686. 


Unsteady-state Heat Transfer to and 
from Gases in Laminar Flow, by T. D. 
Patten, Chartered Mech. Engng., vol. 4, 
no. 6, June 1957. 

Radial Temperature Distribution in an 
Annular Coolant Passage, by Henry 
Barrow, The Engr., May 10, 1957, p. 713. 

Convective Heat Transfer from High 
Temperature Air Inside a Tube, by H. E. 
Zellnick and S. W. Churchill, ASME 
Paper 57-HT-33 (ASME- AIChE Heat 
Transfer Conference, University Park, Pa., 
Aug. 11-15, 1957), Aug. 1957, 8 pp. 

Gas Friction-Heat Transfer Charts for 
Ducted Flows, by S. V. Manson, ASME 
Paper 57-HT-34 (ASME- AIChE Heat 
Transfer Conference, University Park, Pa., 
Aug. 11-15, 1957), ‘Aug. 1957, 3 pp. 

Theoretical and Experimental Investi- 
gation of a Two Dimensional All Super- 
sonic Diffusing Cascade at Mach Number 
2.5, by E. Beder, Propulsion Res. Corp., 
R-155, June 1955, 126 pp. 

Systematic Investigations on Secondary 
Flow Losses in Cascades, by N. Scholz, 
Braunschweig, Techn. Hochschule, Inst. 
fiir Stromungsmechanik, Rep. 56/18a 
(AFOSR-T R-57-39; ASTIA A D 132400), 
March 1957, 61 pp. 

Thermal Conductivity of Fluids: Nitric 
Oxide, by G. N. Richter and B. H. Sage, 
Project Squid, Tech. Rep. CIT-2-P 
(ASTIA AD 139905), Aug. 1957, 18 pp. 
(available only on microcard). 

Thermal Conductivity of Pure Liquids, 


by W. J. Scheffy and E. F. Johnson, 
Project Squid, Tech. Rep. PR-71-M 
(ASTIA AD 142428), Aug. 1957, 19 pp. 


(available only on microcard). 

Initial Inclination of the Mixing Bound- 
ary Separating an Exhausting Supersonic 
Jet from a Supersonic Ambient Stream, 
by Eugene 8S. Love, NACA RM L55J14, 
Jan. 1956, 30 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
121, p. 27, 11/4/57). 

Preliminary Investigation of Perform- 
ance of Variable-throat Extended-plug- 
type Nozzles over Wide Range of Nozzle 
Pressure Ratios, by Carl C. Ciepluch, H. 
George Krull and Fred W. Steffen, NACA 
RM E53J28, Feb. 1954, 35 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs. 121, p. 18, 11/4/57.) 

Some Studies of Axisymmetric Free 
Jets Exhausting from Sonic and Super- 
sonic Nozzles into Still Air and into Super- 
sonic Streams, by Eugene S. Love and 
Carl E. Grigsby, NACA RM L54L31, 
May 1955, 178 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
121, p. 26, 11/4/57.) 

Film Boiling from Vertical Tubes, by 
Y. Y. Hsu and J. W. Westwater, ASME 
Paper 57-HT-24 (ASME-AIChE Heat 
Transfer Conference, University Park, Pa., 
Aug. 11-15, 1957), Aug. 1957, 8 pp. 
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The Correlation of Nucleate Boiling 
Burnout Data, by P. Griffith, ASME 
Paper 57-HT-21 (ASME-AIChE Heat 
Transfer Conference, University Park, Pa., 
Aug. 11-15, 1957), Aug. 1957, 17 pp. 

On the Stability of Boiling Heat Trans- 
fer, by Novak Zuber, ASME Paper 
57-HT-4, (ASME-AIChE Heat Transfer 
Conference, University Park, Pa., Aug. 
11-15, 1957), Aug. 1957, 4 pp. 

Bubble Growth Rates in Boiling, by 
Peter Griffith, ASME Paper 57-HT-2 
(ASME-AIChE Heat Transfer Conference, 
University Park, Pa., Aug. 11-15, 1957), 
Aug. 1957, 6 pp. 

Description of Boiling Project Burnout 
Detector, by M. W. Raymond and J. M. 
Reynolds, Mass. Inst. Tech., Heat Trans- 
fer Lab., Tech. Rep. 11, July 1957, 6 pp. 

Burnout :n Forced Convection Nucleate 
Boiling in Water, by J. M Reynolds, Mass. 
Inst. Tech., Heat Transfer Lab., Tech. 
Rep. 10, July 1957, 17 pp. 

The Correlation of Nucleate Burnout 
Data, by Peter Griffith, Mass. Inst. Tech., 
Tech. Rep. 9, March 1957, 15 pp., 2 figs. 

Some Proposals Regarding the Defini- 
tions of Terms Relating to Various Flow 
Regimes of a Gas, by C. H. E. Warren 
and A. D. Young, Gi. Brit., Aeron. Res. 
Council, Curr. Paper 368 (formerly A.R.C. 
Tech. Rep. 19583), 1958, 6 pp. 1 fig. 

Radiative Properties of High Tempera- 
ture Air, by W. H. Wurster, H. S. Glick 
and C. E. Treanor, Cornell Aeron. Lab., 
Inc., Rep. QM-997-A-1, Sept. 1957, 28 pp., 
15 figs. 

Limited Investigation of Noise Suppres- 
sion by Injection of Water into Exhaust of 
Afterburning Jet Engine, by Max C. 
Kurbjun, NACA Res. Mem. L57L05, Feb. 
1958, 15 pp. 

Theoretical Investigation of Subsonic 
Oscillatory Blade-row Aerodynamics, by 
Frank Lane and Manfred Friedman, 
NACA TN 4136, Feb. 1958, 64 pp. 

The Effect of Inlet Circumferential 
Maldistribution on an Axial Compressor 
Stage, by R. C. Turner, J. Ritchie and C. 
E. Moss, Gt. Brit., Aeron. Res. Council, 
Rep. & Mem. 3066 (formerly A.R.C 
Tech. Rep. 19047; Natl. Gas Turbine 
Estab., Rep. R. 203), 1958, 16 pp. 

Thermodynamic and Transport Prop- 
erties of Dissociated Hydrogen Mixtures, 
by Martin J. Reisfeld, Atomic Energy 
Comm., LA-2123, June 1957, 23 pp. 


Heat Transfer in Manhattan District 
and Atomic Energy Commission Labora- 
tories: A Critical Survey, by W. B. Har- 
rison, Atomic Energy Comm., ORNL-156, 
Oct. 1948, 255 pp. 


Use of a Free Molecular Probe in High 
Speed Rarefied Gas Flow Studies, by J. 
A. Laurmann and D. C. Ipsen, Wright 
Air Dev. Center, Tech. Rep. 57-440 (AS- 
TIA AD 142109), Oct. 1957, 40 pp. 

An Instrument for the Determination of 
Instantaneous Densities in a Flowing Gas, 
by Henry Lopez, Rensselaer Polytech. 
Inst., Res. D'v., TR AE 5708, Jan. 1957, 
49 pp. 

The Thermal Conductivity of Nitrogen 
and — in the Liquid and Gaseous 
States, by H. Ziebland and J. T. A. Bur- 
ton, British J. of Appld. Phys., 
2, Feb. 1958, pp. 52-59. 


An Investigation of Droplet Oscillation 
During Mass Transfer,I. The Conditions 
Necessary, and the Source of a Energy 
for the Oscillations, by D. A. Haydon, 
Proc. Royal Soc., Series a. ca 243, no. 
1235, Feb. 11, 1958, pp. 483-491. 


An Investigation of Droplet Oscillation 
Dong Mass Transfer, II. A Dynamical 


vol. 9, no. 


Investigation of Oscillating Sp ~ 
Droplets, by T. V. Davies and D. 
Haydon, Proc. Royal Soc., Series A,, v: aL 
243, no. 1235, Feb. 11, 1958, pp. 492-499 

Experimental Verification of Nozzle 
Admittance Theory in Simulated Rocket 
Chamber, by Sotirios Lambiris, Project 
Squid, Tech. Rep. PR-74-R (ASTIA AD 
147309), Nov. 1957, 81 pp. (available only 
on microcards). 

Experimental and Theoretical Investi- 
gations of the Flow of Air Through Two 
Stage Compressors, by J. H. Horlock, ‘“t. 
Brit. Aeron. Res. Council, Rep. & Mem. 
3031 (formerly A.R.C. Tech. Rep. 17517), 
36 pp. 

Hydromagnetic Effects on Stagnation 
Point Heat Transfer, by Joseph L. Neu- 
ringer and William Mcellroy, Reputb/ic 
Aviation Corp., Nov. 1957, 8 pp., 1 tab., 4 
figs. 

Theoretical Study of Magnetohydrocy- 
namic Magneto-acoustic, and Magneto- 
elastic Phenomena, by Alfredo Banos 
Jr.. AFOSR-T R-57-35 (ASTIA AD }2- 
6561), May 1957, 30 pp. 

Comparison of the Take-off Noise 
Characteristics of the Caravelle Jet Air- 
liner and of Conventional Propeller Driven 
Airliner, by Laymon N. Miller and Leo L 
Beranek, Bolt Beranek and Newman, Inv., 
May 1957, 44 pp. 

The Mixing of Unbounded Coaxial 
Compressible Streams, by D. Roger 
Willis and Irvin Glassman, Project Squid, 
Tech. Rep. PR-70-P (ASTIA AD 146039 
Nov. 1957, 34 pp. (available only on mi- 
crocards). 

Pressure Temperature Relation for Con- 
stant Area Compressibility of a Gas Con- 
sidering Heat Transfer and Friction with 
Constant Wall Temperature, by Franklin 
P. Durham, Atomic Energy Comm., LA 
2129, Jan. 1957, 81 pp. 

Prediction of the Pressure Loss and 
Density Factors for Two Phase Annular 
Flow with or without Heat Generation, by 
J ». Westmoreland, Atomic Energy 
Comm., KAPL 1792, Feb. 1957, 35 pp. 

A Method for the Calibration of Flexible 
Plate Supersonic Wind Tunnels and Cali- 
bration Results for the 12 Inch Wind 
Tunnel at the Jet Propulsion Laboratory, 
by G. G. Goranson, R. M. Barnett, R. FE. 
Covey and H. N. Riise, Calif. Inst. of 
Techn., Jet Propulsion Lab., Rep. 20-110, 
June 1957, 78 pp. 

Relativistic Magnetohydrodynamics, by 
Edward G. Harris, Physical Rev., vol. 108, 
Dec. 15, 1957, pp. 1357-1360. 

Hydrogen Flammability Data and Ap- 
plication to PWR Loss of Coolant Acci- 
dent, by Z. M. Shapiro and T. R. Mof- 
fette, Atomic Energy Comm., WAPD- 
SC-545, Sept. 1957, 23 pp. 

Formulae and Approximations for Aero- 
dynamic ~— Rates in High Speed 
Flight, by J. Monaghan, Gt. Brit., 
Aeron. Res. eae Current Paper 260 
(formerly A.R.C. Tech. Rep. 18567; 
Roy. Aircr. Estab., TN Aero. 2407), 1957, 
31 pp., 19 figs. 

Analysis of Incompressible, Nonviscous 
Blade-to-blade Flow in Rotating Blade 
Rows, by J. J. Kramer, Trans., ASME, 
vol. 80, no. 2, Feb. 1958, pp. 263-275. 

Investigation of Burnout Heat Flux in 
Rectangular Channels at 2000 Psia, hy 
H.S. Jacket, J. D. Roarty and J. E. Zerbe, 
Trans., ASME, vol. 80, no. 2, Feb. 1958, 
pp. 391-401. 

Turbulent Heat Transfer in the Thermal 
Entrance Region of a Pipe with Uniform 
Heat Flux, by E. M. Sparrow, T. M. Hall- 
man and R. Siegel, Appld. Sci. Res., 
Section A, vol. 7, no. 1, 1957, pp. 37-52. 
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For over twenty years—from 
the early space labs of Annapolis to the launching 
pads at Cape Canaveral — setting the standards for 
the control of high pressure fluids...the regulators 
ijwhich served as the critical control on World War II 
flame throwers and torpedoes...and made missile 
sdevelopment work possible... Grove High Pressure 


Regulators... frequently imitated but never equalled. 


ROVE VALVE and REGULATOR COMPANY @: 


29 Hollis St., Oakland 8, California - 2559 W. Olympic Bivd., Los Angeles 6, California 
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_ Analysis of Heat-driven Oscillations of 
Gas Flows, Part III—Characteristic Equa- 
_ tion for Flame-driven Oscillation of the 
_ Organ-pipe Type, by H. J. Merk, Appld. 
_ Sct. Res., Section A, vol. 7, no. 2-3, 1958, 
pp. 175-191. 
__ Analysis of Heat-driven Oscillations of 
Gas Flows, Part IV—Discussion of the 
-‘Theoretical Results Concerning Flame- 
_ driven Oscillations, by H. J. Merk, Appld. 
Sci. Res., Section A, vol. 7, no. 2-3, 1958, 
pp. 192-204. 
Statistical Mechanics of Transport Proc- 
esses, XI, Equations of Transport in 
_ Multicomponent Systems, by Richard J. 
_ Bearman and John G. Kirkwood, J. 
_ Chem. Phys., vol. 28, no. 1, Jan. 1958, pp. 
136-145. 
_ Theory of Stagnation Point Heat Trans- 
fer in Dissociated Air, by J. A. Fay and 
F. R. Riddell, J. Aeron. Sci., vol. 25, no. 2, 
Feb. 1958, pp. 73-79. 
_ Stagnation Point Heat-transfer Meas- 
urements in Dissociated Air, by P. H. 
; Rose and W. I. Stark, J. Aeron. Sci., vol. 
a 25, no. 2, Feb. 1958, pp. 86-92. 

An Approximation of the Boundary of a 
Supersonic Axisymmetric Jet Exhausting 
into a Supersonic Stream, by E. L. Love, 

ie J. Aeron. Sci., vol. 25, no. 2, Feb. 1958, p. 


Transfer Through Liquid Metals, 
by G. Braudeau, Houille Blanche, vol. 12, 
oa 1957, pp. 81-86 (in French). 

Practical Calculation of Compressible 
Flow in Axial Turbines, by G. Corde:, 
Technik, vol. 12, 1957, pp. 279-284 (in 
i. Engineering Data for Diphenyl Cooled 
4 Nuclear Reactors, by L. W. Froom and 
_ Kermit Anderson, N uclear Sci. and Engng., 
2, 1957, pp. 160-168. 


KEY PROPULSION 
ASSIGNMENTS 
4 To men already seasoned in advanced 


power research, Chance Vought offers 
opportunity to join small-group R & D 
effort in nuclear and space propulsion, 
a and in advanced fluid mechanics. 


Senior Propulsion Specialist. Advanced 
Engineering or Physics Degree with a 
os minimum of 6 years opplicable experi- 


ence. To conduct research and develop- 
ment in the fields of advanced heat 
a transfer and space environmental sys- 
Ae tems; pursue technical studies of specific 
cooling methods for hypersonic vehicles; 
interpret advanced hypersonic cooling 
| progress and advise Chief of Propulsion. 


+ Engineering Specialist. | Advanced A.E. 
or M.E. Degree preferred with a mini- 
mum of 5 years experience. To conduct 
7 applied research in advanced fluid 
mechanics in high supersonic region; per- 
form design analysis of high Mach num- 
7 q ber aircraft and missile inlets and exits, 
-: and technical studies of fluid dynamics 
om and control devices. 


Chief of Propulsion 


Dept. JP-7 


AIRCRAFT 
7 Po DALLAS, TEXAS 


Further Investigations of the Laminar- 
turbulent Transition in a Free Jet (An- 
nular Nozzle), by O. Wehrmann and H. 
Fabian, Hermann Fottinger-Inst. fur Stro- 
munstech., Berlin Tech. U. Rep. (AFOSR 
TR 57-31) (AD 126494), Dec. 1956, 1 vol. 

On the Theory of Binary Fluid Mixtures, 
by F. J. Pearson and G. 8. Rushbrooke, 
Proc. Royal Soc. (Edinburgh), Sect. A, 
1955-57, pp. 305-317. 


Combustion, Fuels 
and Propellants 


Oriole ...A Really Low-Cost Research 
Missile, Missiles and Rockets, vol. 2, Dec. 
1957, pp. 86-87. 

Mechanisms for the Formation of Ions 
in Flames, by H. F. Calcote, Project 
Squid, Tech. Rep. EXP-6-P (ASTIA AD 
145989), Nov. 1957, 44 pp. (available only 
on microcards). 

NACA Conference on Combustion; A 
Compilation of the Papers Presented by 
NACA Staff Members, by V ACA, March 
26, 1947, 77 pp. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 123, p. 18, 1/7/58.) 

Studies in Bomb Calorimetry VIII— 
Design of Calorimeter Systems, by R. A. 
Mott, Fuel, vol. 37, Jan. 1958, pp. 3-18. 

Stability of Propane-air Flames in Vor- 
tex Flow, by A. E. Potter Jr., E. L. Wong 
and A. L. Berlad, NACA TN 4210, Feb. 
1958, 27 pp. 

New Methods for the Measurement of 
Relative Ignitability and Ignition Effi- 
ciency, by Cecilio R. Grande, Picatinny 
Arsenal, Samuel Feltman Ammunition 
Labs., Tech. Rep. 2469, Feb. 1958, 1 vol. 


Turbulent Burning Velocities of Natural 
Gas-air Flames with Pipe Flow Turbu- 
lence, by J. Kenneth Richmond, Joseph 
M. Singer, E. B. Cook, James R. Oxen- 
dine, Joseph Grumer, and David S. Bur- 
gess, Bur. Mines, Div. Explosives Techn., 
Tech. Rep. 1205 (AFOSR-TN-56-238, 
ASTIA AD 88358), Apr. 1956, 16 pp. 

Combustion in Laminar Mixing Regions 
and Boundary Layers, by Donald Allen 
Dooley, Calif. Inst. of Techn., Jet Propul- 
ston Lab., Mem. 20-151, June 1956, 244 pp. 


Shock Waves and Flame Interaction, by 
George Rudinger, Project Squid, Tech. 
Rep. CAL-74-P (ASTIA AD 147699), 
Dec. 1957, 42 pp. (available only on mi- 
crocards). 

The Ignition of Combustible Mixtures 
by Hot Gases, by H. G. Wolfhard, Proj- 
ect Squid, Tech. Rep. BUM-24-P (ASTIA 
AD 147702), Dec. 1957, 27 pp. (available 
only on microcards). 


The Study of Flow and Reaction Rates 
in Turbulent Flames, by R. P. Barbor, 
J. D. Larkin and C. W. Shipman, Project 
Squid, Tech. Rep. DEL-8-P (ASTIA AD 
147703), Dec. 1957, 38 pp. (available only 
on microcards). 


Theoretical Rocket Performance of 
4 Fuel with Several Fluorine-oxygen Mix- 
tures Assuming Equilibrium Composition, 
by Sanford Gordon, NACA Res. Mem. 
E57K22, Feb. 1958, 69 pp. 


Infrared Analysis of Mixtures of Nitro- 
glycerine and Ethylene Glycol Dinitrate, 
by Woodruff Huff, Michael Halik and 
Frank Pristera, Picatinny Arsenal, Sam- 
uel Feltman Ammunition Labs., Tech. Rep. 
2472, Dec. 1957, 9 pp. 

An Approximate Treatment of Hydra- 
zine Decomposition in a Laminar Non- 
isothermal Flow, by Mitchell Gilbert, 
Calif. Inst. Techn., Jet Propulsion Lab., 
Progr. Rep. 20-336, July 1957, 25 pp. 


Combustion of Elemental Boron. Sum- 
mary Report for June—October 1957, 
Experiment Inc., TM-1009, 10 pp. 

Phase Behavior in the Hydrogen Am- 
monia System, by H. H. Reamer and B. 
H. Sage, Calif. Inst. Techn., Jet Propulsion 
Lab., Progr. Rep. 20-237, Dec. 1957, 9 pp. 

Experimental Thermal Conductivities of 
the N.O, =2NO, System, by Kenneth P. 
Coffin and Cleveland O’Neal Jr., NACA 
TN 4209, Feb. 1958, 22 pp. 

Final Report Summarizing Research in 
Rate of High Speed Reactions, by Walter 
H. Wurster, Cornell Aeron. Lab., Inv., 
Rep. AD-959-A-2 (AFOSR-T R-58-8; AS- 
TIA AD 148054), Jan. 1958, 11 p)., 
8 figs. 

Unified Theory of Internal Ballistics, |v 
J. N. Kaper, Quart. J. Mech. and Apyl. 
Math., vol. 11, Part I, Feb. 1958, pp. 98 
111. 

Combustion Efficiency Performance of a 
Mil-F-5624 Type Fuel and Monomethy|- 
naphthalene in a Single Vaporizing-type 
Combustor, by Anthony W. Jones and W'I- 
liam P. Cook, NACA Res. Mem. E51K3), 
Feb. 1952, 26 pp., diagrs., photos., 2 tabs. 
(Declassified from Confidential by ai- 
thority of NACA Res. Abstr. 123, p. 8, 
1/7/58.) 

Turbojet Combustor Efficiency with 
Ceramic-coated Liners and with Mechani- 
cal Control of Fuel Wash on Walls, liy 
Helmut F. Butze and Edmund R. — 
NACA Res. Mem. E52125, Nov. 1952, 
pp. (Declassified from Confidential 
authority of Res. Abstr. 124, 
p. 41, 2/10/58.) 

Effectiveness of a Turbojet Tubular 
Combustor in Screening the Turbine from 
Foreign Objects, by Patrick T. Chiarito, 
NACA Res. Mem. E55E16, July 1955, 20 
pp. (Declassified from Confidential by 
authority of NACA Res. Abstr. 124, 
p. 44, 2/10/58.) 

Spontaneous Ignition Data of Hydro- 
carbons and Aviation Fluids, by E. M. 
Goodger, Cranfield, Coll. of Aeron., Nolte 
68, Sept. 1957, 6 pp., 10 figs. 

Report of Research and Technologic 
Work on Explosives, Explosions, and 
Flames; Fiscal Years 1953 and 1954, by 
Ruth F. Brinkley and Robert W. van 
Dolah, Bur. Mines, Info. Circ. 7804, Nov. 
1957, 127 pp. 

Preliminary Investigation of a Chemical 
Starting Technique for the Acid Gasoline 
Rocket Propellant System, by Glen Hen- 
nings and Gerald Morrell, NACA Res. 
Mem. E52K21, Jan. 1953, 23 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 124, p. 41, 2/10/58. ) 

Methods for Calculating Thrust Aug- 
mentation and Liquid Consumption for 
Various Turbojet-afterburner Fuels, by 
James F. Morris, NACA Res. Mem. k- 
56A23, Oct. 1956, 73 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 124, p. 44, 2/10/58.) 

_A Study of Sampling of Flame Gases, hy 

’. Halpern and F. W. Ruegg, J. Res., Nat. 
Bur Standards, vol. 60, no. 1, Jan. 1958, 
pp. 29-37. 

How to Estimate Solid Propellant 
Rocket Performance, by Saxe Dobrin, 
Aviation Age, vol. 28, no. 8, Feb. 1958, 
pp. 70-73. 


Materials of Construction 


Electrical, Thermoelectric, Hardness 
and Corrosion Properties of Vanadium- 
base Alloys, by H. J. Cleary, Nuclear 
Metals, Inc. (U.S. Atomic Energy Comm.), 
NMI-1161, Sept. 1957, 1 vol. 

Thermal and Electrical Conductivities of 
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Tungsten. Experimental Results on Sin- 
gle Crystals Compared with Theory, by J. 
DeNobel, Physica, vol. 23, March 1957, 
pp. 261-269 (in English). 

Thermal and Electrical Resistivity of 
Some Tungsten Single Crystals at Low 
Temperatures and in Strong Magnetic 
Fields, by J. DeNobel, Physica, vol. 23, 
April 1957, pp. 349-358 (in English). 

The Mechanical and _ Engineering 
Properties of Commercially Available Tita- 
nium Alloys, by H. V. Kinsey, NATO, 
Advisory Group for Aeron. Res. and Dev., 
Rep. 100, April 1957, 15 pp. 

How to Design Radome Structures for 
High Speed Aircraft and Missiles, by R. 
M. Kubow and N. J. Linardos, Aviation 
{ve, vol. 28, no. 8, Feb. 1958, pp. 74-79. 


Instrumentation, Tele- 
metering, Data Recording 


Central Automatic Data Processing 
System, by the Staff of the Lewis Labora- 
tory. Chapter I: General Description, 
by Bert A. Coss. Chapter II: Central 
Recording System, by Robert L. Miller. 
Chapter III: Automatie Voltage Digi- 
tizers, by Leonard Jaffe and Richard 
L. Smith. ChapterIV: Frequency Data, 
by John Ryskamp. Chapter V: Digital 
Automatic Multiple Pressure Recorder, 
by Leonard Jaffe, Arthur J. Gedeon and 
tichard N. Bell. Chapter VI: Playback 
and Control Room Equipment, by John 
Rvskamp, NACA TN 4212, April 1958, 
96 pp. 

New Breed of Microwave Tubes 
Spurs Radar, ECM, Scatter Designs, by 
James Holahan, Aviation Age, vol. 29, 
April 1958, pp.22-23, 157-164. 

Closed Loop Instrument Systems Boost 
Performance, by V. A. Orlando, R. G. 
Jewell and D. W. Dibble, Aviation Age, 
vol.29, April 1958, pp.44—47. 

Commutating Switches Critical for 
High Performance Telemetry, by George 
P. Bentley and Sumner Ackerman, 
Aviation Age, vol. 29, April 1958, pp. 
70-72, 75-78. 

Mechanoreceptors, Gravireceptors, by 
Hubertus O. Strughold, J. Astronautics, 
vol.4, no.4, Winter 1957, pp.61-63. 

Cosmic Ray Instrumentation for the 
IGY Program, by Robert C. Haymes, J. 
Astronautics, vol. 4, no. 4, Winter 1957, 
pp. 64-68. 

Reduced-size Satellite Transmitters 
Coming, by Raymond M. Nolan, Mis- 
siles and Rockets, vol. 3, April 1958, pp. 
152-153. 

An Engineering Approach to Cost 
Reduction as Applied to the Engine 
Analyzer, by J. Frank Morgan and 
William Van Rosenbergh, Sperry Engng. 
Rev., vol. 11, Mar. 1958, pp. 22-26. 

Development of the Calorimeter Heat 
Transfer Gage for Use in Shock Tubes, 
by Peter H. Rose, AVCO Mfg. Corp., 
AVCO Res. Lab., Res. Rep. 17, Feb. 
1958, 49 pp. 

The Satellite Tracking Camera, by 
Jospeh Nunn, Inst. Aeron. Sci., Preprint 
794, Jan. 1958, 16 pp., 2 figs. 

The Army Navy Instrumentation Pro- 
gram, by George W. Hoover, Inst. Aeron. 
Sci., Preprint 815, Jan. 1958, 9 pp 

The Air Force Instrument Development 
Program, by Marvin C. Demler, /nst. 
Aeron. Sci., Preprint 795, Jan. 1958, 15 
pp., 43 figs. 

Temperature and Pressure Corrections 
to be Applied to the Shielded Hot Wire 
Anemometer at Speeds for Which Natural 
Convective Cooling Is Negligible, by C. F. 
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Man’s major problem in his use of energy has 
always been materials. As we advance into more 
complex technologies the need for a more practical ~~ 


method of working with high temperature and high 

stress materials has become compelling. Now, Plasma © 7 
Spray, a radically new device capable of vaporizing © 
and propelling any element or combination of ele- | 
ments, offers the materials engineer new horizons in — 

the design of coating compositions and properties. | 


This new tool extends the range and capabilities of J 


materials processing into refractory metals, 
carbides, borides, silicides and nitrides. With Plasma © = 
Spray practically any material can be spray coated © 7 


on another. 

The Plasma Spray system consists of a standard 
Plasmatron head fitted with a concentric spray 
attachment. With its inert jet of intense heat and con- 
trolled temperatures to 15,000°F, the instrument can 
vaporize, spray and coat even the most reactive and 
refractory of materials. Plasma Spray utilizes the 
assemblage of electrically neutral, partially ionized 
gas with an electric arc created by the Plasmatron. 
Obtaining heat from a constricted and sustained dis- 
charge, plasma temperatures are three to five times 
that of an oxyacetylene flame, and the inert jet offers 
a completely uncontaminated spraying environment 
which does not oxidize either the base or the 
sprayed material. 

Plasma Spray places in the engineer’s hands an infi- 
nite number of combinations of physical properties and 
coating compositions for more practical solutions of 
high temperature, high stress materials requirements. _ 


*ZrB2 coating illustrated abave is from 90-95% of theoretical density 
being sprayed at 15,000°F. Melting point of ZrB2 is 5500°F. Micro- 
hardness of coating is 2000 Kilograms per mm. under 50 gram load. 


Write for price and delivery. 


Plasmadyne 


Corporation 


3839 South Main Street, Santa Ana, California 
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CONTROL ENGINEERS 


(Electric Servo « Valves) 


Move into Large Rocket Engineering eat 


and put yourself way ahead in your field 


Help us to automate millions of horsepower designed into a 
jet-size package—the High-Thrust Rocket Engine. Here are the 
fields : 

The Electrical System includes Ground Support and Check- 
out Equipment which must be operable by military personnel. 
Aboard the missile, engine controls must be carefully isolated 
from other missile systems. Miniaturization is striven for, but 
never at the expense of reliability in extremes of temperature, 
vibration and acceleration. You'll cover all aspects of circuitry, 
deal with every branch of weapons systems. 

Servo-mechanisms offer a broad spectrum—electronic, pneu- 
matic, mechanical, hydraulic. Your analytical ability will be at 
a premium here, to evaluate methods of Mixture Control, Thrust 
Control, and Pressure Control which must compensate for vari- 
ables like changing mass, drag lapse-rate, altered combustion 
efficiency, heat, cold, G, vibration etc. You'll be free of routine 
details, able to apply your training and experience toward a high 
level of professional growth. 

Valves run to 6” diam. and up, with very high pressures and 
flow rates, extremely rapid action, temperatures down to —300°F. 

This is where the real advanced work in controls is being done. 
Join the trailbreakers. Write, giving your background: Mr. 
A. L. Jamieson, Rocketdyne Engineering Personnel, 6633 
Canoga Avenue, Canoga Park, California. — 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIAT 


BUILDERS OF POWER FOR OUTER SPACE 


eee 


Cowdrey, Brit. J. Appld. Phys., vol. 9, 
Mar. 1958, pp. 112-116. 

Strain-gage Principles, by M. A. 
Legette, Instrs. & Automation, vol. 31, 
Mar. 1958, pp. 447-449. 

Strain-gage Instrumentation, by C. \. 
Hathaway, Instrs. & Automation, vol. 31, 
Mar. 1958, pp. 450-454. 

Optical Strain-gage Standard, /nsirs. «& 
Automation, vol. 31, Mar. 1958, p. 455. 

Sensitivity Ratings for Strain Gages, |\y 
Peter. K. Stein, Jnsirs. & Automation, vol. 
31, Mar. 1958, pp. 456-561. 

High-temperature Strain Gages, 
Louis Herezeg and Paul Beckman, /nsts. 


& Automation, vol. 31, Mar. 1958, pp. 
460-461. 
Thermal Compensator for Bourdon 


Gauges, by Otto W. Heise, Jnsirs. + 
Automation, vol. 31, Mar. 1958, p. 473. 

Methods of Testing Thermocouples ar.d 
Thermocouple Materials, by William *. 
Roeser and S. T. Lonberger, Nat. Bur. 
Standards, Circ. 590, Feb. 1958, 21 pp. 

Time Constants and Frequency Re- 
sponse of Coated Hot Wires Used as 
Turbulence Sensing Elements, by Avis 
Borden, David Taylor Model Bocin. Rep. 
952, June 1957, 48 pp. 

Electronics: What’s Coming After the 
Missile Age? by W. R. G. Baker, Prov., 
Inst. Radio Engrs., vol. 46, Mar. 1958, pp. 
534-537. 

Thermoelectric Effects, by Frank !:. 
Jaumot Jr., Proc., Inst. Radio Engrs., vol. 
46, Mar. 1958, pp. 538-554. 

Unusual Propagation at 40 MC from the D 
USSR Satellite, by H. W. Wells, Prov., 
Inst. Radio Engrs., vol. 46, Mar. 1958, p. 
610. 

A Note on Some Signal Characteristics - 
of Sputnik I, by J. D. Kraus and J. 5. 
Albus, Proc., Inst. Radio Engrs., vol. 46, d 
Mar. 1958, pp. 610-611. Es 

Intermetallic Semiconductors, by Henry 
T. Minden, Sylvania Technologist, vol. 11, 
Jan. 1958, pp. 6-9. m 

Apparatus for Producing and Measuring 
High-energy Electrical Discharges, hy 
W. E. Richeson, Rev. Sci. Instr., vol. 29, 9M: 
Feb. 1958, pp. 99-104. vo! 

Novel Circuit for Delivering Known 
Spark Energies, by T. A. Erikson, Rev. py 
Sci. Instr., vol. 29, Feb. 1958, pp. 173-174... JAp 

Piezoelectric Detector for Low-pressure ( 
Shock Waves, by Herbert T. Knight, Rev. §0rc 
Sci. Instr., vol. 29, Feb. 1958, pp. 174-175. §Yl4 

Indium Resistance Thermometer, by f 
G. K. White, 8. B. Woods and F. Anglin, Bez 
Rev. Sci. Instr., vol. 29, Feb. 1958, pp. Wshi 
181-182. 

Backlighting the Shock Arounda Moving RP 
Explosive Charge, by D. J. Hinz and J. §gy; 
Wenig, Aberdeen Proving Ground, Ballistic 
Res. Labs., Mem. Rep. 1110, Oct. 1957, 16 
pp. 

Aeronautical Studies in the Aerobal- 
listic Range, by G. V. Bull, Canadian 
Armament Res. and Dev. Estab., Rep. 
302/57, July 1957, 25 pp. 

Comparison of Ballistic and ‘‘Breguet”’ 
Ranges, by D. Dembrow, Johns Hopkins 
Univ., Appld. Phys. Lab., CF-222, March 
1958, 8 pp., 5 figs. 

Redesign Adapts Accelerometers to 
High Speed Use, by Daniel T. Gundersen, 
Aviation Age, vol. 29, May 1958, pp. 132 
141. 


Cyclops Cores Simplify Earth-satellite 
Circuits, by W. Matthews, R. W. Rochelle, 
C. B. House, R. L. Van Allen, D. H. 
Schaefer and J. C. Schaffert, Electronics, 
vol. 31, no. 9, Feb. 28, 1958, pp. 56-63. 


Approximation Methods for the Caicula- 
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tion of Ultra-high-frequency Radio Wave 
Field Strength, Considering the Effects of 
Local Terrain, by A. N. Kalinin, Radio- 
tekhnika, vol. 12, no. 4, April 1957, 
pp. 13-23 (in Russian). 

A Method of Determining the Thermal 
Conductivity, Electrical Resistivity and 
Total and Spectral Emissivity of a Con- 
ducting Material at High Temperature 
from Electrical and Thermo-optical Meas- 
urements on a Heated Tube in a Vacuum, 
by Walter B. Powell, Calif. Inst. of Techn., 
Jet Propulsion Lab., Progr. Rep. 20-324, 
July 1957, 17 pp. 

Experimental Investigation of the Sur- 
face Impact Pressure Probe Method of 
Measuring Local Skin Friction at Super- 
sonic Speeds, by Charles J. Stalmach Jr., 
Univ. of Texas, Defense Res. Lab., CF 
2675 (DRL-410), Jan. 1958, 53 pp., 7 tabs., 
36 figs. 

A Contact Extensometer for the Record- 
ing of Flight Load Spectra, by O. Svenson, 
NATO, Advisory Group for Aeron. Res. 
an! Dev., Rep. 121, May 1957, 10 pp. 

A Study of Error Effects in Measuring 
Cyclic-temperature MHeat-transfer Co- 
efficients, by David A. Dingee and Joel W. 
Chastain, Battelle Mem. Inst. (U.S. Atomic 
Energy Comm.), BMI-1167, Feb. 1957, 


1 vol. 


Guidance Systems 
and Components 
1958 Missile Guidance Roundup, by 


Heyward E. Canney Jr., Missiles and 
Rockets, vol. 3, no. 2, Feb. 1958, 
pp. 77-83. 


How Industry Solved the Air-bearing 
Gyro Stabilization Problem, by John P. 
Jagy, Missiles and Rockets, vol. 3, no. 2, 
Feb. 1958, pp. 84-86. 

Workhorse of Inertial Guidance, by 
Earl Finkle, Missiles and Rockets, vol. 3, 
no. 2, Feb. 1958, pp. 91-93. 

Tracking Long-range Missiles, by Ray- 
mond M. Nolan, Missiles and Rockets, 
vol. 3, no. 2, Feb. 1958, pp. 125, 127. 


Target Simulator Tests Beam-rider 
Missiles, by C. E. Hendrix, Electronics, 
vol. 31, no. 5, Jan. 31, 1958, p. 32. 

Star Sensors for Automatic Navigation, 
by R. B. Horsfall, Aviation Age, vol. 29, 
April 1958, pp. 150-154. 

Guiding the Terrier, by R. M. Nolan, 
Ordnance, vol. 42, March-April 1958, pp. 
14-915. 

An Analysis of the Optimization of a 
Beam Rider Missile System, by Marvin 
Shinbrot and Grace C. Carpenter, NACA 
I'N 4145, Mar. 1958, 34 pp. 

Proceedings, Guided Missiles Seminar: 


Guidance and Control, Sept. 2, 1956, 
NATO Aav. Group for Aeron. Res. and 
Dev., AGARDograph 21, Sept. 1956, 291 
pp. Contents: Weapons System Phil- 
tsophy, by George H. Clement, pp. 1-10; 
ORO Weapon System Philosophy, by 
Harlan C. Meal, pp. 11-28; New Prin- 
tiples in the Design of Superior Communi- 
tatiuns, Navigation and Missile Guidance 
Systems, by W. P. Lear Sr., pp. 29-38; 
idance Techniques, by Walter L. 
Webster Jr., pp. 29-56; Considerations in 
he Choice of a Missile Guidance and 
Control System, by Robert W. Mayer, pp. 
7 66; Inertial Navigation, by Norman F. 
Parker and Charles P. Greening, pp. 67- 
Aiding the Inertial Navigation Sys- 
m, by William F. Ballhaus and Frederick 
tevens Jr., pp. 87-100; Linear Homing 
Navigation, by Robert K. Roney, pp. 101- 
ll4; Pitfalls in Missile Control, by 
Robert L. Johnson, pp. 115-138; The 
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Phone or Write 


MISSILE 


HARDWARE 
| By 


20 Mechanic Street 


the “Know- How!” 


We have developed new tech- 
niques, new methods, new 
processes that effect production 
economy so necessary to a success- 
ful missile program. 


Here at Newbrook you will find 
men with experience gained 
from doing ...a modern plant 
with up-to-date equipment... 
precision inspection to meet your 
most exacting quality control 
requirements. 


And most important, Newbrook 
specialization results in strict 
reliability ! Let us help you with 
your Missile Hardware problems. 


Phone 45 


SILVER CREEK, NEW YORK 


we have ad 
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MARQUARDT 


Professional Personnel 
Requisition 


AERODYNAMICIST 


Internal Flow Systems 


To formulate and test inlet systems and 
exit nozzles for hypersonic air breath- 
ing propulsion devices. Create new 
designs for advanced and unique en- 
gine configurations. Supervise aero- 
dynamic design of hypersonic free jet 
and wind tunnel research facilities. 


Five years’ experience in internal flow 
and fluid dynamics analysis. Re- 
quires Mechanical or Aeronautical 
Engineering degree with appropriate 
advanced study. 


Small work group offers unusual oppor- 
tunity for professional recognition 
and advancement. 


Contact: Floyd E. Hargiss, Manager 
Professional Personnel Department 

Marquardt Aircraft Co. 

16556 Saticoy Street 

Van Nuys, California 


Marquardt c 


COMPLETE 


TESTING 


FACILITIES 


* Qualification Tests 
* Evaluation Tests 

* Performance Tests 
* Environment Tests 


Manufacturers of Metal Boss 
Seals to Military Specifications 


AIRCRAFT 
EQUIPMENT 
TESTING CO. 


1812 Fleet St., Baltimore 31, Md. 


ORleans 5-8337 ORleans 5-2222 


Control System Design, by F. E. Perry, 
pp. 139-152; Geometrical Stabilization 
Based on Sevodriven Gimbals and Inte- 
grating Gyro Units, by on 8. Draper 
and Roger B. Woodbury, 153-182; 
Sampled Data Systems, by Jo 3 Ragaz- 
zini, pp. 183-228; Digital Techniques in 
Missile Guidance Systems, by Sidney 
Darlington, pp. 229-248; The Use of 
Digital Computer Techniques in Missile 
Design and Control, by D. H. Gridley, pp. 
249-264; The Application of Noise and 
Filter Theories to Guidance Problems, by 
R. J. Parks and Robert M. Stewart, pp. 
265-284; Recent Developments in Fixed 
and Adaptive Filtering, by A. G. Carlton 
and J. W. Folin Jr., pp. 285-300; Practical 
Problems Encountered in Missile Guid- 
ance and Control Design, by R. E. Whif- 
fen, pp. 301-314; Application of Methods 
of Science to the Problem of Reliability, 
by C. Raymond Knight, pp. 315-330; 
Retebiiey of Guided Missiles, by Edwin 

A. Speakman, pp. 331-340; Laboratory 
and Flight Evaluation of Airborne Guid- 
ance Systems, by W. H. Clohessy, pp. 
341-346; Trends in Field Testing of 
Guided Missiles, by Ernest A. Steinhoff, 
pp. 347-358; Low Signal Level Missile 
Instrumentation, by L. G. DeBey, pp. 
359-374; On the Way to Automated 
Processing of Flight Measurement, by 
W. E. Klemperer, pp. 375-386: Paper on 
the Guidance and Control of Missiles, by 
Stephen Waldron, pp. 387-391. 

Attenuation of Radio Waves as a Result 
of Boundary Layer Heating in Missiles 
Traveling at Supersonic Speeds, by Kil- 
burn MacMurraugh, NAVORD Rep. 
5816 (NOTS 1896), Jan. 1958, 18 pp. 

Guidance and Control, by Peter A. 

Castruccio, Aviation Age, vol. 28, Mar. 
1958, pp. 64-68. 
Traveling Wave Slot Makes Novel X-band 
Beacon Antenna, by L. K. DeSize and 
L. J. Kuskowski, Aviation Age, vol. 28, 
Mar. 1958, pp. 188-192. 

Atmospheric Angels Mimic Radar 
Echoes, by Vernon G. Plank, Electronics, 
vol. 31, Mar. 14, 1958, pp. 140-144. 

Simple Plotter Analyzes Radar Noise 
Rapidly, by Daniel J. Zoll, Electronics, vol. 
31, Mar. 14, 1958, pp. 162-164. 

Inertial Guidance, by Walter Wrigley, 
Roger B. Woodbury and John Hovorka, 
Inst. Aeron. Sci., SMF Fund Paper FF-16, 
Jan. 1957, 69 pp. 


Servomechanisms 
and Controls 


The Improvement of Stabilization of 
Control Systems with Bounded Speed of 
the Servomotor by Means of . Memory 
Device, by V. A. Kotelnikov, Automatika 
i Telemekanika, vol. 18, no. 4, 1957, 
pp. 289-295 (in Russian). 

Signal Stabilization of a Control System, 
by R. Oldenburger, Automatika 7 Tele- 
mekanika, vol. 18, no. 5, 1957, pp. 392-397 
(in Russian). 

The Syathesis of Linear Servosystems 
Using the Criterion of Minimum of Prac- 
tically Critical Reproduction Error, by 
K. I. Kurakin, Automatika i Telemekanika, 
vol. 18, no. 5, 1957, pp. 402-408 (in Rus- 
sian). 

An Analog Study of a High-speed ~ 
cording by J. 
Schwartzenberg, Trans., ASME, vol. a, 
no. 2, Feb. 1958, pp. 490-496. 

A Method for Designing a Sampling 
Servo System, by Saburo Kumagai and 
Akira Nishigori, Osaka Univ., Technology 


Reps., vol. 7, March 1957, pp. 19-30. 
Intelligence Sources for Guidance, by 


Arthur 8. Locke, Jnst. 
print 795, Jan. 1958, 15 pp., 43 figs. 
The Development of the Talos Land 


Aeron. Sci., Pre- 


Based System, by D. B. Holmes, Jnst. 
Aeron. Sci., Preprint 804, Jan. 1958, 9 pp., 
9 figs. 

Report on the Endurance of Irreversible 
Servo Controls, by L. Vandenberghe, 
NATO, Adv. Group for Aeron. Res. und 
Dev., Rep. 128, May 1957, 6 pp. (in 
French). 

Improvements of the Transient Re- 
sponse of Contractor Servomechanismns, 
by Yoshikazu Sawaragi and Hajime 
Akashi, Proc., Japan, 6th Natl. Congr. for 
Appld. Mech., 1956, pub. 1957, pp. 441 
464. 


On the Self-excited Oscillations of a 
Non-linear Cascade Control System, by 
Yoshikazu Sawaragi, Yoshifumi Sunahira 
and Akira Someya, Proc., Japan., 6th 
Natl. Congr. for Appld. Mech., 1956, pub. 
1957, pp. 477-482. 


Flight Vehicle Design 
and Testing 
History of German Guided Missiles 


Development, by Th. Benecke and A. W. 
Quick, Brunswick, Germany, Appelh:ns, 


1957, 420 pp. (North Allantic _o 
Organization, Advisory Group for Aeron. 
Res. and Dev., AGARDograph no. 20.) 


Rockets, Missiles and Space Vehicles, 
Interavia, vol. 13, no. 2, Feb. 1958, 
pp. 150-151. 

NACA Hypersonic Rocket and High 
Temperature Jet Facilities, by Paul EP. 
Purser and Aleck C. Bond, NATO, Ad- 
visory Group for Aeron. Res. and Dev., Rep. 
140, July 1957, 13 pp. 

Red Missile Arsenal, Aviation Age, vol. 
29, April 1958, pp. 20-21, 94-96, 98-101. 

Some Variation Problems Connected 
with the Launching of Artificial Satellites 
of the Earth, by D. E. Okhotsimskii and 
T. M. Eneev, J. Brit. Interplanet. Soc., 
vol. 16, Jan.—Feb. 1958, pp. 263-294. 

Rockets, Missiles and Space Vehicles, 
Interavia, vol. 13, March 1958, pp. 243 
245. 

Recoverable Boosters are Studied to 
Cut Manned Space Flight Cost, by D. C. 
Romick, R. A. Belfiglio and F. B. Sand- 
gren, Missiles and Rockets, vol. 3, April 
1958, pp. 95-96, 98, 100. 

Role of Guided Missiles, by Dr. C. C. 
Furnace, Ordnance, vol. 42, March-April 
1958, pp. 817-819. 

Generalized Variational Approach to the 
Optimum Thrust Programming for the 
Vertical Flight of a Rocket, PartI: Neces- 
sary Conditions for the Extremum, by 
Angelo Miele, Zeitschrift Plugwissen- 
schaften, vol. 6, March 1958, pp. 69-77 
(in E nglish). 


Guided Weapons and Aircraft—Some 
Differences in Design and Development, 
by J. E. Serby, J. Roy. Aeron. Soc., vol. 
62, Mar. 1958, pp. 187-202. 

Application of the Electronic Digital 
Computer to the Solution of Flight Dy- 
namics Problems, by Werner Niemz, 
Zeitschrift fiir Flugwissenschaften, vol. 6 
Feb. 1958, pp. 47-52 (in German). 

Missile Market Guide and Directory, 
Missiles and Rockets, vol. 3, No. 5, Mid- 
April 1958, 462 pp. 


Effects of Size of Ballistic Missiles, by 
Oscar Scholze, Raketentechnik und Rawm- 


fahrtforschung, vol. 2, no. 1, 1958, pp. 2-8. 


Entering the ee, by Terence 
Nonweiler, Spaceflight, vol. 1, no. 7, April 
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1958, pp. 238-240. 

A Test and Launching Standard for the 
Vanguard Rockets, by Don Karshan, 
Spaceflight, vol. 1, no. 7, April 1958, pp. 
253-257. 

Motion of a Ballistic Missile Angularly 
Misaligned with the Flight Path upon 
Entering the Atmosphere and Its Effect 
upon Aerodynamic Heating, Aerodynamic 
Loads and Miss Distance, by H. J. Allen, 
NACA TN 4048, Oct. 1957, 66 pp. 

A Study of the Motion and Aerody- 
namic Heating of Missiles Entering the 
Earth’s Atmosphere at High Supersonic 
Speeds, by H. Julian Allen and A. J. 
Eggers Jr.. NACA TN 4047, Oct. 1957, 
61 pp. 

New Methods in Heat Flow Analysis 
fy th ary to Flight Structures, by 

A. Biot, J. Aeron. Sci., vol. 24, Dee. 
1037, pp. 857-872. 

Flight Measurements of Boundary- 
layer Temperature Profiles on a Body of 
Revolution (NACA RM-10) at Mach 
Numbers from 1.2 to 3.5, by A. G. Swan- 
son, J. J. Buglia and L. T. Chauvin, NACA 
TN 4061, July 1957, 40 pp. 

Theory and Apparatus for Meaurement 
of Emissivity for Radiative Cooling of 
Hypersonic Aircraft with Data for Inconel 
and Inconel X, by William J. O’Sullivan 
Jr., and William R. Wade, NACA TN 
4121, Oct. 1957, 48 pp. 

Interception Problems for Surface to 
Air Missiles, by James Hay Stevens, 
Acronautics, vol. 38, Mar. 1958, p. 27. 

Ballistic Missiles, by Terence R. F. 
Nonweiler, Aeronautics, vol. 28, Mar. 1958, 
pp. 28-31. 

Instruments of Interception (Missile 
Catalogue), Aeronautics, vol. 28, Mar. 
1958, pp. 38-40. 

A Note on Goddard’s Problem, by G. 
Leitmann, Astronautica Acta, vol. 3, no. 4, 
1957, pp. 237-240. 


Flight 
pace Flig 


A Note on Space Travel in a Gravita- 
tional Field, by E. Newton Rowland, M 
A., Ph.D., J. Brit. Interplan. Soc., vol. 16, 
Oct.—Dec. 1957, pp. 216-221. 

Effects of Acceleration on the Human 
Frame, by Wg. Cdr. F. Latham, J. Brit. 
Interplan. Soc., vol. 16, Oct.-Dee. 1957, 
pp. 222-224. 

Gravitational Torque on a Satellite 
Vehicle, by Robert E. Roberson, J. Frank- 
lin Inst., vol. 265, Jan. 1958, pp. 13-22. 


The True Nature of the Boiling of Body 
Fluids in Space, by Capt. Julian E. Ward, 
J. Aviation Medicine, vol. 27, Oct. 1956, 
pp. 429-439. 


Astrophysics, Aerophysics 


A Continuously Recording Automatic 
Auroral Radar, by McNamara, 
Canadian J. Physics, vol. 36, Jan. 1958, 
pp. 1-8. 

On the Theory of Propagation of Elec- 
ey age Waves Along a Curved Sur- 
face, James R. Wait, Canadian J. 
Physics, vol. 36, Jan. 1958, pp. 9-17. 


A Large-area Liquid Scintillation Coun- 
ter and Some Measurements on High- 
energy Cosmic-ray Particles, by C. H. 
Millar, E. P. Hineks and G. C. Hanna, 
Canadian J. Physics, vol. 36, Jan. 1958, 
pp. 54-72. 


Further Interactions of the Heavy Nu- 
clei of the Cosmic Radiation, by V. Y. 
Rajopadhye and C. J. Waddington, H. H. 
Wills Phys. Lab., Bristol (England), 
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Philosophical Mag., 8th Series, vol. 3, 


Jan. 1958, pp. 19-32. 


Nuclear Propulsion 


Stability of the Pinch, by Marshall 
Rosenbluth, Atomic Energy Comm., LA- 
2030, June 1956, 36 pp. 

Conceptual Design of an Advanced 
Engineering Test Reactor, by R. G. Mal- 
lon, J. Saldick and R. E. Gibbons, Atomic 
Energy Comm., NYO-4849, March 1957, 
124 pp. 

A Selection Study for an Advanced En- 
gineering Test Reactor, Aeronutronic 
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